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Berechnung von Beulwerten mit Hilfe von 
Mehrstellenoperatoren") 


Von Pretrer LAucHLI und Heinz WALDBURGER, Ziirich 2) 


1. Problemstellung 


Die Differentialgleichung fiir das elastische Beulen der isotropen ebenen 
Platte [1]3) 


ad 
AAu = — Dp (GM 2 Bey Bay Ey tt .) (1) 


wobei d = Plattenstarke und D = Plattensteifigkeit, reduziert sich fiir den 


_ Fall des Schub-Beulens auf die folgende Form [2]: 


yO ey (2) 


po xy xy 


Mit 


erhalten wir 


AAu—2 k= uw, 0 (3) 


(o, = Eulersche Knickspannung eines Plattenstreifens der Breite 1 und der 
Lange b). Der Beulwert k vergleicht die Beulsicherheit der Platte mit derjenigen 


eines entsprechenden Plattenstreifens. 


Zur numerischen Behandlung des durch Gleichung (3) gegebenen Eigen- 


_wertproblems (der erste Eigenwert 2 k 2?/b? entspricht dem gesuchten Beul- 


wert k, der zugehérige Eigenvektor der Beulflache) wahlen wir ein verbessertes 


_ Differenzenverfahren (Mehrstellenverfahren [3]), welches unser Problem in ein 


allgemeines Matrix-Eigenwertproblem 


(A—AB)u=0 | (4) 


uberfiihrt. 


1) Die vorliegende Arbeit wurde erméglicht durch die Schweizerische Volkswirtschaftsstiftung. 


Die Verfasser sind dieser Institution zu bestem Dank verpflichtet. 


2) Institut fiir angewandte Mathematil der ETH. 
3) Die Ziffern in eckigen Klammern verweisen auf das Literaturyerzeichnis, Seite 454. 


— 2 


446 PETER LAucHLI und HEINZ WALDBURGER ZAMP 


Die Differenzen-Operatoren OP, und OP;, welche die Matrizen A und B 
erzeugen, werden derart aus den Differentialausdriicken hergeleitet, dass in 
deren Entwicklung nach dem Taylorschen Satz u und die partiellen Ableitungen 
yon u bis zu Gliedern méglichst hoher Ordnung verschwinden. Dabei wird es 
uns gelingen, einen dem Mehrstellenverfahren entsprechenden héheren Taylor- 
abgleich zu erhalten, ohne indessen gewisse Ableitungen als neue Unbekannte 
einfiihren zu miissen. 


2. Herleitung der Operatoren 


Wir wahlen als OP, den nachstehenden von CoLLatz [3] gegebenen Aus- 
druck fiir AA4u in einem quadratischen Gitter: : 


und bestimmen eine dazu passende Approximation von #,, durch einen Diffe-— 
renzenausdruck, genauer einen Differenzenoperator OP; derart, dass . 


/ 


il 
1 1 ; 
ht 1 2-1 )u,, = (OP,) w+ (ey. (6) 
1 


: 

é 

Dadurch wird erreicht, dass beim Ers~tzen des Eigenwertproblems AAu = M thy . 
durch ein Differenzengleichungssystem die Abbrechfehler der beiden Seiten _ 
der Eigenwertgleichung von derselben Ordnung sind, namlich ~ /8. . 
Zur Herleitung von OP; driicken wir zunachst den Ausdruck 


findet 
Sy = Oty +H (ty, + Uy yay + EC...) , 
also 
1 


if it 
= hl 1 2 } Ugy = 314 thy + WS (tey + Uy slay + HE (0.2) . (7) 


durch w,, und seine partiellen Ableitungen im Sternmittelpunkt aus; man 
1 


Pa oe aa 
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Nun ist aber bekanntlich 


2 Me,» ge 
; 2 
| ° OO) = 4 Way + WE (Oley + teyy)ey + HE (.-.) 5 (8) 
ite 


so dass mit 


‘ — @ O 1 
OP, = : a : ‘ 
oe 


/ 


in der Tat die Beziehung (6) besteht. Damit wird nun 


/ 


L 
(OP,) u —A'(OP;) u = - h# ( 2 | (AAuw — 2 u,,) + AB (...) , (9) 
] 


so dass das gestellte Eigenwertproblem (3) mit 4’ = 2 k 22/b? bis auf Fehler der 
Ordnung h* durch das algebraische Eigenwertproblem 


(OP,) u = 4 (OP)) u (10) 


ersetzt werden kann, welches von der Form (4) ist. 
Es ist allerdings etwas bequemer, statt des Operators OP; den Operator 


; Sle 5 re k 
1 Q --1 
einzufiihren und das Eigenwertproblem 
(OP,) w=A(OP,) u tt LON 


zu lésen; die Beulwerte & stehen dann nach (3) mit den A in der Beziehung 


2, 2 
~ 302 he we) 
Die Operatoren OP, und OP; fiir die Matrizen A und B sind nur dann in der 
in Gleichungen (5) und (11) angegebenen Form giiltig, wenn alle beteiligten 
Punkte innerhalb des Gebietes liegen. Ihre Anpassung an die Randbedingungen 
sowie die Aufstellung der Matrizen sind anhand des nachstehenden Beispiels 
ausfiihrlich beschrieben. 


? 
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3. Beispiel 


Figur 1 
Rechteckige Platte mit Netz. 


' 


Zur Berechnung des Beulwertes k der Platte mit den Seiten a= 2 und 
b = 1 legen wir ein quadratisches Netz mit der Maschenweite h = 0,2 derart 
iiber die Platte, dass das Gitter sich auch mit dem Rand der Platte deckt. Die 
Platte sei an allen Randern frei drehbar (gelenkig) gelagert, was den folgenden 
Randbedingungen entspricht [4]: 


“=U,, =0 fin. 4 = Ound Aa 


“u=u,,=0 fir y= 0 und y =o. 


Figur 2 
Ausschnitt aus dem Netz. 


Mit den Bezeichnungen von Figur 2 kénnen wir die zweite Ableitung w,, 
schreiben als 
U, — 2 Uy + Uz 


Uy x ra h2 


(13) 


Ist wv, ein Randpunkt, also wu, = 0, so wird auch w,, = 0, da die angenommene 
Lagerung der Platte den ausserhalb des Randes liegenden fiktiven Punkten 
gleiche Betrage, aber umgekehrte Vorzeichen zuordnet (vu, = — us). Analoges 
gilt fir w,,,,. 

Unter Verwendung der «Spiegelung am Rand» stellen wir nun die modifi- 
zierten Randoperatoren auf, welche die Vorschriften geben fiir die Gleichungen 
der randnahen Punkte. Figur 3 zeigt die Modifikationen des Operators OP, 
Figur 4 diejenigen des Operators OP,. 


ee ee ee 


P 


ae TY Te 
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Figur 3 
Operator OPg mit Modifikationen. 


ah Red 35 -M 
10 |36 0 17 

Regen ih 
SS 


on: 


Figur 4 
Operator OP 4 mit Modifikationen. 


‘lenoperatoren 


449 


Dadurch, dass die Differenzengleichung in jedem inneren Punkte des 


gleichungen (4), deren absolut kleinster Eigenwert gemiiss 


wD | Vie 


fae Plant ee 


mit der kritischen Schubspannung zusammenhanet. 


_ ZAMP X1/29 


Netzes angeschrieben wird, entsteht die Gesamtheit der Matrix-Eigenwert- 


(14) 
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Die beiden Matrizen A und B sind durch OP, und OP, und das Netz 


bestimmt. 

Eine Betrachtung der Symmetrieeigenschaften, welche die beiden Differen- 
tialoperatoren im rechteckigen Grundgebiet aufweisen und welche sich auch 
auf die beiden Matrizen A und B iibertragen, zeigt, dass die Beulflachen 


facut mom Oma 10 11 12 13 14 15 16 17 78 


psa a 
M55 12 ais 


8 
Oh if 
7 7 
elf 
710 0 
0-11 -2 
i 


1 

2 e a 7 

3 | 7-1 35-11 1 2-10 Da 
4101-1340 1 -2-T if Se 
5-1-2 1035-07 O-fl-27 01 i 

6 | 2-0-2 7-10 36-10 f 2-0-2 11a tf 
7 | 1-2 -10-2 1-0 86-101 -2-10-2 0 7 1 
8) 0 7-2-10 1-03 Of 2 Ue a 
G) 1. O0-N-2 1 0 3-0 tl Oleg ff a 
NE Si es NSU Se 
11 US I ae Ie 7 -10 36 -10 7 Adis 
12 MA AY Ue) 1 -10 35 0 hath 
13 G1 J-8 O-N-2—1 6 Go —10 2 
14 he ee Od 10 ral 10 39) 
15 TE WIS, hE Te 
16 WMO SE INL TP Te ATS 
17 OE SPS lS 
18 a hl AES = 


Figur 5 


Durch den Operator OP, erzeugte Matrix A. 


punktsymmetrisch bzw. -antimetrisch sein miissen. Man kann deshalb die 
Anzahl der Variablen auf die Halfte reduzieren und den symmetrischen und 
den antimetrischen Fall getrennt behandeln (siehe Figur 1), selbstverstandlich 
unter Erhaltung der gewohnlichen Diagonalsymmetrie von A und B. Damit 
wird der Rechenaufwand erheblich herabgesetzt. Figur 5 und Figur 5a zeigen 
die Matrizen im symmetrischen Fall unseres Beispiels von Figur 1. 

Im Falle des Quadrates ware noch eine weitere Reduktion méglich gewe- 
sen; wir haben davon jedoch keinen Gebrauch gemacht. Schliesslich brachte — 
eine darstellungstheoretische Diskussion am ehesten Licht in die Symmetrie- 
betrachtungen; sie lieferte insbesondere sofort die Dimensionen der invarianten . 


Unterréume von A-1B. 
4. Behandlung des Matrix-Eigenwertproblems 


Im folgenden soll noch kurz die numerische Methode skizziert werden, nach 
welcher die Eigenwerte und -vektoren von 


. 
(A —AB)u=0 as) 
: 


auf der ERMETH berechnet wurden. 
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Es wurde die Tatsache beniitzt, dass A und B symmetrisch sind und 
dass ferner A positiv definit ist und sich somit symmetrisch in ein Produkt von 
zwei Dreiecksmatrizen zerlegen lasst: 


A=LR mit L=R; 7,50 fir i>k. 


P2324 5 607 8-9 JOU 1293 1415 16 17 18 

1 Re ee at 

2 Q001 0-1 

8 00001 0-1 

4 90000100 

5 100000 00-1 

6 Bee PE OP De 

7 OG, Sell (ROG 8 GI Oot 

8 TO0000007 08 

9 90710000000-1 af 

10 oer ey ew et 

71 hee ea Ba eas 

12 1000000071098 

13 CU? oOGe oO U GH 

It, SPO POO -Oeney 

15 had 0 OOO hel 

16 ‘10 0.0 OQ 0 7 

17 OO f= 00 8 

18 a ie ea) 

Figur 5a 
Durch den Operator OPg erzeugte Matrix B. 
Statt (15) lésen wir dann das Problem 
I-7A (A—AB)R'* Ru=0 
oder: 
(pE— C)o=0, (16) 

wobei 


pHa); “CHI BRT, = Ru. 


Damit ist die urspriingliche Aufgabe auf ein spezielles Matrix-Eigenwertpro- 
blem mit symmetrischem C zurtickgefiihrt. 

Diese neue Gleichung wird mit einem modifizierten Iterationsverfahren’) 
- angegangen, welches zuerst die absolut gréssten yw, und damit die vor allem 
interessierenden kleinsten / liefert. Es werden standig zwei (im allgemeineren 
Fall®) mehrere) Vektoren v, und v, mitgeftihrt, die wir zur zweispaltigen 
Matrix V = [v,; v,] zusammenfassen. 


4) Die Tatsache, dass B singular sein und somit (16) den Eigenwert #4 = 0 haben kann, ist hier 


_ unerheblich. 
5) Es handelt sich um den Spezialfall eines von H, RuTISHAUSFR angegebenen Verfahrens, 


welches spater publiziert werden wird. 


, 
f 
; 
: 
: 
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Der k-te Iterationsschritt besteht aus folgenden Operationen: 

V®) wird orthonormiert, und zwar orthogonal zu allfalligen schon berech- 
neten Eigenvektoren. Diese neue Matrix nennen wir Y = [¥,; y_|. Nun dreht 
man in der Ebene der y, auf Hauptachsen, das heisst Y ist mit einer solchen 
zweireihigen orthogonalen Matrix U zu multiplizieren, dass die Spalten von 
YU extremale Rayleighsche Quotienten liefern. Wegen der Orthogonalitat von 
Y heisst das aber nichts anderes, als dass in jedem Schritt das «kleine Eigen- 
wertproblem » fiir die zweireihige Matrix Y7 CY zu lésen ist. Die Eigenvektoren 
sind die Spalten von U, die Eigenwerte sind zugleich Naherungen fir die 
Eigenwerte w von (16). Die Iteration von YU mit C ergibt das neue V: 


Vet = CYU., 


Praktisch wird zuerst Z = CY und Y‘Z gebildet und, nach der Berechnung 


von U, V“+1) = ZU. Es wird bei derartigen Problemen oft iibersehen, ist aber 
numerisch von ausserordentlicher Wichtigkeit, dass die Matrix C = L+BR+ 
gar nicht explizit aufgestellt zu werden braucht, sondern es wird in jedem 
Iterationsschritt der Reihe nach von links mit R~+, B und L~! multipliziert. 
Fir die Multiplikation mit R-1+ und mit L~1, sowie auch fiir die Zerlegung 
A = LR am Anfang, wurden bei uns Teile eines Programms fiir die Auflésung 


von symmetrischen Gleichungssystemen nach GAUSS-CHOLESKI verwendet, so — 


dass sich auch die Berechnung von L~! und R-! eriibrigte. 
Als Kriterium fiir das Abbrechen der Iteration kann die absolut grésste 


Anderung in den Komponenten des normierten Vektors y, mit einer vorgege- 


benen Toleranz verglichen werden. 

Am Schluss hat man wieder auf 4 = 1/u und u = R~-1v zuriickzutransfor- 
mieren. 

Das beschriebene Verfahren liefert ohne weiteres Paare von entgegengesetzt 
gleichen Eigenwerten, wie sie in unserer Anwendung auftreten. 


5. Resultate 


Figur 6 zeigt den Beulwert k, und die zugehorige Beulflache fiir den sym- 
metrischen Fall: 


ks = 6,700 


Figur 6 
Beulflache im symmetrischen Fall. 


he eae - 
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Figur 7 zeigt den Beulwert k, und die Beulflache im antimetrischen Fall: 
Ge. oe ky == 6,744 


Figur 7 


Beulflache im antimetrischen Fall. 


Da k, < k,, wird die Platte (Seitenverhiltnis 2:1) also gemass Figur 6 aus- 
beulen. Fiir héhere Seitenverhaltnisse, z. B. 3:1, ist hingegen der dem anti- 
metrischen Fall entsprechende Beulwert &, kritisch. 

Zur Erzeugung der Matrizen A und B fiir beliebige Seitenverhaltnisse a: } 
und fiir beliebige Maschenweite 4 wurde ebenfalls ein Programm aufgestellt, 
welches erlaubt hat, die weiter unten angegebenen Falle durchzurechnen, ohne 
dass dazu die Matrizen A und B von Hand aufgestellt werden mussten. 

Die nachstehende Tabelle zeigt die Resultate der Rechnungen fiir verschie- 
dene Seitenverhaltnisse mit jeweils verschiedener Maschenweite h. 


innere Gitterpunkte h A, ky he | (ee 
Quadvat 
ee 2 af 23,238 14,127 36,000 21,885 
4x 4 Jf 5,807 9,806 7,324 12,369 
ex 6 Gite 2,855 9,451 3,549 11,747 
S56 48 ats 713 9,372 223, 11,618 
Rechteck 2:7 
205 1/, 13,236 8,046 13,555 8,240 
4x 9 he 3,968 6,700 _ 3,994 6,744 
Rechteck 3:17 | 
A 8 Ne 11,848 7,203 11,518 7,002 
4x14 ee Re fi | 3,600 6,079 | 


6. Extrapolation der Beulwerte 


Da aus der Gleichung (9) durch Division mit h* hervorgeht, dass die Diffe- 
renzengleichung (10) bis auf einen Fehler der Ordnung /* mit AAu =i" u,, 
‘iibereinstimmt, kann man annehmen, dass die Eigenwerte von (10) und damit 
‘die Beulwerte (12) mit Fehlern behaftet sind, die im wesentlichen zu #4 propor- 
‘tional sind. Héhere Fehlerglieder sind zu h, h8, usw., proportional. 
Man kann daher eine bessere Approximation fiir den Beulwert finden, 
indem man mit Hilfe der Werte &(#), die man mit verschiedenen Maschen- 


AL aN loa 
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weiten / gefunden hat, mit dem Ansatz 
k(h) =k+cht+daho+-: 


das Interpolationspolynom konstruiert und k(0) berechnet. Beispielsweise 
wurden fiir den kleinsten antimetrischen Beulwert die in der folgenden Tabelle 


angefiihrten Werte k(h) gefunden. 


| Netz | h | k(h) | (k(h) — k)/h* 
2x2 0,333333 21,885 837 
4x4 0,2 12,369 514 
6x6 0,142857 11,747 480 
SIP as CO) aati ital 11,618 466 


Mit dem Ansatz k(h) = k + c h* + d h® findet man aus den drei letzten Tabel- 
lenwerten: k= 11,547: 

und damit die neben den k(h) aufgefiithrten Werte von (k(h) — k)/h4, welche of- 
fenbar den Ansatz rechtfertigen. 

Es soll freilich nicht verschwiegen werden, dass diese Art der Extrapolation 
nicht ohne Gefahren ist. Sie soll jedenfalls nur unternommen werden, wenn die 
Beulwerte fiir mehrere verschiedene Maschenweiten berechnet wurden; ausser- 
dem sollen immer auch die Kontrollgréssen (k(2) — k)/h* berechnet werden; 
bei unruhigem Verlauf derselben ist die Extrapolation zu verwerfen. . 

Es ist schliesslich noch darauf hinzuweisen, dass diese Extrapolation nur 
dann mit Erfolg ausgefiihrt werden kann, wenn die Beulwerte fiir die verschie- 
denen Maschenweiten / mit iiberschiissiger Genauigkeit berechnet wurden. Es 
wirkt sich daher sehr vorteilhaft aus, dass Rechenautomaten mit wesentlich 


mehr Stellen arbeiten, als man bei technischen Problemen normalerweise 
bendotigt. i 
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Summary 


Buckling of a plate leads to an eigenvalue problem for which a very accurat 
finite difference approximation is given. The resulting algebraic eigenvalue prob- 
lem Ax = 1 Bx is solved for specific problems by an improved iteration method 
Numerical results are discussed. 


(Eingegangen: 14. Marz 1960.) 
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Kritische Drehzahlen unter Torsion und Druck 
bei Beriicksichtigung des Eigengewichtes 


Von Horst Lerpuorz, Stuttgart, Deutschland) 


1. Einleitung 


Es sollen einfach besetzte Wellen mit mehrfach beaufschlagten Turbinen- 
scheiben betrachtet werden. Die Wellen seien zusatzlich durch ¥8 auf Torsion 
und durch $$ und q auf Druck beansprucht. 

Das Torsionsmoment Y werde durch die Turbinenscheibe in die Welle 
eingeleitet; 8 greife als Einzelkraft an den Enden der Welle an; q sei eine 
kontinuierlich verteilte Langskraft (zum Beispiel Eigengewicht); $$} und q 
mégen stets die Richtung der unverformten Wellenachse beibehalten; YB sei 
semitangential. Die Welle sei von gleichbleibendem kreisf6rmigem Querschnitt ; 
die Biegesteifigkeit a und die Torsionssteifigkeit c sind daher konstante Gréssen. 
Alle Verformungen, die die Welle erleidet, seien sehr klein. 

Bekanntlich [1]?) werden die kritischen Drehzahlen durch die Torsions- und 

-Druckbeanspruchung der Welle beeinflusst. Untersuchungen hiertiber sind auch 
von H. Z1EGLER [2] durchgefiihrt worden, der insbesondere gezeigt hat, dass das 
Problem fiir die angenommene Belastung konservativ ist und dass die Annahme 
eines semitangentialen Momentes fiir mehrfach beaufschlagte Scheiben be- 
sonders sinnvoll ist. 

Unter den gemachten Voraussetzungen und ohne Beriicksichtigung der | 
Kreiselwirkung gilt, wie man bei H. ZIEGLER [3] nachliest, die Beziehung 

aoe (1) 

wobei @, die kritische Drehzahl, m die Masse der Scheibe und « die Einflusszahl 

fiir die Durchbiegung der Welle an der Stelle ist, wo die Turbinenscheibe sitzt. 
Das ganze Problem besteht also in der Bestimmung von OL. 

: Es sei noch besonders darauf hingewiesen, dass sich die kritische Drehzahl 

‘nur deshalb in der einfachen Beziehung (1) darstellt, weil das Torsionsmoment 

als semitangential vorausgesetzt worden ist. Die Begriindung hierfiir ist in der 
zitierten Arbeit von H. ZIEGLER angegeben. 


1) Institut fiir technische Mechanik, Technische Hochschule Stuttgart. fo 
2) Die Ziffern in den eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 471. 
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Als neu gegeniiber friiheren Arbeiten soll jetzt « auch noch unter Beriick- 
sichtigung von q errechnet werden. 


2. Differentialgleichungen und Randbedingungen 


An Hand von zwei Beispielen soll die Aufstellung der Differentialgleichun- 
gen und der Randbedingungen gezeigt werden. Die Uberlegungen lassen sich 
dann leicht auf andere Falle tibertragen. 


2.1 Die flregende Welle 


Die fliegende Welle ist in Abbildung 1 dargestellt. Sie ist mit dem unteren 
Ende in einem «langen» Lager eingespannt und erfahrt am oberen Ende durch 
eine Horizontalkraft § eine Auslenkung. 


W 
z 
7 |} 
= H' 1 y 
pW) 4 * te 
Y p 
\f 
q(2) 
uf 
a b 
Abbildung 1 
Die fliegende Welle. 
Die Belastungen haben folgende Komponentendarstellung: 
newt» : 
TD Wh (a aee ailialale . 
pa Ce! eae 
2 
g = ( 0, OF Y ) , 
pr ear 


‘S 
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und der Tangenten-Einheitsvektor der Stabachse ist 
t= (*', y= &,, 44). (3) 


Fiir kleine Auslenkungen der Wellenachse hat man fiir die Komponenten 
des inneren Stabmomentes die Beziehungen 


M,= ax", | 
M,=-—ay", (4) 
W =cr=konst., | 


wobei t die infolge der Torsion auftretende Verwindung pro Langeneinheit der 
Welle ist. 

Durch Anwendung des Schnittprinzips und der statischen Gleichgewichts- 
bedingungen erhalt man 


ax" =—Wy'+ = Wey) At = 2) EP xd) = 2] 


l 


+ f a) (e) — xl at, 


ee 
ay" =Wx'— Wx) + Al—2) + PH) --yI 
+ | 9 (©) — y@) ae. 


Es seien noch weitere Bezeichnungen eingefiihrt: 


: ii 
P+ [ 4 (¢) a 
— 


w= ha, b=, 8h) = wat 6) 
a - a : a 
und es werde (3) verwendet. Dann erhalt man aus (5) durch Differentiation : 
fH +wt + ge)t=—*,, | : 
(7) 
i" — wt, + g(z) ty = — hy. | 


Das ist ein lineares, inhomogenes Differentialgleichungssystem 2. Ordnung 
mit teils konstanten und dem variablen Koeffizienten g(z). Abhangige Variable 
sind die Komponenten #, und #, des Tangenteneinheitsvektors, unabhangige 
Variable ist z, welches die Lange der Wellenachse vom unteren Ende her zahlt. 
s 


PY ee OR 
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Als Randbedingungen hat man 
2.11: wegen der Einspannung des unteren Wellenendes: 


P=), ht) use == (0), (8) 
2.12: als Momentenbedingung am oberen Wellenende: 
r if 
a t. a 2 WwW t, 
fir z=T1. (9) 
at’ = - Wt | 
y 2 x 


2.2 Die beiderseits gelagerte Welle 


Wie Abbildung 2 zeigt, soll die Welle beiderseits «lange» Lager haben, SO. 
dass ihre beiden Enden als eingespannt gelten kénnen. 


a) . att 


eq Pe ae ee 


2 oo pee 


sh 


Abbildung 2 
Die beiderseits gelagerte Welle. 


Fiir die an der Welle angreifenden Belastungen gilt wieder (2), fiir die Kom- 
ponenten des inneren Momentes (4). Zusatzlich braucht man jetzt fiir di 
Rechnung noch die an den Lagern auftretenden Auflagerkrafte. Diese habe 
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die Komponenten: 


t— 8 W , , 
A, at A, asi sii 2 1 [y (0) en (z,)] , 


i — 2 W , , 7 
A,= A, l ai Zl ie oe (0) 

as ° W i iM } 
3 A= eg [y'(z,) — y'(0)], 


Zs Ve es hiehe 
B, = ff. i “ >] [x’(0) — x'(z,)] , 


— 


und es ist z, der Wert von z, der angibt, wo die Turbinenscheibe sitzt. 
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(10) 


Wieder durch Anwendung des Schnittprinzips und der statischen Gleich- 


gewichtsbedingungen liest man aus Abbildung 2 ab 


2.21: fiir den Bereich 0 < z <z,: 


x 


I 
ax"=—Wy' + 5 Wy'(0) A,z—Px x | qd: 
0 


z 


+ f a6) be) — meat, 


. (11) 
ay =W : y= Wa 0) 403 =P ay / g az 
+ [ a6) ve) — yea, 
; é 
2.22: fiir den Bereich ey eB 
| 
: ax" = — B,(-2)— Px— fg) ee) — xO, | 
; . 
; - | (12) 
| ay"=— B,(—2)-—Py- rs q(¢) [v(z) — y(¢)] 4 
, Verwendet man auch hier (3), (6), sowie 
A B, 
Eg Se Se me (13) 


hungssysteme 


nd differentiiert man (11) und (12), so hat man die beiden Differentialglei- 
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2.21: fiir den Bereich 0 =z =z,: 


ti + wt + o(2)t, = — a, 
i" — wil +e) t=—4y | 
und 2.22: fiir den Bereich z, = 2 S27: 
ti + e(z) t= B,, 
gle) t, = B a 
i" + g(2) ty = By. | 
In diesem Fall gelten folgende Randbedingungen: 
2.23: wegen der Einspannung der Wellenenden: 
i=0, =0 fir z<=0, z=l, (16) 
2.24: aus Stetigkeitsgriinden : 
l. (z, — 0) = ty (z, at 0) , | 
(17) 
ty (% — 0) = ty (2, + 0), | 
2.25: wegen der Unverschieblichkeit des oberen Wellenendes: 
| t,(z)dz+ | U(2z)dz=0, (j=%,4). (18) 


0 Bg 


3. Die allgemeine Lésung der Differentialgleichungen 


Von den bisher aufgestellten Differentialgleichungen sind (7) und (14) vom 


gleichen Typ. Wie man ihre Lésungen erhilt, soll am Beispiel des Systems (7) 
gezeigt werden: ‘ 


Man fiihrt die komplexen Gréssen 


o=t,+2 a | 

(19) 

h=h,+th, | : 

ein und erhalt statt (7): ; 
oe" —iwo’+go=—h. i (20) 

Als neue Bezeichnungen seien ; 
u=—* 

=— 5, : 

(21) 

oo we ; 

Vion cea | 
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benutzt, dann wird aus (20): 


e" +2iue—(w+V)o=—h, (22) 
was sich mit der Transformation 
O =F ertu= (23) 
auf 
r" —Vr=—heirz (24) 


zurtickfiihren lasst. Die zu (24) gehérende homogene Gleichung 
ee Vel) (25) 


sei die «Basisgleichung» genannt, denn ihre linear unabhangigen Fundamental- 
l6sungen 7, und 7, beherrschen die weitere Rechnung. Der Charakter von 7, 
und 7, ist bestimmt, wegen V und der Bedeutung von g, durch die Art der 
Langsbelastung q(z). 

Unter Verwendung der Wronskischen Determinante 


: (26) 


die in diesem Fall eine Konstante ist, erhalt man durch Variation der komplexen 


_ Integrationskonstanten Y% und 8 fiir die inhomogene Differentialgleichung (24) 


die allgemeine Lésung 


r=Ur,+ Bret ae ( e: 4, 6°"? dz — r | 1, eiue 7 : (27) 
0 0 


| Hieraus folgt, wegen (23), als allgemeine Lésung der Differentialgleichung (22):. 


o= Jatin (rf cue da=, | ey 3) ca oie tl Pet 
0 


0 


| Das ist dann aber auch die gesuchte Lésung des Systems (7), wenn man noch 


die durch (19) und (21) festgesetzte Bedeutung von g, # und w beachtet. 


Ganz entsprechend erhalt man mit 


a= 4, +14, | (29) 


als allgemeine Lésung von (14): 


ee SERS en se 


= lens ont (n[nere—n]n a) Crit ean (30) 
0 0 
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Als letzte der bisher betrachteten Differentialgleichungen kénnen die 
Gleichungen des Systems (15) mit 9 und mit 


auf die Form 

@ ele =p ee) 
gebracht werden, die vom Typ der inhomogenen Basisgleichung (24) ist. Be- 
zeichnet man die Fundamentallésungen von 


7" + g(2)r=0 (33) 
mit 7f und 7§, und die zu (33) gehérende Wronskische Determinante mit W%, 


so bekommt man daher, ganz entsprechend wie bei (24), fiir (32), durch Varia- 
tion der Konstanten, die allgemeine Lésung 


eG Beit gpe(—ni [nee | ds) (34) 
20 0 


3.1 Die Berechnung von a 


Wie (1) zeigt, kennt man die gesuchte kritische Drehzahl w,, bei gegebener 


Scheibenmasse m, vollkommen, wenn es gelungen ist, die Einflusszahl « fiir die — 


Durchbiegung der Welle an der Stelle z = z, anzugeben. 
Der Bedeutung von @ entsprechend gibt 


v(z,) = fe dz (35) 


gerade die Durchbiegung der Welle in vektorieller Form unter den in den Ab- 
bildungen 1 und 2 gezeigten Belastungen. 


Wie schon H. ZIEGLER [4] gezeigt hat und wie es auch hier die weitere 


Rechnung bestatigen wird, ist es dem Wirken eines semitangentialen Torsions- — 


momentes zu verdanken, dass die Welle sich unter einer beliebig gerichteten 
horizontalen Kraft § in Richtung der Kraft § verschiebt. Das heisst aber: 


Auslenkung v(z,) und Kraft § haben gleiche Wirkungslinie. Wahlt man also — 
fiir die Komponenten von § Krafte vom Betrag 1, so wird auch v(z,) Kompo- 
nenten haben, die betragsmissig gleich sind, und dieser gemeinsame Betrag — 


ist das gesuchte «. 
Die weitere Rechnung hat demnach so zu verlaufen, dass man 


v(z,) = A o dz 


unter der besonderen Bedingung § = 1 + i (was sich in entsprechenden Werten 
fiir 4, « und # ausdriickt) errechnet und den Realteil (oder den gleichgrossen 


Vol. XI, 1960 Kritische Drehzahlen unter Torsion und Druck 463 


Imaginarteil) von v(z,) ermittelt. Dann ist 
a = Re(n(z,)) fir §=1+;. (36) 


Fiir 9 stehen die Lésungen (28) oder (30) zur Verfiigung 


l 
(wav man mit (34) rechnen, so ware v(z,) = — i: 0 7 } 
in denen aber noch Integrationskonstanten vorkommen. Die Anwendung der 
in Absatz 2 gegebenen Randbedingungen gestattet jedoch die Bestimmung 


dieser Konstanten. Es ist nur noch nétig, auch die Randbedin 


gungen 
komplex zu schreiben. Aus (8) wird 


o=0 fir z=0; (37) 
aus (9) wird 
ae a 
o=>twe fir z=]; (38) 
aus (16) wird 
o=0 fir z=0 und z=1]; (39) 
aus (17) wird 
@ (2, — 0) = @ (2, + 0); (40) 
und aus (18) wird 
25 l 
/ o(2) dz a o(z) dz=0. (41) 
0 25 


Hat man mit Hilfe der Randbedingungen (37), (38) oder (39) bis (41) die 
Integrationskonstanten ermittelt, so steht der Berechnung von p(z,) und dar- 
iiber hinaus von « nichts mehr im Wege. Die wirkliche Durchfiihrung der bisher 
nur angedeuteten Rechnung soll am Beispiel der fliegenden Welle erfolgen. 


4. Kritische Drehzahl fiir die fliegende Welle 
Fiir diesen Fall muss man (28), (37) und (38) benutzen. Durch Anwendung 


- der Randbedingungen (37) und (38) auf die allgemeine Lésung (28) von 9 er- 


Se ee 


halt man fiir die Integrationskonstanten: 
l l 


7; (L) A ¥, 6% dz — 1,(1) / 7, ett® dz 
2,0 


Ww ¥4(0) 72(0) — 730) 740) ! (42) 
l I) 
ri) [ ree dz =n) fre ae 
0 h é 0 
B= 7,(0) W 73 (2) %2(0) — 75 (2) 74(0) 
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Mit den Abkiirzungen 


72(0) 74(2) Ft 72(0) 73(2) 
aE oa ove , 2 0 
w(;) w(;) its 
OLR OL) 
6 ae een /(0; 5) — 0 
w(;) w(,) 
und mit (42) folgt dann aus (28): 
I 1 
a= ee [+ freee ‘ Vp C8? CER Gat eae if (eae 
0 0 
i I 
+ dy, 1%, er Ve C "A AZO aia e a ¢ 1, ae (44) 
0 0 


Zz z 
+4 rsa bi eiuz dz — ~, e~i¥Z [r owas 
L 2 2 ue 2 
0 0 


womit man v(z,) nach (35) errechnen kann. ‘ 
Vor Ausfiithrung der nach (35) vorgeschriebenen Integration sollen weitere — 
Abkiirzungen eingefiihrt werden: 
1 l 
See [ rsinu edz, Si = f rsinuzdz, 
0 0 
l l 


10 
Sua [ rycosms de, Si = fr, cosuzdz. 


0 0 
Hiermit erhalt man: 
1 l | 
[netrds [ rods = (Se SY + Si 51) +2 (52 SB — 58 Sm, 


0 0 
U I 

A ettede | eine dz = (S12 4 (Sioa, 

é é (46) 
1 1 


[or eins de | ry e'"= dz = (SO)® + (S10)2 
0 


0 
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l 1 


a 


—iuz iuz a 10 cl0 10 , 
[ raertes dz [7,02 de = (Sit Sig + Si Sd) — i (Si 5 — Si im, 
0 


0 


l 


1 z s 
ic ( Takeda ova dz — ie e~iuz / ty oma] dz = 
ra 
J : 


/ . 
0 0 0 
I 1 l 


ee S20 RE: . p 
= / r, cosuz SS}. dz + / r, sinu z S2 dz — | r, cosu z Si dz 
0 0 


at 


0 


+65 
| (46) 
1 
= / rz sinu z SP dz + ¢ (Si° So — Si0 si) . 
0 
Benutzt man ausserdem die Bezeichnungen 

Lyo(l) = Sy, Soo + Sys Sop, Ly(l) = (St2)® + (Su)? 
LBM) = Ste Sos — Sts Sse» La(2) = (Sze)? + (S25)? 
1 l 1 


L(!) =| r,cosuz Sxdz+ | rzsinuz S2dz— | 7, cosu z S3 dz (47) 
0 0 0 
el. 
eal r, sinu z S32 dz, 
0 


so lasst sich mit (44) bei Verwendung von (45 bis 47), die Ausrechnung von 
(35) schreiben: 


v(z,) = ae [Lyo(l) +4 LEO] + ae Ly) + dy L2() 


| (48) 
+ Dag [Lyo(2) — 7 LE] — LO +2 LED}. 
Es ist aber immer, wie man an Hand von (43) nachpriifen kann, 
ay, —0.+1=0, (49) 


so dass bei (48) wegen (49) die imaginaren Terme aus der geschweiften Klam- 
mer herausfallen und man 


(z,) = {Lyall) (ara + baal + ae £1) + bun Zell) - LO} (50) 


hat. Daraus ersieht man, dass das Wellenende sich tatsiachlich unter einer 
beliebig gerichteten horizontalen Einheitskraft um die Strecke a in Richtung 
der Kraft verschiebt. 


ZAMP X1/30 
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Aus (50) erhalt man die gesuchte Einflusszahl «, indem man fiir h = 1/a 
setzt. Das ergibt: 


& =F {Lrall) (aaa + Pred + a2 Lill) + Or Lal) — LO}. G1) 


Die Gleichung (51) kann fiir alle méglichen Lastverteilungsgesetze (2) 
angewendet werden. 

Der Fall g = 0 wurde bereits von H. Z1EGLER behandelt. Hier besteht die 
Druckbelastung nur aus der Einzelkraft P, die an den Enden der Welle an- 
greift. Mit P TP 


2 
b=—-Vagts —- + aay = konst. (52) 


heissen daher die Fundamentallésungen der Basisgleichung (25): 


7, = SiN Uy Z | 
(53) 
15 = COSU, 2 - | 
Das gibt: 
7,(0)=0, 7(0)=1, ri (2) =vecosng?, 75 (1) =— vpsindgl , 
Ww 0 
='—= 0); ) W fe thes Vp COS Ug L , 
so dass man nach (43) 
My=—1, a.=—tgrl, by=0, db, =—0 (54) 
erhalt. 
Setzt man (54) in (51) ein, so ist 
= 1 
a = — —-{— Lyall) — tg vp £4) — LQ}. (55) 


Verwendet man fiir die Ausrechnung von L(2), L,(1) und L,,(/) die Formeln 
(45-47) und fiir 7, und 7, (53), so ergibt sich: 


Lf) = > [(2 + =") sinv,/ cosu 1] — eee cosv,/sinul — vo i | 


a 2u : é 
EQ) = > [2 ao > — sin?v, 1 — 2 a COS, 1 cosu L 


Avgu .. Pe) 
—2 mae SiN Vp / sin uw i ; 
ll) = a [Sin Vg / cosul — sinv, 1 cos) /] . 
Dies, in (55) eingesetzt, fiihrt nach einigen Zwischenrechnungen zu 
- 1 fuet+uwa 2ua sinul 
sik pled an ov ? sinu 
iP | ee 23 Ae P cosayt 22 
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Der Ausdruck (57) stimmt aber genau mit einer fiir diesen Lastfall auf anderem 
Wege von H. ZIEGLER [5] berechneten Formel iiberein, wenn man den Unter- 
schied in der Bezeichnungsweise beseitigt : 


den Bezeichnungen 1 ; 
g ; in dieser 
bei H. ZIEGLER entsprechen Arbeit 
w —_—_______— —2u 
oO —_____— 205 
b a Pla. 


Der Fall g(z) =q=konst. soll jetzt neu in diese Arbeit aufgenommen 
werden. 


Als Fundamentallésungen der Basisgleichung (25) hat man fiir ihn: 


2 3 
ry = 3"? Sirs (= ae} ; 


e (58) 
1%. = 3°? J ais (} 3%) ; 
mit 
3 = a B23 — Bus z , | 
- - (59) 
Sod 1 eee ig ( 
2 a sa a FS 4a’ ge yes | 


ip 
Ay, Ay9, O44, O49, L(l), Ly (1), L.(2) und L,,(/) zu verwenden. Dann hat man alles 
fiir die Bestimmung von « nach (51) zur Verfiigung. 

Die Rechnung muss im konkreten Fall numerisch erfolgen, da sich wegen 
der Kompliziertheit der Funktionen die nach (45) und (46) vorgeschriebenen 
Integrale wohl kaum geschlossen ausrechnen lassen. Die Durchfiithrung der 
Rechnung wird zur Erlauterung der Formeln an einem Zahlenbeispiel vorge- 
nommen. Fiir die Bestimmung der Funktionen /J,,,; und /_,,3 konnten die 
Tafeln von Watson [6] benutzt werden, da die Argumente fiir das gewahlte 
Beispiel alle positiv reell bleiben. 

Gegeben seien: 


; : ey 0 
Die Funktionen (58) sind in (26, 43, 45-47) zur Berechnung von W, w( 


g=10*kpjm, T=2-10"%kpm, a=10'kpm?, /=3m. 


Es ist die Konstante W wegen 


A eee ae a Me NY 


a, 2 3/2 
—* = B83 Jos (= 3"), 


dy 2 43/2 
@ hice — Pe 3 Jas (53 ) 
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gleich Se ie 
W = B28 33? (Jus Jas + Jas J 233) = BYP zeae (61) 
Mit den gegebenen Werten errechnet man nach (59) und (21): 
3-10? 40-108 ~ fp A0t gees 
Soa Oe qo 190 13 Oe gee 
z= 13+ 10-* (10-*)-28 — (10-%) #8 z = 0,215 (13 — 2) , 
2-109? 
a aoe ae 0,316. 


Die weitere Rechnung wurde in Tabelle I durchgefiihrt. Mit den Werten 
der Tabelle I erhalt man fiir 7 = 3 durch numerische Differentiation: 


71 (3) = ve (0,095 — 8 - 0,352 + 8 - 0,738 — 0,835) = 0,196 , 


r,(3) =, (— 0,612 + 8 - 0,428 + 8 - 0,056 — 0,307) = 0,246 . 


Tabelle I 
[el ese 
3 2,795 2,580 2,305 2,150 1,935 1,720 1,505 
alah 1,673 1,606 1,539 1,466 1,393 nissibil 1,229 
Bale 4,657 4,173 3,652 3,151 2,703 2,253 1,856 
08 = INFO oval 2,78 2,44 PRAM) = 1,80 1350 2A 
J113(0) — 0,094 | +0,059 | +0,229 | +0,396 | +0,530 | +0,637 | +0,701 
J—113(0) — 0,429 | —0,381 ; —0,278 | —0,128 | +0,040 | +0,234 | +0,421 


1,(2) = 318+ Jx13(6) | —0,160 | +0,095 | +0,352 | +0,581 | +0,738 | +0,835 | +0,862 
r,(2) =3!2- J_1i3(c) | — 0,718 | =0,612 | —0,428 | —0,188 | +-0,056 | 40,307 | 10,517 
i ee OM Cee es i |. 


Damit und mit 7,(3) und 7,(3) aus Tabelle I erhalt man nach (43) 


0 ‘ 
w( = — 0,160 - 0,246 + 0,718 - 0,196 = 0,101 , 


i 
hae aE See 
nee EE — 1,749, 
by = ee 08155 
bie a oe — 0,390 


ne 


Pe ee a eT 


+ 
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wo 


und nach (61) 
3 
Vs = 9,245 + 0,866 = 0,178 . 


Die Berechnung der in (47) festgelegten Werte erfolgt durch numerische 
Integration in Tabelle II, aus der man 


St = 0464, SI —0,474,- "SP o — 1,340) SP =0,513; 


u l 


i %, COSU z Sy, dz = — 0,057, / 7, sinu zS,,dz = — 0,074, 
0 


0 


I U 


/ 7, cosu z S,.dz = — 0,601, af r, sinuz S,. dz = — 0,189 
d ny 
herausliest, womit man eben nach (47) 
L,,(3) = — 1,340 - 0,464 — 0,474 - 0,513 = — 0,865 , 


3 0,4642 + 0,4742 = 0,440, 


L, (3) = 
L, (3)= 1,340? + 0,513? = 2,059 , 
L (3)= 


tS 0,601 + 0,189 — 0,057 — 0,074 = 0,659 
erhalt. 
Jetzt steht alles zur Verfiigung, um nach (51) errechnen zu kénnen: 


- 1 
% = 9773 10F [(— 1,393 + 0,390) - 0,865 + 1,749 - 0,440 + 0,315 - 2,059 


— 0,659] = 1,293 - 10-4. 


Zum Vergleich sei noch nach Formel (41) gerechnet und dazu angenommen, 
dass die Langskraft nicht verteilt sei, sondern als Einzelkraft P = 3 - 10° kp 
an den Enden der Welle angreift. Alle anderen Daten sollen unverandert 
_ bleiben. 

Dann hat man nach (52): 
Le 40 + 108 
oa ioe 4-101 


—13-10-?, = 3,606-10-, 


und es ist bereits bekannt: « = — 0,316, u? = 10: 10-*. Hiermit gibt (57): 


ks 1 13-10-2 + 10-107? LO Foe 

gees se | 3,606 - 10-2 Seah koe ) 
2+ 0,316 pe 0016 e524 613.10 
3-103 ~ ee (3,606 - 10-1 - 3) 3] se : 


Horst LEIPHOLZ 
Tabelle II 
Zz (a) 1 2 | 3 
Uz 0 — 0,316 — 0,632 — 0,948 
COSU Z 1 0,950 0,807 0,583 
sinu Zz 0) — 0,301 — 0,591 — 0,812 
¥, COSU Z — 0,160 0,090 0,284 0,339 
7, SinuU z 0) — 0,029 — 0,208 — 0,472 
Yq COSU Z — 0,718 — 0,581 — 0,345 — 0,110 
Yq SiNU Z 0 0,184 0,253 0,153 
see al r, cosuz dz 0 ~ 0,035 0,152 0,464 
0 
sf = / yy sinu 2d 0 -0,015 | -0,134 | -0,474 
0 
sa = i Vy cosu z dz 0 — 0,649 =f — 1,340 
0 
a 
Sy =) 7, sinu 2 dz 0 0,092 0,310 0,513 
0 
Sip %_ COSU Z 0 0,020 ~ 0,052 — 0,051 
20 ; 
Si, % Sinu z 0 — 0,003 — 0,034 — 0,073 
Soe" 7 COSU Z 0 ~0,058 — 0,316 ~ 0,454 
S® 9, sinus 0 — 0,003 — 0,064 0,242 
x 
i Oe %> COSU z dz 0 0,010 — 0,006 — 0,057 
0 x 
Si, %2 sinu z dz 0 — 0,002 ~ 0,020 — 0,074 


20 
S5.% Cosu z dz 


20 . 
So,% Sinu z dz 
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Nach (1) ist das Verhialtnis der zugehérigen Drehzahlen («,, fiir die konti- 
nuierliche Langskraft ¢, @,p fiir dic entsprechende Einzelkraft P): 


Org 7/1913 _ 1.217 
py p 1,293 R ‘ 
woraus man ersieht, dass man bei Beriicksichtigung der kontinuierlichen Ver- 


teilung der Last eine merkliche Erhéhung der kritischen Drehzahl erhalten 
kann. 
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Summary 


The paper deals with the critical angular velocities of continuously loaded shafts. 

A method is given for the solution of the problem and an example is presented 

which shows that the critical angular velocities increase, when, instead of a 

single concentrated load on the ends, the shaft carries a continuous load of the 
same magnitude. 


(Eingegangen: 16. Marz 1960.) 


Steady-State Thermal Stresses in an Elastic Cone’) 


By Roxuro Mux! and Ett STERNBERG, Providence, R. I., U.S.A.) 


1. Introduction 


The stress analysis of circular cones, and of bodies bounded by two coaxial 
conical surfaces with a common vertex, has received repeated attention in the 
linear theory of elasticity. Thus MIcHELL [2]*) established the solution appro- 
priate to a semi-infinite solid cone under a concentrated force of arbitrary 
orientation, applied to the vertex. L. FGpp [3] disposed of the analogous pure 


1) The results communicated in this paper were obtained in the course of an investigation 
conducted under Contracts Nonr 562(20) and Nonr 562(25) of Brown University with the Office 
; of Naval Research in Washington, D. C. 
2) Division of Applied Mathematics, Brown University. 
3) Numbers in brackets refer to References, page 490. 
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torsion problem. NEUBER [4] extended MIcHELL’s results to the case of a conical 
shell of linearly flaring thickness and, in addition, treated the flexure of such 
4 shell due to a concentrated couple at the apex. Some of the solutions contained 
in [4] were recently re-derived in a different manner by Knops [5]. 

Certain axisymmetric problems concerning solid and hollow cones subjected 
to surface tractions which are distributed over the conical boundary were 
considered by Morcan [6, 7]. In [6] the lateral loading is assumed proportional 
to an arbitrary real power of the distance from the vertex, while [7] presup- 
poses piecewise continuous surface tractions which vanish identically beyond 
a fixed distance from the vertex. 

In the present thermoelastic investigation we seek the thermal stresses 
produced in a solid, homogeneous and isotropic, elastic cone of semi-infinite 
extent, which is free from external loads but is exposed to a stationary discon- 
tinuous surface temperature: the tip of the cone (up to a given distance from 
the apex) is held at a constant temperature, while the remainder of the boun- 
dary is maintained at another uniform temperature. 

Throughout the interior of the body the temperature field is required to 
obey the equation of heat-conduction. In the absence of thermoelastic coupling, 
the stress inducing temperature field is therefore characterized as the solution 
of an independent Dirichlet problem. Since the solution to this heat-conduction 
problem was deduccd separately in a previous publication we may proceed 
directly to the determination of the associated thermal stresses. 

Evidently, the thermal-stress problem to be considered presently bears a 
close mathematical relation to the ordinary elastostatic problems discussed by 
MorGAN [7]. The method of attack sketched in [7] is a generalization of a pro- 
cedure employed by TRANTER [8] in connection with similar two-dimensional 
equilibrium problems for an elastic wedge. MoRGAN’s scheme rests on the use 
of the Mellin transform together with Love’s stress function for torsionless 
rotational symmetry. In contrast to TRANTER’s work, MorGAN’s treatment is 
purely formal. He leaves the problem, except in the degenerate case of the 
half-space, after recording complex inversion integrals for the desired compo- 
nents of stress (and without specifying the appropriate choice of the corre- 
sponding path of integration). The chief complications inherent in the analysis, 
however, merely begin to arise at this stage; they stem from the need to secure 
a representation of the solution which is amenable to a complete numerical 
evaluation. 

Instead of applying Love’s stress function, we approach the thermoelastic 
problem at hand on the basis of the axisymmetric version of PAPKOVICH’S 
general solution to the elastostatic field equations, which is readily modified to 
account for the presence of a temperature field*). With the aid of the Mellin. 

7 


*) See [9], 


where this scheme was used in the stead -state t i i 
ee y-state thermoelastic study of an elastic 


. 
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transform and of the corresponding inversion theorem we thus deduce complex 
integral representations for the required temperature stresses. These complex 
integrals, in turn, are transformed into real integrals involving Legendre func- 
tions of the first kind of complex degree. The validity of the solution in this 
final form is examined; its numerical evaluation is discussed in detail and is 
illustrated by a specific example. The method adopted here lends itself equally 
well to a complete isothermal stress analysis of an elastic cone which is sub- 
jected to piecewise continuous surface loads. In fact, the latter problem is 
distinctly less troublesome than the one under present consideration. 


2. Formulation of Problem. Stress Functions 


We recall at this place the general formulation of boundary-value problems 
in the classical theory of steady-state thermoelasticity. With reference to 
rectangular cartesian coordinates x; (¢ = 1, 2,3), and in the usual indicial 
notation®), the governing system of field equations may be taken in the 
following form: 


T,,=0, (1) 
il 2(1+ 7) a 
0S ae ep eked aay ne ea (2) 
204; 
Tip = MY Mag + Met TSS DP tee — (+9) 0 7}. (3) 


Here T(x, %2, Xs), U;(*1, X2, ¥) and T;;(%1, %2, ¥3) denote the temperature field, 

the components of the displacement vector, and the components of the stress 

tensor, respectively. The constants «, w, and y, in this order, stand for the 

coefficient of thermal expansion, the shear modulus, and Polsson’s ratio, 
whereas 6;; designates KRONECKER’s delta. Equation (1) is the equation of 
heat-conduction; (2) and (3) represent the displacement-equations of equili- 

brium (in the absence of body forces) and the stress-displacement relations. 

The foregoing field equations, which must hold throughout the space region 

_D occupied by the medium, are to be supplemented by appropriate boundary 

conditions. If the temperature itself is prescribed along the boundary B, and 
_ provided the surface tractions vanish identically, one has on B 


Dee Rigs ie 8s) (4) 
t,4=0, (5) 


“where T* is the given distribution of the surface temperature and J; denotes 
the components of the unit normal to the boundary. 


: 5) Summation over repeated subscipts is implied. Subscripts preceded by a comma indicate 
partial differentiation with respect to the corresponding coordinate. 


Sees NYS 
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Once the independent Dirichlet problem characterized by (1), (4) has been 
solved, the temperature field T(x, %2, %3) in (2), (3) is known. The integration 
of (2), (3), subject to (4), may then be approached conveniently with the aid 
of an obvious modification of PaPKOVICH’s [10] general solution to the iso- 
thermal elastostatic field equations which was introduced previously in [9]. 
Thus, the complete integral of (2) admits the representation 


26 u,=9; — (3-47) vi + 4) Yi» (6) 
2(1+»v)pma 
ap ae ae Lf Oe (7) 


Consequently, the determination of the thermal stresses (and displacements) is 
reduced to finding stress functions gy, y; which satisfy (7) and — in the sense of 
(3), (6) — generate stresses conforming to (5). 

With a view toward the particular thermoelastic problem described in the 
preceding section, it is expedient to introduce spherical coordinates (7, 6, ¢) by 
means of the mapping 


% =rsin6 cosé; x= 7 smi sind, 4 7c0s 0; 

(8) 

O'S p< co 0 0 So 0 ee aero 

Suppose the axis of the cone in question coincides with the x,-axis, let its 

apex be at the origin, and let 0) (0 < 0, <a/2) be its semi-opening angle. 
Further, it is helpful to introduce the auxiliary variables 


p=cos?, g=sind; p)=cos0,,° ¢5=sind, . (9) 


The region D is in this instance then given by 0<1r< oo, {p< p51 
(0 S py < 1), while the boundary B corresponds to p = y. 

By virtue of the rotationally symmetric character of the problem, the 
temperature field is independent of the longitude 4, so that the Laplace 
equation (1) here becomes 

VL, p) =O; (10) 
with 
0? Js 0) 1 0 ) 
pA re 
V Or? = yr Ov at re Op l(a ie mak: (is 


On the other hand, the boundary condition (4) to be met at present, appears as 


Lr, Po) =T (O<7r<a), Tir,p)=0 (a<r<oo). (12) 
In addition we require that 
T(r,p) +0 as roo (13) 
and that the temperature be bounded in D + B. 
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Two alternative forms of the solution to the foregoing heat-conduction 
problem were established in [1]. We cite here merely the integral version 


fp B~2Ey P_ypo,is(P) 
F(o, ~) = a p ap = =124%6 ae 
(2, 2) 2Vo J Sea Pa aaslPol dé (14) 
where 
0= 7, BlE,e) = cos (E loge), » (E, 9) = sin(E logo) , (15) 


and P_jj,;; stands for the Legendre function of the first kind of degree 
— 1/2 +76, which is real-valued for real arguments. 

In the case of axisymmetry about the x-axis, the general solution (6), (7) 
of the equilibrium equations (2) remains complete if g and yw; are restricted by 


P=" ~), Y=¥P), W=p=0, (16) 
whence (7) now reduce to 


Vip(r,p)==E tH 7 pray, p) =0, (17) 


in which the Laplace operator ’? is once more given by (11). Equations (6) 
when referred to spherical coordinates, pass over into 


*e 


(18) 
Ss eae, e 


HOR ame as 


The spherical displacement component w,, of course, vanishes identically. 
On substituting from (18) into the spherical counterpart of (3), one arrives 
at formulae relating the spherical components of stress to the stress functions 
yr, p) and yp(r, p). These formulae, after some simplifications based on (17), 
read: 


2(1+) po orp | 
aa eee Las LF or? 
ary ane 
(19) 
SU SB ae SCL Op 


Pap — — ae v2 Op? v2 Op y Or | 


eS aS ye 
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2.(Lsh+) Hit my PROS pean er 


ae te 1—y v2 Op y OY 
1 dy eee be is 
OG OZ 09 
be TSI E op Opor 
0? ow p op 
-pat +l 2v) 4+ 2(1 "= el 
The remaining shear stresses T,, and Ty 9 vanish. 
Condition (5) is met for the surface of the cone if 
Ta9(%, Po) =9, Trol% bo) =9 (0<41< 00). (20) 
To these boundary conditions we adjoin the regularity conditions 
Trrr Teer Tego Tra ~9 aS YOO, (21) 


as well as the requirement that the entire stress distribution be bounded in 

D+ B,i.e.for0 <7 < ow, f) SP <1. The problem under consideration thus | 
reduces to finding solutions y and w of (17) such that the stresses (19) obey — 
(20) and (21). 


3. Application of the Mellin Transform®) 


Let T(p, s), p(p, s), and #(p, s) be the Mellin transforms with respect to 7 
of T(r, p),r-* yr, p), andr y(r, p) , s being the transform parameter. Thus 


Tip.) = [Ti d) rar, | 
0. 8) = | ptr, B) ar, (22) 


Further, suppose 
¥ “ a 0 
a Os Peep ie tas i) — = Of) as ry-+oco (s, = OF Bi2ar 


°) See Dorrscu [11] for an exposition of the theory of the Mellin transform. 


ay a ee Rr 
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In addition, assume that the four functions appearing in (23) remain bounded 
as ry > 0. Then (22) and integration by parts imply the identities 


lo, 2) oo 


“yee . “-e 
as rs-2 dy = —(s — 2), / “or i dr = (s = 1) (s —2) P* 
0 0 

: . (24) 
a ae “ r Oy : 

hiraet 1dr=—(s—1)$, | —¥ wdr=s(s—1)¢ 


0 0 


provided s is confined to the strip 0 < Re(s) < s, of the complex s-plane. We 
now formally apply the Mellin transform to (17), after multiplying the second of 


(17) by r. Bearing in mind (11), (22), and (24), we obtain in this manner the 
Legendre equations 


(=p) FE] + (s=3)(s- 2) = EHS P 
isles] oo-s im 
py [O-*) B] + 6-2 6-No=0. 


Next let #,(p, s), %(p, s) be the Mellin transforms of r-! u,(7, £), 7-1 u(r, p), 
and denote by 1,,(p, s), Tga(p, s), etc. the transforms of the corresponding 
components of stress. Explicitly, 


(26) 
inelb, 8) =f trelt b) ar, ete 
0 
Equations (18), (19), because of (22), (24), (25), (26), now yield 
2pu, =—(s—2)p—Pist+2(1—2»)|p | 
} ; (27) 
P ay ay ~ 
2 ig =—9| ae +P ae — | (3 — 4) et 
4 = EE ee Fy (s—1) (8 -2)9 | 
ie thoes ae | 
@. | (28) 
x dp 
ton = — (6-2) +P ap 
Ae 
~(s—1) 6-1-2») ppt+-GB-2) AZ, 
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5, ee | 
= (12) GSP ee + (28) 

) ag 

p= 4 [6-1 2-0-2) 6-1 G+ 641-29) 0 FB]. | 

Finally, the boundary conditions (20) in the transform domain become 
Te9(Po 8) =O, T, a(Po 8) = 0. (29) 


The Mellin transform of the particular temperature field (14) is given by’) 


r Th je Poa 
Nb roe 


Taking account of (30) and bearing in mind that the thermal displacements and 
stresses sought must be regular along the axis of the cone, i. e. for # = 1, we 
take the solution of (25) in the form 


(30) 


A (1+) wa Ty a’ P,_(P) 
AOE s) = A(s) P,_a(P) (1 = y) Ss (2 oe 3) P._,(Po) ? (31). 


v(P, 8) = Bis) P._o(P) - 


The arbitrary functions A(s) and B(s) appearing in (31) are to be determined 
consistent with the transformed boundary conditions. To this end we substitute 
from (31) into the second and last of (28), and subsequently invoke (29). The 
solution of the resulting pair of simultaneous equations in A(s) and B(s) may be — 
simplified with the aid of the recursion relations$) 

g Py () Stag! lp P.(p) 2 P._4(0)] , j 
(32)*) 
(2s 53 1) p P.(p) = (s cig 1) P.i1(P) a s P._(P) : 
Setting 
(po, s) = 2 (1 — ») B§ — (s — 1) (s — 2) @, 
Abo, s) = —(8 ~ 1)? (5 — 2° hy [P.-a o)? | 
—(8— I) (5 ~ 2) 48 P,-alba) Pr alba) | 


a [2 (1 =) (8 =f) (s = 2) 0g) pL eee 
*) See Equation (10) of [1]. 
8) See, for example, [12], p. 161. 


®) Primes attached to Legendre functions are used henceforth to indicate differentiation with 
tespect to the argument shown. 


(33). 
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and making repeated use of (32 ), one may put the final formulae for ¢ Q, » into 


the form 
g—— tree Te wb Pe) 
l—» s(s—1yP s—1(Po) 
2(l+% v) wa Ty as 1+ qj (s — 2) © P,_3(P 
+ S (Po, Ss). [ Sel Fis is 
Tes {Po(s — 1) (s — 2) (s— 2») P. Ne 
+ 2 (1 — ») [1 +°q3 (s — 2)] P’_, (p,)} P. 
pa tiwahas—7 Ph) (a 
Bais s) (= >) Balbo) 
2 Po (s on 


+ Apa Pol — 2) Praalt) + BP alba. | 

A(po, 

Inserting (34) in (28), and once again using the recursion relations (32), we 
ultimately reach the following expressions for the transformed spherical stress 
components, in which 


E=2(1+»)p (35) 
is the modulus of elasticity. 
» _ EaTy,a’ : : 
tp = Sop ay Bop (2+ Bs PaO) — 2 PO) 
Po (s 


— Po) (b (s — 2) 2 bo (— 1) (8 — 2) Pals) 
— [2 (1 —») #8 — as (s — 1)] PL albo)} P.-al0) 

+ 9 {bo (s — 2) [2 ~ 5 (—1)] Pals) 

—2(1—» + a5 (s — 2)] Ph alpa)} PL_a(6))}. 


Be Ty a* 


B00 = ap a (Papy (a  - 2) - p(s 0 Bald) (36) 


SS eT Se Se fh 
: ; = 


$— ai(P 


— p (i — P*) a 

aay (6 (8 = 1) (6 — 2) {bo (6-1) (6-2) P,-al60 
+ [2 (1 —») 68-4 (6 — 1) (6 — 3)) Pl oP} Bo) 

— {ho (s — 1) (5 — 2) [2 (1 —») 9? — 9 5 — 2/9) Bab) 
+ [2 (1 =») 28-2 (1 =») (3-9) (6-1) 


~ 98 f* (s — 1) (= 2)] PLalho)} Pi_a(0))} 
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x Bede of : 1 2 2 1)2] P 2 pe 
TE36 ays (Po, s) \ P._1(Po) ([28 % (s ) ] ,-1(P) a % p 1(P)) 


u aa ((1 — 2) p(s — 2 (s — 2) [po (s — 2) P._o(b0) 
+ 3 Pl_o(b)] P.-a(b) — {Bo (s — 1) (s — 2) [2 (1 —) 
CT (s—3+2p) 


+ (1— 2») Bq? (s — 1) (s —2)] Pr alba)} Pi_a(2))h, 


| (36) 
tpg = at | - py (OO - Pal) + 3-29) Pal) 
Po (s — 1) P(e oe 
sae: ~A(Po 8) [(s — 2) {Bo (s — 1)? (s — 2) P._a(Bo) 
+ [2(1—») £5 + 98 (s — 1)] Ph_albo)} P,-2(0) 
— P{2 (L—») bo (5 — 2) P,-albo) + (2 (1 —») 
— gi (s 1) (6 ~2)] PL alba)} Pr_a(0))}. | 


The analogous formulae for %, and %, are omitted here since we shall not pursue 
further the explicit determination of the thermal displacements. 


In order to pass from the transformed stress functions and stresses to their 
antecedents in the physical domain, we appeal to the inversion theorem for the 


Mellin transform. By virtue of (22) and (26) we are thus led to 


¢+1400 
2 
vb) =a; | Obs) rds 
€+400 ag 
Y a 
WNP (Pp, s) r~* ds 
‘ ¢ +100 
ACE p) = Qni Z,,(p, s) rs ds , -etc. (38) 
The integrands in (37), (38) are fully accounted for through (34), . This 


completes the formal solution in complex integral form of the REESE 
problem at hand. The determination of the parameter c which enters through 


the limits of integration in (37), (38), will be discussed in the succeeding 
section. 


¢ 
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4. Real Integral Representation of the Thermal Stresses 


The solution just obtained is evidently not affected by the particular value 
of c as long as the position of the line of integration Re(s) = ¢ is varied within 
one and the same strip of regularity of the integrands in (37), (38). The selection 
of the appropriate strip of regularity, in turn, is dictated by conditions (21), 
along with the requirement that the stresses remain bounded as 7 + 0. 

If we take c = 3/2, the line of integration is given by 


Bey er 
Reet te (= bo < Fic eo) . (39) 


This choice of c is suggested by the observation that the Legendre function 
P._y, as well as the auxiliary functions m and A occurring in (34), (36), are all 
real-valued when s assumes the values (39), as may be inferred from the Mehler- 
Dirichlet integral?) 


_ 2 ; cos[(s + 1/2) 4] da 


r snes Te (Oren Searle (40) 
F = 


P.(p) 
Further, if Re(s) = 3/2, we note on the basis of the recursion formulae (32), 
that all Legendre functions entering (34), (36) may be expressed as linear 
combinations of the real-valued functions P,_, and P/_,, with coefficients which 
are rational functions of s. Thus 


(SS) 2) 
(¢= 2) P..sP)-= 6 (8 = 2). F.(P) +47 Pelf), (41) 
P!_,(p) = (s — 1) P,_2(6) + b Py_o(d) - | 


Also, the identities (24), underlying the developments in Section 3, are evi- 
dently valid on Re(s) = 3/2, provided s, in (23) obeys the inequality s, 2 3/2. 

We show next that the integrands in (37), (38), regarded as functions of s, 
are free from singularities if (39) holds, provided 


tl 
> 
a 
| 
= 
fy 
t 
= 
| 
& 
nw 
Ie 
7 
we 
= 
os 


1 
(“Coe aaa a (0<0 <3). -l<v<z. 


It is clear from the structure of (34), (36) that the singularities in question 
could arise only from the zeros of the functions 


A(s) = s (s — 1) w(bo, S) P,s(Po),  2(s) =s (Po, 8) A(Po, $) « (42) 


It is easy to see with the aid of the first of (33) that w(ppo, s) cannot vanish when 


sis restricted by (39), the same being true of the factors s and s — 1 in (42). 


10) See [13], p. 315. Recall that P,(1) = 1. 


ZAMP X1/31 
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The roots of P.(f,), considered as a function of s, are known to be real’). By 


(40), Pyo(Po) is positive for 0 S py < 1, so that P._,(po), and hence A(s), has 


no zeros on the line of integration (39). 

To confirm that Q(s) is also free from zeros we have yet to show that 
A(po, 8) does not vanish on Re(s) = 3/2. For this purpose we recall first certain 
properties of the associated Legendre functions of the first kind, which will also 
be needed at a later stage. For arbitrary complex degree s and positive 
integral order m, we have the relations”) 


ap™ ? 


Pri) = (-1)" Fete ea 


m (S1ye Risse mag 
I s (p) = - am im I (s — m + 1) 


Pop), 


[1+ O(9q7)] as 6-0, 


provided J* denotes the Gamma-function, and the recursion formula 
Prr**(p) +2 (m+ 1) © Bmp) + (s—m) (s+ m+ 1) Pm) =0- (45) 


Finally, the appropriate generalization of (40) reads 


—_— a) 
Poy tee ae gm cos[(s + 1/2) A] da 
Ve (p) Fi iz Tiajf2 = m) : (cosA in cos §)m+1/2 O<0<z, Re(m) << 1/2) . 
0 
(46) 


We now substitute for s from (39) into the second of (33), write A(fp, €) in 


P™(p) = (—1)" gm &lP) (43) 


(44) 


place of A[po, s(&)], and prove that A(fo, €) is negative. To this end we differen- i 
tiate A(po, €) with respect to 6, and find with the aid of (43), (45), and of the © 


Legendre eta that 


OA (Po, 
aha ! le oor 12 +i¢(Po) [CL ate aS Rice are ae atPe) 


: , ane)y 


+ 7) P* 1248(Po)] (<0) S2/2)% 
According to (46), 


P™ 2 +1e(Po) >0O (m= -—1, —3), (0 <0 5055 >| (48) — 


and, in view of the first of (44), these inequalities remain valid also for m = 1,3. 


11) See [14], p. 91. 
42) See ee p- 140 et seq. 
a) Here pl S124 and P? 


—1/2+i¢ are associated Legendre functions of order one and three, 
respectively. 


| 
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Consequently, (47) yields 


DA(bs, § | 
: a Pa (0< 6 <4), (49) 


while from (43) and the second of (44), 
1 2 
Apo, &) >— > (1+) (4+ 7)" <0 as 6,0. (50) 


It follows from (49), (50) that A(po, €) < 0 for 0 < 6, < 2/2 and — cw < & <oo. 
We are now in a position to conclude that the integrands in (37), (38) are 
indeed regular for c = 3/2, as was claimed earlier. Instead of pursuing the 
justification of this choice of ¢ further at present, we continue to proceed 
formally for the time being and subject the integrals (37), (38) to the change 
of the variable of integration (39). This results in a real integral representation 
of the solution, the validity of which will be examined subsequently. 

The highly cumbersome computation just alluded to may be outlined as 
follows. We first eliminate P._,, P,_3, P’_, from (34), (36) be means of (41), and 
refer (37), (38) to the new (real) sabsable of integration & defined by (39). The 
corresponding range of integration becomes — co < € < oo. We then decom- 
pose the transformed integrands into their real and imaginary parts, which are 
found to be odd and even functions of &, respectively!*). Consequently, the real 
parts do not contribute to the final integrals, which must of course be purely 
imaginary since y, y, and the physical components of stress, have to be real- 
valued. Moreover, since the imaginary parts are even in , the range of inte- 
gration may be reduced to the interval 0 < § < ow. 

With a view toward listing the final results in a reasonably compact form, 
it is expedient to introduce certain abridged notations. Let 


nO=F+5, 60 --Hee ree | (51) 
where o, B(&, 0), and y(é, ge) are given by (15). Writing 
P(E, b) = P_yeisselP) » (52) 
let 
XE pb p)= pee, YEP) =pepy Yl be) = YE bo bd) Se 


in which P’(é, p) stands for the derivative of P(é,) with respect to the 
_ argument #. Next, set 


x(&, 0; Po) a [Po - 2 % Y,(&, Po] B(é, 0) ia 2 Po eye; 0) 5) 
al, 0, Po) = [2 bo — 9% YolE, Po)] (E, @) — 2 Po BLE. @) - 


14) Note that P_4/2 , ;<(P)> and hence P“ 4/9 ,;e(P), is an even function of € since P_,_, = P,. 


(54) 
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Finally, with reference to the definitions (33) of «(fp, s) and A(fp, s), set 


wb) =o (bo s+i8), AE PI=A (do 5 +78), (55) 


and define new functions A,(&, po), 4o(§, po) through 


o(E, Po) A(E, Po) 


o(€, Po) » Adlé, Po) = ~The pF 


Ae Po = TRE bP 


(5 +i£)P W2+ik (bo) 


From (55), (56), with the aid of (33), the first of (41), as well as (51), (53), follow | 


o(E, po) = 2 (1 — ») £5 + 96 HE) 5 
Ax(€, Po) = @(&, bo) (£8 (6) — Po 96 YolE, Bo) + 96 YE(E, Pod] » 


(57) 


A2(&, bo) = (8, bo) {Po WP(E) + 95 (6) YolS, bo) + Pol4 (1 — ») 
+ 99 1(6)] Yo(E, bo)} - 


The desired real integral representation for the stress functions and stresses 
may now be put into the following form in which the auxiliary functions 
n, 6, X, Y, Yo,x,0,@,A,, and A, are understood to be referred to the corre- 
sponding arguments appearing in the defining equations (51), (53), (54), 
and (57). 


(0,0) _ aro 7 Be ee : 
Baty ar ane ewe sae [eX {1 — 248) 0 +453 
0 
2g g BY 2c y 
—=8¢8 is 
+ 9 (1 + 9 1) ¥}] - [262 X {bon —2 (1-9) ) a ¥) 
ee i ee = t 
Oe BY 


{Pon — Yo} + £22 Gonint+20—9)] 


—2(1—») (1+ @n) ie 


(9, P) a pained 
fo Pear oil? / 2 i a ac +2 * [be B ihe pe q ¥o}I} ag 
0 : 


ee: 


en i ara tae et mt will t s min 


PA tp oe 


= 
e9) 
On 
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a te fs 
TrrlQ,P) ¥ 2 4 2 = Dida) ay \ 
“EaTy v7) 7 o8!2 / \Ay |o{2 BPA Gil— 2 po). Y} 
0 


— {2 P32 X — g YY] — 2 [4B {hy g? VY — a8 1X Yo 

— (A= PV VY} + C{2popyrX— poy l(y+2(1+} Y 
+ Pm [95 (n + 2) + 2 (1 — ») 22] X Y, 

+2@[@n+1— ¥ ¥o}]\ dg, 

To9(0, P) i 4 


Bo T, 270°? / 
0 


— x {i bX + (8B) Y}]— 3 [2B Gan Y— ab 1 X Yo 


4: [a {5b X + [P§ — p+ 23 7] Y} 


+ 2 (L—») (Bj — p*) Y Yo} + C {pop 7? X + don (9? (y — 2) 
+ 2(1—y») pf) Y¥ — pm [go (n—1) + 2(1—») po) XY 


— [98 8 + 2(1 —») 23 p*)] ¥ Yo}]} ae, 


(59) 
| : ;. 1 2 2 2 
peeiee” * - ee / a: lot (gon + PO) X + [PN +1 
; 2: 
0 


— pj P*) Y} — x {b [96 4 + 5] X + 96 Y}] 
+ Po [2 B{(1 — 2») pog?n Y— (1-2) hd X Yo 


—2(1—» BY Y + E{1—2 pyopy?X 
+2697 (8+ (—» (1-—2@))VY+1—2rgGgpyuxX Y, 


= [an (1+ @) — 21 =») (QB gy + 1 = 2 95)] ¥ YO} aE, 


00, 2) Z a : 2(n+1)X—g@pY}—x 2 X — Br Y}) 
Pat, re 220" | (A, [a {b5 (4 + ) q P } {B35 0 


— 20 (2B {by 9? X — gy X Yo— b [gin +2 (1—»)) Y Yo} 
— E{pyn2?X —2(1—») pop n Y —nlg—2 (1 —») po] X VY 


— p (aq + 2(1—»)] ¥ Yo})} ae. 
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5. Convergence Considerations. Verification of Solution 


We now examine the validity of the formal solution in real integral form 
deduced in the preceding section. With a view toward establishing the conver- 
gence of the improper integrals in (58), (59), we consider first the asymptotic 
behavior of the corresponding integrands for large values of the variable of 
integration £. In this connection we recall the asymptotic expansions’) 


PO Veg reesay allt a)! (mae zl 
ager ay L(+ 2) 9+ (m+ 2) a] + Oe 


AS hse OO. “(eS Wisse ee Ole (60) 


NSAP, 6 eee i 
cee 


aS |S \->00 4, 1... teal, =m << ae 4S) ate 


From (60), (52), (43) follows 


P(g, p) = PED [14 58, + 0], - 


jzxéq | + BE 
ia ee 1) 3p ‘ 61) 
FE (é; p) Vee gi [1 8 Eg a3 Og *)| ( ) 
as E-oo («S@0a-¢€, e>0). 
Since, according to (43) and the second of (44), ' 
P&=1, Peéi=—3 (e+ 4), (62) 


(61), (53) yield 


X(E,P. bo) =//# expl-E (6 — 0)] 1 + OF), | 
X¢é, 1, Pp) = 2m GE exp (—é 6) [1 OG | 
VE,#, 6) =-|/ = expl-£0,-O1+0E), | (63) 
YEE 1, 6) = —|/7$" em exp(—£ 6) 11 + OF) | 
as &->0co Caves = Fei | 


4 
j 
45) See [12], pp. 147, 162, 47. : 
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In view of the continuity with respect to pat p=1 of X(E, p, H) and Y¢é, P, Po), 
the estimates (63) may be combined into 


ACE, P, bo) =O {E1? exp[—E (6) — 8)}}, 
YE. p, bo) = O {&5? exp[—E (6) — 6)]} | 


(64) 
as &->0o (0<950<%). | 
In addition we shall require the more refined estimate 
YolE, bo) = —F-[1- oP, + 01€)] i 
as &->0o (eS =5, e> 0), oie 


which is also derivable from (53) and (61). Finally, by (15), (51), (54), the first 
of (57), and (65), 


E(&, 0) = — sin (loge) + O(E) , 


rr) IS 


(5, Q, Po) = 2 [go cos(€ loge) — pp sin (€ loge)| € + O(1) , 
o(F, 0, bo) = —2 [bo cos (E loge) + gq sin (€ loge)] + O(E) , (66) 
OE, Po) = 9 & + O(1) 


as &— oo (0<e<x, e055, e> 0), 


while (65), (66), and the last two of (57) lead to 


A,(E, Ho) = 98 €* + O(8), As(&, Ho) = —% & + OE") 


FS (67) 
as — co (es <5, e> 0). 


The integrands in (58), (59) are continuous functions of € for 0 < § < co 


since A,(&, fy) and A,(&, £,) cannot vanish when 0 < 6) S 2/2, as is clear from 


(56) and the discussion following (42). Let £(€, @, , po) and gi, @, A, Po) denote 
any of the integrands in (58) and (59), respectively. The asyinptotic results 


(64) to (67), together with (15) and (51) now yield, as § > o, 


HE, 0, P, Po) = O {E-*? exp[—E (89 — 8)]} , 
g(é, 0, P. Po) = O {7 exp[—E (99 — 9)]} (68) 
(0<g<o, 0505%, 0<%=3). | 
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Therefore all of the improper integrals in question converge uniformly and 
absolutely for 0O<9< co, 0S95%—e(e> 09), O10 peste 2s 

To prove that the stresses (59) conform to the appropriate field equations, 
it suffices to show that @, y, and the spherical components of stress satisfy (17), 
(19). Substituting from (58), (59) into (17), (19), and formally interchanging the 
order of differentiation and integration, one is led to improper integrals the 
uniform and absolute convergence of which for 0 <@ <co,0 56 S6,—é- 
(c > 0), 0 < 0) <2/2 is assured by (64) to (67); this justifies the reversal of 
the two limit processes. In performing the required differentiations under the 
integral signs one needs to make use of the identities 


= Y, eo 2p Y+nX, (69) 
which are obtained from (53) with the aid of (52), (51), and of LEGENDRE’s 
differential equation. 

Both (17) and (19) involve the temperature field 7T(oe, #). The integral 
representation (14), while convenient for the numerical evaluation of the 
temperature distribution, is not suited for the purpose at hand. What is needed 
at present is the alternative representation 


lo-e) 


Te, ?)= zatve | Gb 0) X—oP YI AE, (70) 


which matches the structure of the integrals in (58), (59), and may be established 
as follows. It was shown in [1] that (14) is equivalent to the complex integral 


Berrie le | 
TON | Baie # (71) 


where c may be chosen anywhere in the interval 0 < c < 1, provided 0< 0)<a. 
When 0 < 9, <2/2, as has been assumed in the present investigation, one 
finds by precisely the same reasoning that the admissible range of c may be 
extended to0 <¢ S 3/2. 
If in (71) we choose c = 1/2, as specified in [1], there results (14); if, on the — 
other hand, we take c = 3/2, employ the first of (41), and make use of auxiliary | 
notation introduced previously in (53), (54), and (57), we reach (70). With the 
aid of (70), and by the procedure outlined above, the field equations (17), (19) . 
are eventually found to be satisfied identically throughout the interior of the — 
cone. The underlying computations are, of course, highly voluminous. It 
should be pointed out that the stresses (59) could have been deduced from the 
generating stress functions (58) directly by invoking (19) and without recourse 
to (36), (38). The fact that the physical components of stress (59) were actually 
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obtained through an independent inversion of their antecedents (36) in the 
transform domain provides an additional check on the correctness of the final 
results (59). 

In order to verify the boundary conditions (20) one has to go to the limit 
as Pp > py (9 > 4%) in t,4(0, f) and 7,,(0, p). Noting from (9), (53) that for 
P=bo ; = : 

eae! ee, Ct ae a (72) 


one finds after an elementary computation that the integrands in the second 
and the last of (59) are identically zero when p = 9, whence Teg and T, 
indeed vanish along the boundary 6 = 6). 

We consider next the regularity conditions (21). In this connection we 
observe on the basis of (15), (51), (53), (54), and (57) that each of the definite 
integrals entering (59) may be written as a sum of integrals of the type 

0-3? / F(é, p, Po) cos (& loge) dé , ee G(E, p, po) sin (E loge) d&, (73) 
0 0 
where F'(&, p, py) and G(€, p, po) are integrable functions of & for fp) << <1 
(0 S py) < 1). Hence, appealing to the Riemann-Lebesque theorem!*), we gather 
that all stresses tend to zero as 9 = r/a > ov. 

Finally, we turn to the requirement that the stress distribution be bounded 
“in the closed region 0X r<w,050< 6, (0 < 0 <2/2). In view of the 
convergence properties established already, we have to show first in this 
connection that all four stresses in (59) remain finite for 0 < @ < co, 6= 6). 
This condition is evidently met by Tg, T,g which were found to vanish for 
pf = po, consistent with the boundary conditions (20). Let £,r(, ©, Po) and 
&3(5, @, Po) denote the respective integrands of t,, and Tyg in (59), evaluated 
at p= p,. Then (72), together with (15), (51), and the asymptotic estimates 
(65) to (67), yield the conclusion 


i ] 
BerlE, 0. Po) = OE), Bg glEe, ba) =e) + OF) 


(74) 

as &->0o (0<o<o, 0<%&S5). 

Consequently, the integral representing t,, is uniformly and absolutely con- 

vergent for p=), 0<@< ow. Also, by DIRICHLET's test!”), the integral 

appropriate to T4(@, Po) converges uniformly for @ in any sub-interval of 

0 < @ < & which does not include the point g = 1. At g=1, Ty g(a, bo) has 
a finite jump-discontinuity which will be examined in the next section. 


16) See TircHMARSH [15], p. 403. 
17) See, for example, BRomwicu [16], p. 482. 
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To complete the validation of the solution under consideration we have yet 
to demonstrate that the components of stress remain bounded as 9 = v/a > 0. 
It is apparent from (59) that this is so provided all four integrals in (59) are 
O(o?’2) as 9 > 0. Unfortunately, an attempt to establish this order of magnitude 
meets with considerable analytical difficulties. The numerical results presented 
in what follows, however, indicate that the stresses are in fact bounded at the 
vertex of the cone. | 


6. Numerical Evaluation of Solution. Discussion 


The integrands in the solution (59) for the desired thermal stresses are seen 
to involve only elementary functions, with the exception of X(&, p, pp) and 
Y(E, p, fo), defined by (53). The transcendental functions X, Y, in turn, © 
contain the Legendre function of complex degree P(&, £) = P_1/2,;<(p) and its 
derivative with respect to p, P’(&, p). A computationally convenient integral 
representation for P(é,~) follows from the Mehler-Dirichlet integral (40) 
through integration by parts. Thus 


A 1 
Ul ool ee [cos @ (cosA — cos 6)1/2 
: 0 
75)18 
= (n + 2) (cosa — cos6)8®] cosh (EA) da (3) % 


— 
(O= 0) Ost agea, 
where 7 is given in (51). Similarly, taking m = 1, s = — 1/2 + 1 in (43) and 
using the first of (44), as well as (46), we reach 


ad 6 
apeag apical is J (cos — cos6)1 cosh (é 2) da 


cree (76) 
(O00 <0 Pools 
When ¢ = 0, Equations (75), (76) yield in particular 
P(0, p) = Pal) = = K(sin 4), 

\ 

(77)**) 
Pp 0, = [pe ee 2 . o = 0 . 0 | 
(0, 2) =P!) = —e [E[sin 7) — cost K(sing)], | 


18 A 

é ) oa of P(E, p) appear to be in preparation in the Soviet Union [17]. Observe that the 
integrand in (75), unlike the integrand in (40), has no singularity. Equation (75) may alternatively 
be deduced by setting m = — 2 in (45) and using (46), 


19) See also [12], p. 174 ; : 
reg [12], p - The second of (77) follows from the first by differentiation with 
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in which K and E are the complete elliptic integrals of the first and second 
kind, respectively. The integral representations (75), (76) are valid for 
0 < 6 <q; the corresponding values for @ = 0 (p = 1) appear in (62). Evi- 
dently, X(€, £, Ho) and Y(E, f, p,) are computable on the basis of (53), (62) 
(75), (76), and (77). hoe 
We recall from Section 5 that the integral representing T,,,(0, 9) is uniformly 
convergent in 0 < @ < oo, whence the boundary values of T,, are continuous 
in this range. In contrast, as mentioned in Section 5, Ts4(0, Po) has a jump- 
discontinuity at @ = 1, wehre the surface temperature T (0, Po) exhibits a dis- 
continuity of the same type. We now determine the surface discontinuity in 
ts 4- To this end we observe with the aid of (59), (51),. (54), (57), and (aa, after 
some manipulation, that 
a Ty (rvlO: Po) + ty (0, Poll = 


1 “A P wt V2 (78) 
ere k | /cer2arnes/ en | 
0 


0 


Bearing in mind (51), (15), and noting the identities 


r cos (c €) 
ee eee | 
0 
= (79)?°) 
r Esin(c € : 
| —“B aaa d& = * sgn(c) exp(—b)|c|) (b> 0). | 
0 
one finds that 
= liora< I 
1 @ 
ery Cab = h(1—o) =< (80). 
edie | Ditor oss 
where / is the Heaviside unit step-function. From (78), (80) 
resles Pa) _ (L+») bo [OLY ge tnlos Pd 
ie Te ae Leh aS ee 93/2 } A, ae Bok %) 


By (65), (66), (67), the integral in (81) converges uniformly for 0 < @ < oo 
and hence represents a function which depends continuously on @. Conse- 
quently, the discontinuity of t,4(@, Py) at @ = 1 is characterized by 


Tg (1+ » Po) — TE (1— , Jo) = EaTy, (82) 
as was to be anticipated in view of (12) and the plane-strain theory of steady- 


20) See [18], pp. 127, 128. 
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state thermoelasticity2). Equation (81), in which the discontinuous part of 
Ty 3(0, py) has been isolated explicitly, provides, of course, a more efficient . 
means for computing the surface values of t,, than does the corresponding 
discontinuous integral in (59). 

A procedure for the numerical evaluation of the thermal stresses may be 
outlined as follows. Let the opening angle 9) (0 < 0) <2/2) of the cone be 
given, along with the physical parameters E, y, and «. Once P(E, p), P(E; Do) 
and P’(é,p) have been tabulated on the basis of (75), (76), and (77), the 
values of X(é, p, fy), Y(E, P, fo), and Yo(&, Po) follow at once from (53) and 
(62). At this stage, the auxiliary functions A(€, @), y(§, @), n(§), C(E, @), #(E, @ Po), — 
o(E, 0, Po), (&, Po), Ar(§ bo), Aa(§, Ho) are directly computable from the ~ 
defining equations (15), (51), (54), (57), with the aid of the previously deter- 
mined values of Y). Thus the integrands in (59) and (81) may now be tabulated 
and one may proceed to the numerical evaluation of the corresponding im- 
proper integrals. 

In connection with the evaluation of the improper integrals just referred to, — 
it is helpful to have available asymptotic estimates of the truncation errors thus © 
encountered. Since the stresses on the axis and along the boundary of the cone 
are of particular interest, we shall confine our attention presently to the 
truncation errors arising when 6=0 (p= 1) and 0= 0, (f= py). Suppose — 
that t4(0, £; po, &o) and ti(0, P; fo, 9) denote the right-hand members in © 
(59) if the lower limit of integration (zero) is replaced with & (& > 0). Ap- — 
plying integration by parts to (59), after making use once again of (63), (66), 
(67), we ultimately reach the following estimates, which are valid as & > co — 
fore = 05. 0o:(s > 0): 


; Pa \fk : 6 o 
EQ, Ui bo Go) = Bo |/ 2 POEL HCa e+ Po 7lEo el exp (—g, 04) 


a OE explo byl 


. as q [9 B(So, @) — 1 0» 
Thal 1s bo bd) = — a a-Si Po More) — Voge Hn 8) exp (6, 6) 


+ 0 [85% exp(—E5 64), (83) 
i 
(0, Pai Po So) = oa {~g, [B Go 0) + 22 y (&,0)] 
+ oe loge Ci (| & logo |) 


+ [loge| [3 — $i ([éologe )]} + 015%) , 


1) See, for example, [19], p. 428, Equation (d). 
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where Ci and Si are the cosine and sine integrals defined by 
‘a : 


cifay=— [ SS aa, sifey =f 24 ar. (84)2) 


0 


In an analogous manner one obtains for the truncation error appropriate to 
the improper integral in (81) the asymptotic estimate 


wl ¥2 2 v(So @) | 
/ oot d= [ es loge Ci(|& loge |) + 0G, (85) 


WS “fg co “le So<oo, «> 0). 


For the purpose of illustrating the method of computation outlined in this 
section, the stresses on the axis 0 = 0 (f = 1) and at the surface of the cone 
6 = 0 (P = fo) were actually calculated for an opening angle of sixty degrees 


(0) = 2/6), Pp = V3/2, go = 1/2) and for vy = 1/4. To this end, Srpson’s rule 
was applied to the evaluation of P(&, 5) and P(g, #9) from (75), (76) for 
& = 0 (0-2) 20, with 50 subdivisions of the range of integration 0 <4 < 2/6. 
The computation of 1,,, Tg at p = 1 and of t,, at p = py was based on (59), 
while (81) was employed to calculate ts, at $= 9. In each instance the 
requisite improper integral was truncated at € = €&) = 20 and evaluated numeri- 
cally (once again by Srmpson’s rule) for g@ = 0-1 (0-1) 2:5, with 100 subdivisions 
of the interval of integration 0 < € < 20. The stress values thus obtained were 
corrected by means of the truncation-error estimates (83), (85). All of the 
numerical work described here was carried out on an I.B.M. 650 electronic , 
computer. The graphs displayed in Figures 1 to 4 are self-explanatory. 

In conclusion we merely mention that the solution to the thermoelastic 
problem of the cone established in this investigation, when 6) = 2/2, must be 


0-08 


10 net 


r/a —<— 


— 0-04 
Figure 1 
Variation of T,, along axis of cone. 0) = 7/6,” = 1/4. 


22) Numerical tables for Ci(z) and Si(z) appear in [20, 21]. 
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Figure 2 
Variation of tg9 along axis of cone. 0) = 2/6, vy = 1/4. 
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Figure 3 


Variation of t,, along surface of cone. 0 = 7/6, »y = 1/4. 
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Figure 4 
Variation of Td along surface of cone. O,) = m/6,v = 1/4. 
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reducible to the corresponding results for the half- -space reached in [9]. In 
particular it is easily seen from (59) that for 0) = 2/2 (Py) = 0) the normal stress 
T,,(Q, p) vanishes identically along the axis 0 = 0 (p = 1), in agreement with [9]. 
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Zusammenfassung 


Diese Arbeit befasst sich mit der Bestimmung der stationaren Warmespan- 
nungen, die in einem elastischen Kreiskegel durch eine unstetige Verteilung der 
Oberflachentemperatur hervorgerufen werden. Die Temperatur der Mantelflache 
wird oberhalb und unterhalb einer festen Entfernung von der Kegelspitze konstant 
angenommen. Das zugehérige Warmeleitungsproblem wurde in einer fritheren Ver- 
éffentlichung unabhangig gelést. In der vorliegenden Arbeit wird eine strenge 
Lésung fiir die gesuchten Warmespannungen mit Hilfe der Mellin-Transformation 
hergeleitet. Das hier benutzte Verfahren ist auch anwendbar auf das gewohnliche 
Randwertproblem fiir einen elastischen Kegel, der durch eine entsprechende un- 
stetige Oberflachenbelastung beansprucht wird. 


(Received: December 24, 1959.) 


Large Deflections of Circular Plates under Uniform 


and Concentrated Central Loads *) 


By Epwarp SAIBEL 2) and Irapj TapjBaKusH®), Troy, N.Y., U.S.A. 


1. Introduction 


In this paper we present a perturbation solution for the problem of large 
deflections of thin circular plates subjected to the combined action of uniform 
lateral pressure and a concentrated central load. The plate is assumed to be 
clamped at the boundary with no radial displacement, Figure 1. 

In recent years the problem of large deflections of circular plates based upon 
the von Karman [1]*) equations has received considerable attention. Various 
approximate solutions have been obtained for the case of uniform lateral loads 
by S. Way [2], CHIEN [3], BromBErc [4] and Retss and KELLER [5]. In some 
simpler instances the questions of convergence of the approximate solutions 
have been successfully treated [4,5]. 

For the present problem a simple solution beyond the linear classical 
theory seems to be lacking. Here we employ the well-known method of pertur- 
bation to obtain some nonlinear results inherent in the problem. 


1) This work was supported in part by the Office of Ordnance Research, U.S. Army. 
*) Dept. of Mechanics, Rensselaer Polytechnic Institute, Troy, N.Y. 


8) Formerly Rensselaer Polytechnic Institute, presently Inst. Math. Sciences, New York 
University. 


4) Numbers in brackets refer to references, page 502. 
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Nomenclature 


With the exception of symbols defined in the body of the paper the follow- 


ing nomenclature is used: 


r, 9 


w 


plane polar co-ordinates 
transverse deflection 
Airy stress function 
Younc’s modulus 

plate thickness 

plate radius 

Potsson’s ratio 

lateral pressure intensity 


a eS 


Figure 1 


Clamped circular plate subjected to uniform pressure and concentrated center force. , 


2. Formulation 


The von Kérmén equations for rotationally symmetric deformations of a 
solid circular plate assume the following form in polar co-ordinates. 


EW p 1 d (dg dw 
ym eat a7 Sir dy (4 ov) (1) 
E da /[dw\2 
Ab) + 55a (ar) = 9 
where J ans 4 
a *y de _ar> 


ZAMP XI/32 
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Each of equations (1) may be integrated once and the constants of integration 
evaluated by the condition of symmetry at the center, 7= 0. The resulting 
equations may then be simplified to yield. 


AV! =x | a6) dE +Vy, 
0 


y 2) 
32 x2 y"+V2=0, 
where we have introduced the dimensionless quantities 
eae Wie pay te 
a W(x) = 2 [12 (1 — »)] hn? 
3 (1 — d y dw 
v(x) eae is = if ie , V(x) = [12 (1 Fay y?) |e Th dy , (3) 


Ce hes at Ea iene ne 


and primes denote differentiation with respect to x. 

Now if the actual load on the plate consists of uniform load #, and a concen- 
trated central load P, we propose to express the pressure intensity function ~ 
appearing in (1) as 


2 
Pit) = Paar 20 (4) 
where 6(&) is the symbolic Dirac-delta function possessing the properties 


6(é)=0 for €40, 
A 
| 6G) d= 1 for A +0. 


=a 


It is thus seen that equation (4) provides a satisfactory definition of pressure 
intensity, since everywhere it reduces to that of a uniform load p, except au 


the center, where it remains undefined. Furthermore, it satisfies the load 
pressure relationship 


a 


T = total load = 2a r pry ar. 
0 


If one thinks of the Dirac-delta function as an even function of its argument 
then 


A 
[36 de=S, 
0 


| 
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and therefore, 
1 


1 
T= ma? | plr(x)) dv—aa*p,+2P | d(x) dx 
é 


0 
=2@p +P. 


Equation (4) with the aid of (3) can be written in the dimensionless form 


q=4, [1+ 2B d(x)}, (5) 
where 
 at[12 (1 — »)]32 ae 
u Pe) oe p,, B= tap, (6) 


Here #, the ratio of the two applied loads, may assume positive or negative 
values depending upon the relative direction of the application of the loads. 
Substituting (5) into (2) we arrive at the following system of equations 


V"=4q, (1+ px) +x*Vy, 
eae 
32 y" +2-2V2=0. | 


Boundary conditions at the center, x = 0, and edge of the plate, x = 1, must 
be specified to complete the formulation. From symmetry and regularity 
condition at the center we must have 


V(0) = y0) = 0. (8) 
If the plate is clamped at the edge with no radial displacement, we have 


V(1) = W(1) =0, | 


yl) =F" v0). | 


The radial and circumferential membrane stresses of” and o@' are given in 
terms of the new stress function p by 


ney os, 
“30 are ee 


(9) 


mn 
0, 


(10) 
Ene ote. 
Gn, as B a= y?) ae (2 y yx ) . 


The radial and circumferential bending stresses are given in terms of k by 


b WA, IB, Ae ee me V «x 
Hai pee’ ee (11) 
a eee eV Lay 4. 


a fiz@a— y?) ]3/2 a 
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3. Perturbation Solution 


The above nonlinear equations and the accompanying boundary conditions 
may be solved by a variety of approximate techniques. For the present pur- 
poses a perturbation process will be employed. We choose for the parameter 
of perturbation the maximum dimensionless center deflection 


W,, = W(0) , (12) 
and expand 
(GV, W) = 2) Ae Var Wo) Vn 
Wis Ln8 5.) ae (13) 
y= D/P Win 
n = 2,4,6,. 


where W, and V, are related by W,’ = «1 V, and A, are constants which define 
the load as an odd function of the center deflection. From (12) and (13) it is 


evident that W,0=1, W.0\=0,) 7=35 7 (14) 


The boundary conditions (8) and (9) carry over directly to the new functions 
V,, y, and W,. 

Substituting expansions (13) into equations (7) we arrive at a set of ordinary 
linear differential equations the first of which is 


Vi=A +p), 

Wi =e Ve | 

: (15) 
V,(0) = Vi(1) = W,) = 0; | 
W,(0) =1. 


Solution of this system leads to the following result 


il 
Vi= (528 +B xlogx—=x), 


» 
a7 
| 
| 
Jn 
| 
a 
——~ 
— 
a 


which is the linear classical solution of the problem (6). The succeeding V,,’s | 
and W,,’s satisfy linear equations of the form 


je at 
Was eer 


Vi=4, (1+ Bx) +26, 
| (17) 


V,(0) = V,(1) = W,(0) = W,(1) = 0 
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where 
G, = 50. ay as BOM on 


R=1, 3, 5 


A formal solution of this system, after some simplification, leads to 


“ 1 
Vn=4yVi— f fle, 81) Gyle) db + x f fC, & 1) GE) a8, 
0 0 
dy = dy [x og x G(x) dx, 
where 
f(x, €; w) = € 1 (u— x 4). 


Similarly the differential equations for p,’s are of the form 


Rh. ~y us = A 3 1 =- 
Py 39 a i =0 , y, (0) = y, (1) a ae y,(1) aa 0, 
where 
n—-1 
Hips SD MeV es ice 24,6, Se 


Be'4, 6, 


Solving system (21) we obtain 


ie iG E;1)H aes ft, e257) H,(6) a 


For ” = 2, we have H, = V?. Then from (22) 
RR a i. i 
oe ee hr [42(- oo x? log? x + 5 log ne a + ay 2) 


3] 2 log x = vse 
+B(-=* gets og x + 


CL: LAP ge Re: ) 
ee he ea a 
+( ea sag 


where 
t= Fy 16— 11» 5-30 
jae ae ee, ial oe ls 


Then from (18) we obtain for A, 


As 1 47 Jon | ee ~=) 
As = 64 |B (3% — a spray (22+ 4 “1 ~ 64 


2 ayy ‘) 
eta oe 


S01 


(18) 


(19) 


(22) 
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At this stage of the solution the load deflection relationship at the center 
of the plate assumes the form 


T=ak*p, (1+ 6)=P(1+ 8") 


8a Eh‘ (1+ B) W2). 
a® [12 (1 — »)]3 (Ay Win + As Ve) 


p-0 1 


45 

: 

i= r 

35 f 

30 ; 

" 

25 ' 
ee) 

bars i 

20 i 

5 

15 Bt 
\ 10 p-l 
5 oe 
ios B= ip 

0 


Figure 2 


Dimensionless uniform pressure g,. Dimensionless center deflection. 


As f > 0 we obtain the case of uniform load 


b= az — ae [Yn + Tyg (ty) ML, 


and as f + oo, a simple limiting process shows that 


Ges a “i ee [¥,, i = ( ; Oy al w3,| f 


Both of the above limiting cases have been previously obtained by variou 
investigators, e.g. reference (3) 
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The numerical work for explicit evaluation of the succeeding A’s becomes 
very laborious beyond this stage and can best be done by means of modern 
digital computers. The question of convergence of the series 


m 


Vu = LS We 


is a difficult and important one. It is evident that for convergence one must 
limit the range of W,, to 


W,> lim (= a)". 
n / 


In the engineering spirit and noting that the von Karman equations are 
valid for deflections of the order of the thickness of the plate, we have confined 
W,, to the range 0 < W,, <7. 

In Figure 2 is depicted the maximum dimensionless center deflection in 
terms of the dimensionless uniform loading for various values of the ratio of 
concentrated load to distributed load. 
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Zusammenfassung 


Die gewohnlichen von Karman-Foppl-Gleichungen fiir grosse Durchbiegungen 
einer kreisf6rmigen Platte unter kombinierter, gleichmassiger und konzentrierter 


_Belastung werden in Potenzreihen gelost. 


_ (Received: March 11, 1960.) 


— Se 


504 : ZAMP 


ee a 
Kurze Mitteilungen — Brief Reports — Communications bréves 


On the Calculation of the Minimum Wave plus Vortex Drag 
of Supersonic Wings’) 


By Raut R. HUNZIKER, Melbourne Beach, Florida, U.S. A.?) 


1. Introduction 


The calculation of the minimum wave plus vortex drag for a wing of given 
planform and carrying a given total lift at supersonic speeds poses the solution of 
the following mixed harmonic problem [1, 2]8) (see Figure 1) 


Foe aif 
ee ae at Ane 2 ; 
-1/ I; D i y u 
A 
B ¢ Ba Ee eae 
' 
! 
: 
Figure 1 ; 
Conformal Mapping of (D on upper plane 9. ; 
O*p 0*p 
yt ae aa 0, and p(y, 2) = y(—y, 2) over a domain D with | af 
symmetry axis z, | ; 


p=, over the curved piece J’, of the boundary joining the points A, D, C, (2) 
Op : 
yes 1, over the straight segment I’, (-1 < y <1, z = 0) joining A, B, C. 


Ayeoys p ; 3 

= The results presented in this paper were obtained in the course of research for CoNVAIR 
A peceion of General Dynamics Corporation, San Diego, Calif. 

‘) Department of Mathematics, Brevard Engineering College, Melbourne, Florida 

) Numbers in brackets refer to References, page 507. 
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The boundary ’= I’, Uv I’, of D corresponds to the orthogonal projection for 
2 = Ovof the forward and rear Mach envelopes intersection on the plane y, z trans- 
versal to the direction of motion. The segment J’, corresponds to the projection of 
the wing planform of span 2b. 

The mixed harmonic problem defined by Equations (1), (2), (3) can be solved 
by determining a complex function 


F(s) = p(y, 2) + 7 ply, 2) , (4) 


which is analytic regular over D + J’ with a vanishing real part over J’, and with 
an Imaginary part equal to 


vy) =—(a+y) over F,. (5) 


Condition (5) follows from a Cauchy-Riemann equation. 

This holomorphic problem is a particular case of a more general one, related to 
a contour consisting of m pairs of connected pieces, which was treated by S1GNo- 
RINI [3] who showed that certain m — 1 definite relations are necessary and sufficient 
for the existence of the solution. The same problem was treated later by KELDYSH 
and SEDov [4] and more recently by Komatu [5]. 

Komatu [5] has given an integral representation of the solution for the case of 
a single pair of arcs of a circumference bearing respectively the prescribed boundary 
values of the real and imaginary parts of an analytic function to be found. Although 
the kernels contained in this representation are expressed by elementary functions, 
an integral formula which is more tractable for the effective calculation of the 
boundary values of the solution gm over J, can be derived. The continuity and 
boundedness of the boundary values (2), (5) assures the existence of the solution 

_F(s) without any restriction as in the case m > 1[5]. The uniqueness of F(s) 

follows from the assertion that if a function is analytic and bounded in the domain 
@ and has real and imaginary parts which vanish respectively over the pieces I’ 
and I’, of the boundary contour I’ it must vanish identically. The analyticity of 
F(s) in the symmetric domain @ and the symmetry condition of Equations (1) 
imply, by ScHwarz’s reflexion lemma [6] the antisymmetry of p(y, 2) with respect 
to y = 0 and the identical vanishing of y over B D. Therefore a = 0 is taken in 
Equation (5). ‘ 


2. Conformal Mapping of 


Let 
w(s) = u(y, 2) + 7 0(y, 2) (6) 


be the complex analytic function which maps conformally the domain © on the 
upper half plane @ of the complex variable w. a ad 
Point A, which is not shown in Figure 1 is chosen as the point of infinity of 
the w-plane. Since the two plane domains (D and @ can be mapped on a circle 
‘|p| <1 they can be mapped on one another. Since the domain 9 is mapped by 
the inverse analytic function s = s(w) on the domain @® of the s-plane, the solution 
EF (s) is translated into a function F(s(w)) = p(w) +1 p(w) analytic and bounded in 9. 
‘As is known several methods for the effective construction of the mapping 
function of D on the unit circle have been proposed [7] of which THEODORSEN’S [8] 
has been the most generally used. A comparative review of the various proposed 
methods using high speed computing machines is given in Reference 8[]. 
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The mapping function of the circle |p| <1 on the upper half plane Pv = 0) 
is a linear function containing three real parameters which can be determined so 


as to have I’, mapped on v= 0, —- co <u S0. 


3. The Singular Integral Equation and the Calculation of p(s.) over I, 


The function 
y i , 
K(w) = yo. a= KR kG ”) 
w 


is analytic in the upper half plane 9 and 


Ia O) - kK; = —_, v=0, u< 0; , 
V\~| (8) 


Kp = — Go Oy Oy 0 


Let us introduce now the function 


2(w) = F(w) K(w) = G(w) + 1 H(w) , 
where 
Gw)=pKe—ypK,, Hw) =~ K+ yp Ke. (9)% 


By Equations (2), (5) and (7) to (9), it follows that the values of © on the real 
axis v = 0 are given by : 
G(u) = vis u<0O 


Viel’ 


and because mg = 0 on J, which is mapped on u > 0, 


Guy 0 n= 0: 
similarly it follows 


fe : 
Hw) = ay Oe ee . 
Ve 
Let 2 be a closed rectangular contour consisting of a part Q’ and two segments 
of Ti), Fay s uSe, e <u <y and a semi-circumference x of radius ¢ 
with center at the singular point w = 0 of K(w). Since F(s) is by hypothesis analytic — 


regular over D+ J, the unknown function F(w) will be also analytic regular 
within and over the contour 2 and the singular integral equation 


, 


1 f 2(w) 


Tt W— Wo 
Q 


| 
L(w) = dw , 10) 


where w and W» are on 2 will relate the boundary values of the function Q(w) 
which is analytic regular when w is in the domain of contour Q [6]. This singular 
integral equation is the necessary and sufficient condition for (w) to express the 
boundary values ©(w») through which it is defined as analytic regular in the 
domain interior to Q. 

By the maximum-modulus theorem and since F (w) is analytic regular within 
and over 2 it follows | Q(w) | < M | w |” which shows that Q(w) is of the order 
of w” as w tends to infinity. Then the integral over Q’ vanishes for y > 7 

: 
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Since the integral over x vanishes with ¢ and by taking 


arn : the imaginary part of the 
limit of the integral for y + oo it follows i i 


where 
o(4) =x 4Bw(x), x=-—u, O0S%<o0. 


The function (— +9) is expressed as the principal value of an improper integral 
containing a real parameter x, and it can be shown to be uniformly convergent and 
therefore continuous in the interval 0 < #% Ss 00 [9, 10, 11). 

Since the function (uw) vanishes at the points C, and A, the function g(s) over 
is continuous and vanishes at C = (1, 0) and 4 = (— 1, 0) as required by the boun- 
dary condition (2). 

Since the Cauchy principal value exists, the numerical computation of g(— ¥5) 
at a number of points +) which are transformed onto the axis y allows the numerical 
calculation of 


1 
| oo, y) dy. (11) 
0 


4. Conclusions 


The minimum wave plus vortex drag coefficient C,,, of a wing of fixed 
planform and lift coefficient C; can then be calculated by 


: ce, 


Dop. — x 
16 b? | @(0, y) dy 
j 


_ where S,, is the total area of the wing of span 2b [1]. 

This procedure can be coded for a computer and applied to wings of general: 
planform for which a symmetric domain © is defined. 

As compared to another procedure that has been suggested [2] the functional 
theoretical basis of the present method assures that the numerical values of (0, v) 
and its integral (11) correspond to the unique solution of the mixed harmonic 
problem (1) (2) (3). 
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Résumé 


Les valeurs de contour de la solution d’un probleme mixte harmonique sur l’are 
portant des données de NEUMANN sont obtenus par une représentation conforme 
et par évaluation d’une intégrale généralisée. Cette méthode permet de calculer le — 
minimum du coéfficient de trainée d’une aile supersonique. 
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Bemerkunégen zur glatten Interpolation . 
' 


Von HeEInz RUTISHAUSER, ETH, Ziirich. 


Die Interpolation in einer durch m Funktionswerte y, fiir gleichabstandige — 
Argumentwerte 4, = %) + k h gegebenen Funktionstabelle wird gewohnlich folgen- 
dermassen ausgefiihrt: Fiir gegebenes x wird zundchst festgestellt, in welchem der : 
Intervalle I, = (+, ¥,4,) das Argument liegt. Alsdann wird die Funktion # durch — 
ein Polygon P, vom Grade 27 + 1 ersetzt, welches mit f an den 27+ 2 Stellen 
ey Ne-r+1> +++¥pyr41 Ubereinstimmt, und dann der Wert von P,(x) anstelle 
von f(v) genommen (siehe Figur 1 fiir y = 2), Die numerische Berechnung ol 
P,,.(x) erfolgt zum Beispiel mit der Besselschen Interpolationsformel. ' 

Es ist nun schwer einzusehen, warum man von P,(¥) Ubereinstimmung mit f(z) 
an den Stellen ¥,_,, ... %, 4,41 fordert, wahrenddem doch die Werte von P,,(¥) nur 
im Intervall I, gebraucht werden. In der Tat wird dafiir ein hoher Preis bezahlt, 
indem die durch dieses Vorgehen definierte Approximationsfunktion fiir f(x), die 
aus Polynomstiicken (P, in I,,) besteht, bereits eine unstetige erste Ableitung 
besitzt. 

Es scheint daher zweckmassiger, zu fordern, dass die Approximationsfunktion 
als Ganzes s-mal stetig differenzierbar sei, dafiir aber zu gestatten, dass die einzelnen 
Interpolationspolynome P,(x) (vom Grade 27+ 1) nur in 27 — 25+ 2 Punkten 
Vp_pts +++ pip gi, Mit f(%) iibereinstimmen. 
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509 
Wir betrachten die GauSsche Interpolationsformel der Ordnung 27 + 1 
t) 2 a a A Pan es IY x 
fess (1) 
eas 2r+1 
* a a :) On 1/2 


wobei ¢ = (~ — x,)/h. Die entsprechende Formel fiir das vorangehende Intervall 
T,,_, erhalt man, indem man & durch k — 1 und ¢ durch t+ 1 ersetzt. Sie kann 


qi 
re F 


£ 


Xq-2 Xa Xx Xa X42 X03 


Figur 1 


Lage der Stiitzstellen relativ zum Interpolationsintervall I; . 


bekanntlich wie folgt umgeformt werden (mit der gleichen Bedeutung fiir x wie 
oben): 


t eae ee 
ey 24) — f+ (;) 0, —aja + ( 2 ) On + 


i+ 1 3 sane 27 
( 3 ) tae os tf Or 


(Bee 2741 
* é Y + Mt On 1/2 " 


Damit lautet die Differenz nach einigen Umformungen: 


Axle) = Pla) — Pa) = [¢+ (3)- ("5 ')] 4 


1 cok fe Ov 4 
etsy. )-C]# @ 
ttr—1 ttyry—1 - gaol or feet) aa 
+t [oy + | 2r age} oh ® Wa paca) ee 
Da dabei alle eckigen Klammern verschwinden, erhalten wir 


Orie, 27 +2 3 
ce Fgh Oa (3) 
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Wenn 62’*? + 0, hat daher P,(¥) — P,_1(#) eine nichtverschwindende erste Ablei- 
tung an der Ubergangsstelle + = x,. Dies mag unbedeutend erscheinen, wenn Cree 
klein ist; jedoch hat der Verfasser diese Unstetigkeit der ersten Ableitung der Ap- 
proximationsfunktion an den Ubergangsstellen schon oft als sehr st6rend empfunden. 

Zur Verbesserung dieser Situation addieren wir eine Korrektur zu dem Inter- 


polationspolynom fiir das Intervall [;, 


Q,(#) = P,(x) + elt) 02733, (4) 


(mit der gleichen Bedeutung von ¢ wie oben) und wahlen das Polynom g(t) vom 
Grade 2 + 10, dass Q,(*) an das entsprechend gebaute Polynom fiir das Intervall 
I;,_1 mit s stetigen Ableitungen anschliesst; das heisst, dass die Differenz 


By beet Qr+1 euat 2741 
4,@) = [00 + (2,°,1)| Hab [PC+D4+(, 4 a) aGk ©) 
an der Stelle x,, mitsamt ihren s ersten Ableitungen verschwindet. Man erhalt nach 
(3), (4): 

re bay 2741 a Qr4d 
4,@) = [90 + (2,°,",)| ak [eG +9460, 4 aan © 


Um die gestellte Forderung allgemein zu erfiillen, miissen wir offenbar verlan- 
gen, dass gleichzeitig : 


[eo + (04 G)] ane [e + te peal 


mitsamt den s ersten Ableitungen an der Stelle ¢= 0 verschwinden. Das erfordert, 


wenn wir 
ON at : 
(5) ke ot 5 mite oe, (6) 


gp (0) = (1) = —Cry (uw = 0,1... 5). (7) 


Da @ ein Polynom vom Grade 27 + 1 sein soll, ist g durch (7) noch nicht festgelegt, 
wenn s < y. Wir kénnen daher weitere Forderungen stellen, zum Beispiel 


bezeichnen: 


Om) =0 fir n=kh-r+s, R—r4is #1, ck by Sse (8) 


Dies sind 27 — 2s zusatzliche Nullstellen (0 und 1 sind unter (7) bereits beriick- 
sichtigt). 
Da nun 


t+y te) eee 
2r+1 (27 + 1)! 
eme ungerade Funktion ist, werden die c,,, fiir gerades yu verschwinden und auch die 


Nullstellen (8) zu ¢ = 1/2 symmetrisch legen. Dies bedeutet, dass das durch (7) 


und (8) eindeutig bestimmt ; : bea . 
duc, Fetches 8 mte Polynom (¢) antisymmetrisch beziiglich ¢ = 1/2 ist, 


p(t) + o(1—t) =0. (9) 
Wir machen nun den Ansatz 


ht | 
fled PAR kT, (10) 


vc OC ee : 
i Ce ee eee ee ee ee 


we eS 
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wobei wir vermerken, dass infolge der Symmetrie von ( a of 
2¥—2s+2 
mmetrisch sein muss. Da 


) beziiglich 
t= 1/2 auch das Polynom y(t) vom Grade 25 — 1 antisy 
anderseits 


et? pai sited Ce er elt Soh te Bhs heh ee stad) t+r—=s 
2r+4 (2¥+ 1) (2%) (2r— 1) (27— 2... (@r—25 4 3) set re | 


und at) 


SSS COR Bi SS diem Uae9) CH) Haff ee 5-1) 
v(t) Blea mts fxs. s) py) + (27 +1) ChE CEE a 


fiir ¢ = 0 mitsamt den s ersten Ableitungen verschwinden sollen und ( y + i ) 


: 2Zy¥—2s+4+2 
bei ¢ = 0 eine einfache Nullstelle hat, folgt: 

Der Inhalt der eckigen Klammer bei (12) muss fiir t = 0 mitsamt seinen s — 1 ersten 
Ableitungen verschwinden. Zusammen mit y(t) + y(1 — t) = 0 ist das Polynom w(t) 
dadurch eindeutig bestimmt (bei festem 1, s). 


Damit liefert dann 


y— Ss 
Onl) = Bia) + vO (, a ome) ein (13) 
eine Interpolationsformel fiir das Intervall J, die mit dem entsprechend gebildeten 
Polynom fiir die Intervalle J,_, und J,,,, mit s stetigen Ableitungen anschliesst. 
(Dabei bedeutet P,, das exakte Interpolationspolynom vom Grad 27 -+ 1, welches 
mit f,(¥) an den 2r-+ 2 Stellen ¥,_,,...% 4,41 tibereinstimmt; es kann nach 
irgendeiner der Formeln von Gauss, BESSEL oder NEwTon berechnet werden.) 

Dieser Erfolg konnte offenbar durch Abanderung nur des letzten (CSGAN -Gliedes 
der Besselschen Interpolationsformel erzielt werden, und darum sinkt auch der 
Grad der Approximation nur um 1. 

Die bisherigen Betrachtungen zeigen, dass man solche Korrekturterme fiir be- 
liebiges y (wobei 2y + 2 die Anzahl der Stiitzstellen und 27+ 1 den Grad des 
Interpolationspolynoms bedeuten) und fiir jedes s Sy (s = Anzahl der stetigen 
Ableitungen beim Ubergang von einem Intervall zum nachsten) bilden kann, doch 
werden die Formeln mit steigendem s rasch komplizierter. Wir behandeln im fol- 
genden die Falles = 1 und s = 2. 


Konkrete Beispiele 


s = 1 (stetiger Anschluss der ersten Ableitung): 


Hier muss offenbar 


(+7) 
Ad 527 
fiir t = 0 verschwinden; also ist 
v y 
p(0) = oe und y(l)=+ ee: 


oder 
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Das Korrekturglied wird also hier: 
‘ Pima feed see 
(=) Fo ae 
Nun kann man dieses aber mit dem letzten Glied der Besselschen Formel 
me. i fat pert 
(f= 24 i 2r k+1/2 


zusammenfassen und erhalt: 


Satz: Durch Weglassen des Nenners 27 + 1 beim letzten Ghied 


ae 1 ae BEGET AF 
( Bi BGs 2 k+1/2 


ZAMP 


der Besselschen Interpolationsformel vom Grade 2 7 + 1 erhalten wir einen stetigen 


Anschluss der ersten Ableitung in diesem Intervall zum nachsten. 
Beispielsweise erhalt man fiir y = 2 die folgende 6-Punkte-Formel: 


t 2 t— 1/2 /t\ «3 
O,(*) = hr + #6, . 4/0 == (;) Oe a -7 a3 (,) On 1/2 
bor dl 1\ (¢+ 1) 35 
ot: ( 4 ) # Sapo it (:— Z| ( a ) Braye 


wobei Oe = 1/2 (6;” ao oa) bedeutet. 


Auf die Funktion e* mit den Stitzstellen —3, —2, —1, 0, 1, 2, 3, ... 


wendet, erhalt man folgende interpolierte Werte (und ihre Differenzen) : 


a) mit normaler Besselscher Formel b) mit Formel (16) 

vom Grade 5 
0,8 2,23148 0,8 2,23849 
D.BBirayegy34432 010) let b68 0,85 2,35074 a oat 
0,9 PAGIBS ir cP ie AGEN! 0,9 2.46793 (eerie 
GOS Wa stIso 050) cea 0,95. 2,50035 5 alee 
1,0 Z1T865 So tea 1,0 2,71828 ais 
1,05  2,86168 57-18 698 1,05. 6 2.853640 oe eee 
1,1 301206 20 | aes is 2.00824 | = ee 
Lee AS16972 57 761 1,15) 3452208 248 
1,2 3, 33499 527 1,2 3,31594 pO, 


Obwohl die Werte unter b) keineswegs genauer sind als die unter a), kann der 
glattere Verlauf der interpolierenden Funktion bedeutende Vorteile bieten. 


s = 2. Hier ist offenbar zu fordern, dass 


= BE) ta) ae 
” (27+ 1) 27 (ary 


an der Stelle ¢ = 0 mitsamt der ersten 


4904 
523 
5a 
743 
924 
945 
960 


Ableitung verschwindet, was zusammen mit 


ee ee 
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der Antisymmetrie das Polynom y(t) vom Grade 3 eindeutig bestimmt. Man findet: 


¥ 1 et (li 
= _ = Dace 
AO garg 2) (! ; pee 


: y 1244/2 fist y — 2 
i= 7 |e? —. 7 — wes d 
v(t) r( Forth aes (oye 


(17) 


Auch hier kann man die Korrektur wieder mit dem letzten Besselschen Term 
zusammenfiigen. Dieser ist 


bo Pb U2 ara 
( 2r eae R+U2 


und also der Koeffizient des 0,” 4 y-Gliedes in der modifizierten Besselschen Inter- 
polationsformel 


(: = 3) (‘ +r — 4 [(272—2r4 1) (¢2 =f v) pert (18) 


ee + 1/29 


wahrenddem die iibrigen Glieder unverdndert bleiben. Damit erhalt man zum 
Beispiel mit y = 3 eine 8-Punkteformel (mit der iiblichen Bezeichnung von #) 


; (t 2 t+ 2 a) 
On(*) = fy + tO, 442 + (5) # Op aja bo + ( 6 jn Of 41/2 
i: (19) 
1 t+ 1 
+ (#= 3) ( i )u3se- ig thee ee 


die eine durchgehend zweimal stetig differenzierbare Approximationsfunktion 
definiert. 


Summary 


The usual interpolation formulae for equidistant abscissae (NEWTON—GREGORY, 
BEssEL) define a piecewise polynomial approximation function whose first deriva- 
tive is generally discontinuous at every meshpoint. It is shown how these formulae 
can be modified (without using higher derivatives of the given function) such that 
the piecewise polynomial approximation function has s continuous derivatives 
(where s is a given integer). 


(Eingegangen: 2. Juli 1960.) 
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Eine Methode zur numerischen Losung einfacher nichtlinearer Rand- 
wertaufgaben mit Anwendung auf ein Randwertproblem 
der Himmelsmechanik 


Von F. SCHLAEPFER, Ziirich?) 


Es sei eine Funktion x(#) gesucht, welche die nichtlineare Ditfferentialgleichung 
a*x 
dt? 


befriedigt, wobei /(x, #) eine vorgegebene Funktion ist, und die inhomogenen Rand- 
bedingungen 


= — f(x, 2) (1) 


e(—1) =a, waA#1)=6 (2) 
sind’). 
Mit der Einfiihrung der zur homogenen Randwertaufgabe gehérigen Greenschen 
Funktion K(é, t) wird das Problem 4quivalent mit der Integralgleichung 


x(t) = x(t) + [Kt ©) feces) de (3) 


wie man durch auf der Hand liegende Erweiterungen der Ausfithrungen in [1]?), 
Seite 306ff. erkennt. 
Dabei ist 
a+b =) 


BGA) = 2 ee i (4) 


also die Gleichung der Geraden durch die gegebenen Randpunkte und 


| (1 = 7) {tf —) 


. nig ees 
2 fiir 1 Se 


| La 


“IA 


oS: | 


K(t, t) = (5) 


2} 


Se ee el. 


2 Kennt man eine Naherungslosung **(¢) fiir x(t), so kann man diese auf der rechten 
patie von (3) einsetzen und bekommt dann auf der linken Seite eine bessere Naherung 
«**(t), die man wiederum als neues *(¢) auf die rechte Seite bringt. Auf diese Weise 


1) Institut fiir angewandte Mathematik, ETH. 
5 : 2 a ; 
) Das hier beschriebene Verfahren lasst sich ohne Schwierigkeiten auf Gleichungssysteme 


ax, 
dz a —hi(%, Na, +26 5%, t) 


und Randbedingungen 
4i{—1) = a, 241) = b,5 f= eee 


fe eSecinyplen Ebenso lassen sich, sofern bereits brauchbare Naherungslésungen zur Verfiigung 
stehen, auch Faille, wo f noch von der ersten Ableitung abhangt, behandeln 


RB oo ; ie , 
) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 517. 


oe 


i] 
f 
5 
5 


. 
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verfahrt man so lange, bis sich zwei auteinanderfolgende Naherungen fiir alle ¢ im 
Intervall nur noch um einen kleinen vorgegebenen Betrag unterscheiden und erhalt 
so die gesuchte Lésung. Verfiigt man iiber keine spezielle erste Naherung, so kann 
man eine solche immer bekommen, indem man das kubische Polynom wahlt, das bei 
¢—= — lund? = + 1 mit derexakten Lésung in Funktionswert und zweiter Ableitung, 
die sich durch (1), (2) berechnen lasst, iibereinstimmt. 

Die Hauptschwierigkeit besteht in der numerischen Berechnung des Integrals 
in (3). Da K(é, t) bei ¢= rt eine unstetige erste Ableitung besitzt, lassen sich die 
gewohnlichen Newton-Cotes-Formeln (vgl z. B. [2], Seite 139ff.) nicht anwenden. 
Wir miissen deshalb spezielle Formeln herleiten, indem wir nicht wie normalerweise 
den ganzen Integranden, sondern nur f(x, ¢) allein durch Polynome n-ten Grades 
in ¢ approximieren. 

Eine Naherungslésung x*(¢) sei an m+ 1 verschiedenen Stiitzstellen Let = 0, 
1,...,, deren Lage im Intervall (—-1s#,<5+1) prinzipiell beliebig ist, bekannt, 
und es sei 

f= fe"), 4), t=0,1,...,0. (6) 


Bezeichnen wir mit /, ,.(¢) das Polynom n-ten Grades, fiir das 
Lx) = 6,, firialle,.¢= 0,1, ... , 2 


gilt (Lagrange-Polynom), so ist 


nn 


D te Inelt) = Pal) (7) 


k=0 


das Polynom vom Héchstgrad », das an den Stellen ¢; mit /(+*, ¢) iibereinstimmt. Die 
verbesserten Werte «**(t,) berechnen wir mit Hilfe von Formel (3), aus der 


+1 
4**(E.) = ¥o(t;) + [Kt Teale ate 0 sacl es oer, (3a) 
ae 
folgt. Dabei approximieren wir das Integral auf der rechten Seite durch [vgl. (7) | 


+1 +1 


| Blt) Pale) de = She | Kis 2) ban) = Se Ae (3) 
a k=0 oor k=0 


Die Koeffizienten 
i= 


Cheats att 
; (9) 
k=0,1,...,n | 


+1 
Aix =] K(é;, T) La r(T) dt , 
-1 


sind von der Funktion f(x, ¢) vollkommen unabhangig; sie sind einzig gegeben durch 


die gewiinschte Ordnung » des Rechenverfahrens und durch die Lage der Stiitz- 


tellen ?;. at 
& Nach Formel (5) ist der Kern K(¢,, t) in den beiden Teilintervallen (—1 St S7;) 


und (¢; <t <1) je eine lineare Funktion in t, der Integrand in (9) also je ein 
eo 


o 
2 


Polynom (n + 1)-ten Grades. Aus diesem Grunde kénnen die A, analytisch be- 
rechnet werden. 


nenas-ants 
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Es seien noch einige Eigenschaften der A;, erwahnt. Setzt man in (8) P,(t) =1 
ein, so erhalt man unter Benttzung von (5) 


pe eo 


n 1—#? 

> Ain = i K(t;, t) dt arr t= Opie ee (10) 

k=l ite 
Sind die ¢; symmetrisch zum Ursprung verteilt, das heisst ¢,_, = —t, fiir alle 1, 
so gilt 

Ain = An_i, n-k- (11) r 

Weitere Vereinfachungen ergeben sich, wenn Stiitzstellen mit Intervallsgrenzen — 
zusammenfallen. 

Fir 4.=—1  folgt. A,6 = 05 t= Ql eee (12) 

Fir #,= 1) folgt A;, = 072 te= O01) ee 


Aus (3a) und (8) bekommt man schliesslich die endgiiltige Iterationsformel 
n 
x**(t) = x(t) +) Ante» += 0.4, 2pm: (13) 
k=0 


Der erste Teil rechts in (13) ist nach (4) die lineare Interpolation zwischen den 
gegebenen Randpunkten, und der zweite Teil lasst sich als Produkt einer Matrix 
(A;,) und dem Vektor der Funktionswerte /, auffassen. 
Beispiel einer Matrix (A ;,): 
Fiir » = 4 und die 5 aquidistanten Stiitzstellen 
1 i 


i cs rae A _=0, =>, 


=, fy = 1 


liefert Formel (9) die folgende Matrix (A,,), welche alle Eigenschaften (10) bis (12) 
besitzt 
0 27 16 7 | 
0 332 256 U32Z 
| 0 222 416 22, 
| 0 132 256 SW 
0 7 16 25 


SVS) SS 


Nach Beendigung der Iteration kann man auf einfache Weise die in vielen Fallen 
wichtige erste Ableitung von x(t) an der Stelle ¢ = —1 berechnen. Aus (1) und (2) 
folgt namlich 
+1 
dx 1 
F (- =F fa Ce a 


=1 


b—a4 
y > 


(14) 


wie man durch partielle Integration verifiziert. Das Integral in (14) kann man mit 
gewohnlicher Quadratur berechnen. 
Anwendun§gsbeispiel 


; Sind zu zwei gegebenen Zeiten die heliozentrischen Standorte eines Himmels- 
kérpers bekannt, so lauft die Bahnbestimmung auf ein Randwertproblem vom 


a ee 
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Typus (1), (2) heraus. Bucerrus (vgl. [3]) hat vorgeschlagen, hier die Integral- 
gleichungsmethode zu verwenden. Auf der ERMETH (Elektronische Rechen- 
maschine an der ETH) wurde folgendes Beispiel gerechnet: : 
Gegeben sind zwei Standorte eines Himmelskérpers verschwindender Masse: 
t > V 
2. April 0,94299 0,47748 
1. Juni 0,35763 1,60295, 


wobei die Koordinaten in Gaufschen Einheiten gegeben sind. Gesucht ist der 
Standort am 22. April. Falls nur die Anziebung durch die Sonne beriicksichtigt 
wird, geniigt der Ortsvektor r = (x, y) der Differentialgleichung 
dr r ——— 
eae? is eae Sead foe 
wobei 
R= 2 9592 10-s7 


Das Zeitintervall vom 2. April bis 1. Juni wurde durch 7 aquidistante Stiitzstellen 
in 6 Teile eingeteilt. Man erhielt fiir die Koordinaten am 22. April 


# = 0,79285,, v= 091025. 


Es waren 4 Iterationsschritte nétig, wobei das erwahnte kubische Polynom als erste 
Approximation diente. Die exakten Werte betragen 


% = 0,79284, y= 091022. 


Es liegt auf der Hand, dass die beschriebene Methode auch angewendet werden 
kann, wenn in einer vorgeschriebenen Zeit ein gegebenes Ziel im Weltraum erreicht 
werden soll. 

Die Mitteilung ist ein kurzer Auszug aus einer Diplomarbeit fiir angewandte 
Mathematik, in der auch Vorschlage zur Verwendung der Gau8schen Quadratur 
gemacht wurden. 
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Summary 


A method is described for solving the boundary-value problem connected with 
the differential equation 
SS ee 
oe = te) 
The problem is replaced by an integral equation with piecewise linear Kernel. The 
integral equation is solved by successive approximation and numerical integration. 
Numerical results are given working on an example of celestial mechanics. 


(Eingegangen: 25. Juni 1960). 
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Varia — Miscellaneous — Divers 


Réunion de Printemps de la Société Suisse de Physique, 
a Winterthour les 7 et 8 Mai 1960 


Compte rendu des communications de physique et de mathématiques appliquées 


Uber das Impulsverhalten von Tunneldioden. Von E. BALDINGER und 
U. SpycHEr, Basel?) 


Anfangs 1958 hat Esaki [1]) als erster dariiber berichtet, dass bei hochdotierten ~ 


Flichendioden eine Kennlinie mit einem Gebiet negativen Widerstandes erreicht 
werden kann. Zur Erklarung der grundsatzlichen Wirkungsweise sei auf die Ar- 
beiten von Sommers [2], SKLAR [3] und Hatt [4] verwiesen. Stellt man einen 
Arbeitspunkt im Bereich des negativen Widerstandes ein (Figur 1), so lasst sich die 


Gr 
Oe i 
Ss ‘ 
4} i ae 
ea +40°C 
3 ii ee . 
es -16°% 
= an ee 
“S 
——" SS 
= 1 
rar 
0 1 1 n J 
0 100 200 300 400 
U|mV] —> 


Figur 1 


Statische Kennlinie einer Tunneldiode 


‘Tunneldiode als Verstarker, Oszillator oder Mischstufe bis ins Gebiet der kMHz 
verwenden [5, 6]. Fiir kleine Signale gilt dabei das Ersatzschema der Figur 2. 


R 


Figur 2 
Ersatzschema fiir kleine Signale. In abgestimmten Schaltungen ist in Serie zu r noch die Induk- 
tivitat von Zuleitung und Gehause zu beriicksichtigen (~ 2 nH) 


z Anstalt fiir angewandte Physik der Universitat 
5 5 . - . oo fi 
) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 521. 


“ 


— 


eee y 
Ne ee 
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Dieses Ersatzschema kann als Naherung fiir die ganze Diodenkennlinie verwendet 
werden, wobei allerdings R sehr stark und C schwach von der Spannung abhangig 
sind. Ist der Widerstand des steuernden Stromkreises grésser als Ro (iavadem klein 
sten Wert des negativen Widerstandes der Dioden-Kennlinie), so erhalt man zwei 
stabile Zustande 4 bzw. B (Figur 3). Messungen im Hinblick auf eine Anwendung 


100 200 300 100 
U|mV| —=— 


Figur 3 


Aussteuerung in Kippschaltungen. 


als Diskriminator ergeben bei konstanter Temperatur eine Reproduzierbarkeit des 
ersten Kipp-Punktes innerhalb von 0,2 °/5,. Wie man aus Figur 1 ersieht, besteht 
bei den ausgemessenen Dioden hingegen eine starke Temperaturabhangigkeit. 
Die Aussteuerung mit einem Stromimpuls, dessen Amplitude knapp tiber J, 
liegt (Lastgerade R, in Figur 3), ergibt einen Spannungsverlauf an der Diode nach 
Figur 4. Der zeitliche Verlauf dieser Spannung lasst sich qualitativ verstehen: nach 


t<~— 


Schaltung: 
R, K.0. 


a Aufgepragter Stromimpuls J = 7 mA, 
b Spannung tiber der Diode U = 410 mV, 
t = 0,1 us/cm. 


Figur 4 


einem Spannungssprung an y mit der Anstiegszeit des Triggersignals muss die 
Kapazitat auf die Kippspannung aufgeladen werden. Dann erfolgt das Umkippen 


mit der spannungsabhangigen Zeitkonstante RC. Die Riickflanke des Impulses 


ergibt zundchst wieder einen Spannungssprung liber 7. Anschliessend wird die 
Kapazitat entladen (die aufgespeicherten Ladungstrager verschwinden durch 
Rekombination, Diffusions- und Feldstrom). Die Tunneldiode verhalt sich vorerst 


wie eine normale p-n-Diode, deren Eigenschaften von HEINLEIN und anderen be- 


schrieben worden sind [7, 8]. Sobald das Gebiet negativen Widerstandes wieder 


_ erreicht ist, fallt die Spannung rasch auf Null zuruck. 
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Die Dauer des langsamen Aufladens von C ist durch die Stromamplitude gegeben. — 
Anstiegs- und Abfallzeiten sind somit in grossem Masse durch das Aufladen bzw. 
Entladen der Kapazitat im Gebiete positiven Widerstandes bestimmt, was insbe- 
sondere bei bistabilen und monostabilen Schaltungen zu beriicksichtigen ist. Die 
besten Resultate erreicht man durch gleichstrommassige Vorspannung der Diode 
iiber einen Lastwiderstand, der die Kennlinie méglichst nahe an den Kippunkten 
schneidet (FR, in Figur 3). . a 

Die Messungen wurden an einer Diode mit relativ grosser Zeitkonstante aus- — 
gefiihrt. Eine Diode mit hoher Grenzfrequenz wird grundsatzlich dieselben Er-) 
scheinungen zeigen. Will man mit Tunneldioden hohe Schaltgeschwindigkeiten 
erreichen, so ist in erster Linie den beschriebenen Effekten Aufmerksamkeit zu — 
schenken. ; 

Bestimmung der Parameter des Ersatzschemas (Figur 2). Die Bestimmung von 
y, R und C erfolgte in einer Briickenanordnung [9], verbunden mit einer ie 
leitung (Figur 5a). Unter der Voraussetzung, dass die Briicke ohne Tunneldiode 


Ry:R i 

p'=p+ 1-2 ‘ 
R,+R, 

; 

ie R ; 


a 


Figur 5 
a Briicke zur Bestimmung der Parameter; 6 Wechselstromersatzschema von Diode und Briicke. 


gleichstrommassig abgeglichen ist, lassen sich Gleichstrom und -spannung der 
Diode direkt aus Spannungsmessungen bestimmen. Es gilt 


U;= —U,, i,= U, >. 


Die Lecherleitung wirkt fiir alle in dieser Anordnung erzeugten Schwingungen als 


bf 
¥ 
H 
¥ 
. 
Induktivitaét. Figur 5b zeigt das Wechselstromersatzschema einschliesslich Diode. 


Im negativen Ast regt sich eine Schwingung w an, wenn 


/ 


Jt 
eee ers 
WS ard (Ce 


Fiir die Kreisfrequenz ergibt sich die Bestimmungsgleichung 


5 1 
or = = es — 2 = 
aye: (1 ined mie dey 
1+ w? R? a4 


Die Messung von o liefert anhand der beiden Beziehungen die Grésse C sowie eine 
Kontrolle der aus den statischen Messungen bestimmten Werte von 7 und R. 
Wir danken der Firma CLevite und Herrn LEUENBERGER fiir die Uberlassung 


einiger Tunneldioden und dem Schweizerischen Nationalfonds fiir die finanzielle 
Unterstiitzung dieser Untersuchungen. 
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Signalfluss und Analogiedarstellung der Neutronenédkonomie beim 
Abbrand von Reaktorbrennstoffen. Von W.HAtLG, Ziirich. Erschienen als 
Kurze Mitteilung, ZAMP 77, 433 (1960). 


Automatische Anlage zur Registrierung digitaler Messwerte. Von 
E. BALpDINGER, W. HAG und P. WEHRLI’) 


1. Arbeitsweise der Registrievanlage 


Es ist gegenwartig in vielen Zweigen der Forschung iiblich, fiir die Auswertung 
von Messdaten elektronische Digitalrechengerate zu verwenden. Eine solche Aus- 
wertung drangt sich besonders dann auf, wenn die Rechnungsprogramme hierfiir 
vielfaltig und zeitraubend werden oder wenn die Anzahl der Messdaten sehr gross 
ist. In allen diesen Fallen ist es wiinschenswert, die Messresultate in einer Form zu 
speichern, welche deren unmittelbare Ubermittlung in das Gedachtnis der Rechen- 
maschine erlaubt. Auf diese Weise sind Fehlerquellen, die bei der Ubertragung der 
Daten aus einem Messprotokoll in eine dem Rechner gefallige Form entstehen, 
vermeidbar. 

Die nachstehend in Figur 1 beschriebene Apparatur stellt eine Méglichkeit dar, 
diesen Anforderungen zu geniigen. Sie soll hauptsdchlich zur Registrierung kern- 
physikalischer Messergebnisse verwendet werden. Als erste Anwendungen sind die 
Aufzeichnungen der Messwerte von Halbwertszeitanalysen und Neutronendiffrak- | 
tionsexperimenten sowie von Neutronenflussverteilungen in einem unterkritischen 
System vorgesehen. Die Aufzeichnung geschieht mittels eines Blattschreibers, wie 
er im europdischen Fernschreibnetz verwendet wird. Angebaut an den Blatt- 
schreiber ist ein Empfangslocher, der gleichzeitig im 5-Lochcode einen Papier- 
streifen perforiert. Ein weiterer Lochstreifensender erméglicht die Kontrolle und 
eventuelle Duplizierung solcher Streifen. Wahrend der Schreibpausen bei automa- 
‘tischem Betrieb kénnen erginzende Informationen mit Hilfe der Handtastatur 
‘dem Messprotokoll beigefiigt werden. Der Blattschreiber bendtigt an seinem zwel- 
drahtigen Eingang fiir jedes zu schreibende Symbol*) eine Impulsfolge, die durch 
ein Umsetzgerat, den Zeichensender, aus den Gleichstromsymbolen der Mess- 
stellen erzeugt wird. Die Schreibgeschwindigkeit ist durch die Arbeitsweise des 
Zeichensenders, welcher wahrend der Erzeugung eines Symbols und nach Ablauf 
desselben fiir die Auslésung gesperrt bleibt, auf 5 Zeichen*) pro Sekunde beschrankt. 


8) Institut fiir angewandte Physik der Universitat Basel, Eidg. Institut fiir Reaktorforschung, 


Wiurenlingen. . 
a} Unter der Bezeichnung Symbol oder Zeichen werden sowohl Buchstaben als auch Ziffern 


‘ : A] sore T ii , 
‘sowie Hilfszeichen wie Zwischenraum, Zeilenvorschub, Wagenrticklauf usw. verstanden. 
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Ein am Ende eines Zeichens erzeugter Fortschaltungsimpuls kann zur Weiter-@ 
schaltung der Schrittschalter in der Haupt- oder Unterzentrale oder in der ange- 
schalteten MeBstelle beniitzt werden. ‘K an 

Die Hauptzentrale enthalt im wesentlichen neben einigen Hilfsrelais einen 
Sucher, der den Zeichensender entweder mit einer EinzelmeBstelle oder einer Unter- 
zentrale verbindet, sowie einen Schrittschalter. Ist eine MeBstelle oder eine Unter- | 


igs ee gs 


‘, 
M138 : A y 24h-Uhr R 


M21 niew2 


oS 


Py emacs el Ze: ice EE Haupizentrale 
Einzel- 
Messielle i 


p 
Piles Read-Out ; 
a Verstarker : 
Hana HOES tes . 
(EES cosas Pe 


Zeichen- Fernschr| 
Sender +F-Locher 7 


Figur 1 
Blockdiagramm der Anlage. Mit M; , sind die verschiedenen Mef stellen der 3 eingezeichneten — 
Experimente bezeichnet. 


zentrale ablesebereit, so setzt sich, durch ein entsprechendes Signal betatigt, der 
Sucher so lange in Bewegung, bis die MeBstelle gefunden ist, worauf der Antriebs- 
motor des Fernschreibers anlauft. Hierauf wird mit einiger Verzégerung ein fest 
eingestelltes Ableseprogramm mit Hilfe eines Schrittschalters gestartet. Dieses 
Programm besteht vorerst aus den fiir Wagenriicklauf, Zeilenvorschub und Um- 
stellung auf Ziffernempfang notwendigen Zeichen, worauf eine Kennzeichnung der 
angeschalteten Mefstelle erfolgt. Da es oft wiinschenswert ist, die Zeit der Auf- 
zeichnung festzuhalten, wird anschliessend eine mit der Hauptzentrale verbundene: 
24-Stunden-Uhr abgelesen, bevor die Daten der betreffenden MeBstelle heraus- 


in einem Transistor-Leistungsverstarker (Read-out-Verstarker) soweit verstarkt, 
dass ein einwandfreies Ansprechen der Relais im Zeichensender gewahrleistet wird 
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Diese Schaltung erméglicht eine Reduktion der Anzahl von Leistungsverstarkern 
auf insgesamt 10, entsprechend den Ziffern 0 bis 9. Die Ausgange der Mef stellen 
erlauben den Anschluss eines bis 7u 300m langen Ve . 
Unterzentrale und Hauptzentrale. 

Im Vollausbau kénnen 10 Abonnenten an die Hauptzentrale angeschlossen 
werden. Die gegenwirtig erstellten Unterzentralen sind fiir 6 MeBstellen ausgelegt. 
Unter der Annahme, dass bei jeder dieser MeRstellen eine 6stellige Zahl abgelesen 
werden soll, erfordert die Registrierung aller Informationen einer 6fachen MeBstelle, 
einschliesslich aller Hilfssymbole wie Zwischenraume, Wagenriicklauf und Zeilen- 
vorschub, ferner einer allfalligen Verzégerung, verursacht durch die Schaltung der 
Minutenuhr, 15 Sekunden‘). 

Ein Teil der MeBstellen wird in der Regel aus einem Geber (zum Beispiel Parti- 
kelzahler) mit nachgeschaltetem Untersetzer bestehen, Gegenwiirtig sind 4 fiir den 
Exponentialversuch und 3 weitere zur Halbwertszeitanalyse im Betrieb. Die Steue- 
rung der Experimente geschieht durch Impuls2 eines Programmgebers mit einer 
Zeitverzégerung von einigen Mikrosekunden. 


rbindungskabels zwischen 


2. Beschreibung der Untersetzer 


Um die Warmeproduktion und das Volumen der hier beschriebenen automati- 
schen Registrieranlage in bescheidenen Grenzen zu halten, wurden zum Bau der 
elektronischen Schaltungen ausschliesslich Transistoren verwendet, wobei die iiber- 
wiegende Zahl auf die Untersetzereinheiten entfallen. Wie aus dem Blockdiagramm 
Figur 2 zu ersehen ist, besteht jeder Untersetzer aus der Eingangsschaltung, welche 


vom Kdo.-US. von 


Unierzentrale 


Sieverimpuls 9 7uIP 


. Sfeuersignal von 
fur weitere US. 


| Unterzentrale 


6 Dekaden 
—»— — Fndverstark|->— Unferzentrale 


Blockdiagramm eines Untersetzers. Kdo-US = Kommando-Untersetzer. 


Schritt- 
schalter 


Read- out 


pingange ” Verstérker 


Vorwan/ 


Figur 2 


die fiir den automatischen Betrieb notwendigen Hilfsschaltungen zur Steuerung 
enthalt, fe-ner aus 6 Dekaden, Endverstarker, Schrittwahler und einem stabilisier- 
ten Transistor-Netzgerat. Der Schrittwahler dient zum Herauslesen der in den 


Dekaden gespeicherten Zahlen. Die Ablesung erfolgt auf ein Kommandosignal der 
- Zentrale. Gleichzeitig stellt der Schrittwahler die Nullstellung der Untersetzer her 


- 


_s 


und liefert beim Zuriicklaufen in die Ruhestellung ein Hilfssignal, das die AblL sung 
der nachsten Einheit einleitet. Zur Kennzeichnung einer Untersetzer-Einheit 


5) Fir den Weiterausbau der Anlage ist die Verwendung eines Blattschreibers oder Streifen- 
lochers mit Paralleleingang vorgesehen. Dadurch wird die Registriergeschwindigkeit auf 10 Zeichen 
resp. 33 Zeichen pro Sekunde erhoht. 
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modgen die folgenden Angaben dienen: 

Auflésungsvermoégen: 0,9 us 

Maximale Umgebungstemperatur: 40°C 

Zulassige Schwankung des Kraftnetzes: 190 bis 240 V 

Leistungsaufnahme: 13 Watt 

Eingangsimpedanz: 330 Ohm 

Eingangssignale: + oder — (umschaltbar) mit minimal 1 V,,, 

«Preset time» oder «Preset count» einstellbar. : 


Limiter Schmilt-Trigger Gate 


zum Abschalten 
weiterer US. 


Vereinfachtes Schema der Eingangsschaltung. 


. 
t 
{ 
. 
Figur 3 } 
1 
: 
: 


Im fol : er Ree 
oe. a ee wird die einzelne Dekade und die Eingangsschaltung etwas naher 


a) Dekade: 
y Jede Dekade besteht aus 4 dualen Stufen mit 2 Riickfithrungen zur Unterdriik- 
ee ee eee Stellen (Code: 1, 2,4, 2). Um mit Transistoren der Type 
ein Autlosungsvermogen von iiber 10-6 i u i 
storen nicht tibersteuert warden AN ES 
aN aha den einzelnen dualen Stufen und den beiden Riickfithrungen sind zur_ 
ee apres: verwendet, welche eine wesentliche Verbesserung der 
eee me = iene os Temperaturabhangigkeit gestatten. Die visuelle Anzeige erfolgt 
alee ae ce ype Philips DM1 60 im erwahnten Code. Die Umrechnung von 
, 2, 4, 2 in die Ziffern 0 bis 9 wird in einer fiir alle Dekaden und alle Unter- 


setzer gemeinsamen Einheit vorgenomm Tamm E 1gur 1 im «Read 


* 
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b) Eingangsschaltung 

Mit Hilfe eines von einem Steuergerat gelieferten Dauersignales lasst sich der 
Untersetzer ein- bzw. ausschalten. Nach Wunsch kann das Ausschalten ebenfalls 
durch einen kurzzeitigen Impuls geschehen, zum Beispiel dadurch, dass einer der 
Untersetzer (Kommando-Untersetzer) eine einstellbare und vorgewahlte Ziffer 
erreicht hat. Alle Untersetzer sind gleichartig und austauschbar und durch geeignete 
Wahl der Verbindungen kann irgendeiner zum « Kommando-Untersetzer» bestimmt 
werden. 

Die zu zahlenden Impulse werden zuniachst im Schmitt-Trigger normiert und 
gelangen tiber das «Gate» zum Eingang der 1. Dekade. Die Verwendung von Tran- 
sistoren der Type OC170 gewahrleistet auch bei us-Impulsen ein einwandfreies 
Arbeiten der Schmitt-Trigger-Schaltung. 

Abschliessend sei noch bemerkt, dass sich eine solche transistorisierte Anlage 
leicht von einem kleinen Akkumulator im Pufferbetrieb speisen ldsst, wodurch 
Stérungen durch Unterbriiche des Kraftnetzes vermieden werden. 

Ein Teil der Entwicklungskosten wurde aus Mitteln der Kommission fiir Atom- 
wissenschaft des Bundes finanziert, und wir méchten es nicht unterlassen, fiir diese 
Hilfe unseren besten Dank auszusprechen. 


An Experimental Investigation on Potential Flow through Spiral 
Casings. By N. ArHaNassIADIs, Ziirich §) 

The flow through a water turbine spiral casing, outside a cylindrical surface Sy 
surrounding the guide vanes, can be considered as a potential one, if: 

a) the fluid is incompressible and frictionless, 

b) the flow is continuous and stationary, 

_c) there is not any vorticity in the inlet pipe and 
d) we eliminate the fixed blades outside Sp. 

It is known from potential theory, that, in a region, an harmonic function e.g. 
the velocity potential ¢, is determined if c, = 0¢/0n throughout its surrounding 
surface is determined; x is the normal direction to this surface. Therefore, assuming 
that the inlet pipe of the casing is straight and of a constant cross-section to infinity, 
the potential flow through the casing is determined if the radial component C,, of 
the velocity along S, is determined. To find the function ¢ (C = grad ¢), through-. 
out the casing, we have to solve an integral equation. The problem is not easily 
treated by analysis, especially in a three-dimensional casing ; thus, in ourinvestigation 
we used the electric analogy method. 

Our experimental device is shown in Figure 1; the conductor was tap water. To 
identify the boundary conditions between ¢ and the electric potential V, we 
replaced a cross-section in the inlet pipe with an electrode (P-electrode) and re- 
presented the surface S, by a collector of electrodes. The current J, through each 
one of the electrodes should be regulated in such a way that Q, ~ o J, (k = 1, 
2,..., N), where Q,, is the volume flow-rate passing on at the K-Electrode, and o 
the specific resistance of the conductor. By 4 or 5 iterations of this procedure, we 
obtain the desired distribution of J, along Q,; thus. we establish in the model an 
electric field corresponding to ¢. In other words, we have solved the appropriate 
integral equation. To avoid any secondary current between successive electrodes, 
and to minimize as much as possible any disturbance of the electric field near 
S,, the electrodes are placed inside Sy as shown in Figure 2. 


6) ETH, Institut fiir Aerodynamik. 
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P-Electrode 


Figure 1 


Electric analogy device 
1 Oscillator (regulated at 300 Hz); 2 Transformer (1:1); 3 Electric Analogy Model (two- or three- 
dimensional); 4 Potentiometers to regulate the current through the electrodes; 5 Electronic Volt- 


meter (for the above operation); 6 Potential pick-up (probe); 7 Zero-indicator unit (Cathode-ray — 


oscilloscope); 8 Decade potentiometer; 9 Pick-up transformer and filter; 10 Condensers for a 
complete balancing of Wheatstone-bridge. 


Figure 2 


From our experiments, we obtained a remarkable result valid for two- as well 
as for three-dimensional casing forms as follows. If Cy, = const. along So, or 
has not extremely large variations along it, the distribution of the moment of 


: 
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momentum per unit mass of the fluid along’a meridian section of the casing 
remains constant, i.e. in cylindrical coordinates (el zie 


4° Vlg = Cz, Yolo = const (1) 


with 6 = const. (9 < 22) and independently of z; c, is the tangential component 
of the velocity. Thus, if the casing form and c,, along Sy are given, we can calculate 
¢,,(9) and consequently c, for every 0 < 22. Equation (1) is well known for a 
two-dimensional pure logarithmic field, where both c,°¥ = const. and ¢,*7 = const. 
throughout. We found here that the same is true for every spiral casing form, but 
only on a meridian section and for @ < 22; in general c,-v changes from section 
to section. 

It has been experimentally verified that the runner of a turbine as well as the 
guide vanes act as rectifiers of the flow through them; thus, if the casing is sym- 
metrical to a plane perpendicular to the z-axis (axis of the runner), ¢,, = const. 
along So. Figure 3 shows the (c,)-distribution in such a casing for c,, = const. The 
solid (¢,),,-curves are drawn according to equation (1); we can see that the measuring 
points verify this equation up to @ ® 320°. 


WL Plane of Measurements 


0 a/2 It 3n/2 2X 
Figure 3 


Conclusion: By the experimentally determined equation (1), we can design 
spiral casing forms to obtain c,, = ¢,, tg% along SG where % is given and ¢,,, 
c, = const. It is also possible to modify existing casing forms and fixed blade 
= : 
shapes outside So. 
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Zur Theorie der Luftkissen-Fahrzeuge. Von H. BAUMANN, eee ’ 
Ein Luftkissen-Fahrzeug wird von einem Kissen getragen, das en pee er 
erhéhtem Druck stehende Luft gebildet wird. Dieses kann auf verschiec ene Arten 
erzeugt werden, von denen hier nur auf die Methode des Luftvorhanges meee 
werden soll. — Ein langs der Peripherie aus dem Boden austretender, run pa 
geschlossener Luftvorhang halt im Raum unter dem Fahrzeug einen Uberdruc 
ne aoe im Vorhang kann als ebenes Problem behandelt werden (ebener 
Strahl). Mit Hilfe des Impulssatzes lasst sich dann ein Zusammenhang zwischen 
Druckdifferenz am Vorhang und Schwebehohe naherungsweise angeben, wenn 
Dicke, Geschwindigkeit und Winkel des Vorhangs beim Austritt aus dem Fahrzeug 
bekannt sind. Andererseits ist es auch méglich, das genaue theoretische Berech- 
nungsverfahren fiir freie Strahlen hier anzuwenden (Hodograph, konforme Abbil- 
dung). Die Ubereinstimmung der Ergebnisse der beiden Rechnungsarten ist gut, 
und auch der Vergleich mit Messungen zeigt, dass die einfache Anwendung ae 
Impulssatzes zu Resultaten fiihrt, die der Wirklichkeit recht genau entsprechen‘). 
Veranderte Bedingungen treten auf, wenn ein Luftkissen-Fahrzeug sich in Vor- 
wartsfahrt befindet. Auf der Aussenseite des langs der Vorderfront verlaufenden 
Luftvorhangs wirkt dann zusatzlich der Fahrtstaudruck. Daher wird die durch 
diesen Vorhang aufzubringende Druckdifferenz vermindert, was eine Einsparung 
von fiir den Auftrieb erforderlicher Luftmenge und Leistung bedeutet. 
Betrachtet man ein rechteckiges Fahrzeug, das ausschliesslich iiber Wasser- 
flachen verkehren soll, so kann man sich die beiden seitlichen Luftvorhange zum 
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Abhangigkeit der Dicke s, des vorderen Vorhanges, der Luftmenge Q und der theoretischen Leistung 

L, fiir den Auftrieb von der Vorwartsgeschwindigkeit V. Das angenommene Luftkissen-Fahrzeug hat 

folgende Daten: Lange L = 25 m, Breite B = 8 m, Dicke des hinteren Luftvorhanges Sy = 20 cm, 

Austrittswinkel m = 45°, Masse M = 20000 kg, Ansaugung der Luft von oben. Die Schwebehiéhe h 
betragt 1 m. 


*) Institut fiir Aerodynamik, ETH. 


8) J. AcKERET und H. Baumann, Uber Luftkissen-Fahrzeuge. Schweiz. Bauztg., 78. Jg. 
137-141 (1960). 
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Beispiel durch Bleche ersetzt denken, die wenig ins Wasser eintauchen. Diese 
k6nnen sehr diinn sein und ergeben keinen grossen Widerstand. Infolge der gerin- 
geren Druckdifferenz am vorderen Vorhang lasst sich entweder bei gleichbleibender 
Vorhangdicke dessen Geschwindigkeit verkleinern oder bei gleichbleibender Ge- 
schwindigkeit kann er diinner gemacht werden. Sobald die Vorwartsgeschwindigkeit 
so gross wird, dass der Staudruck der Flachenbelastung des Fahrzeuges entspricht, 
ist der vordere Vorhang iiberhaupt iiberfliissig. Leistung und Luftmenge sinken 
in diesem Fall auf die Halfte gegeniiber den Werten im Stillstand. In der Figur sind 
diese Gréssen sowie die Dicke des vorderen Vorhanges fiir ein angenommenes Fahr- 
zeug in Funktion der Vorwartsgeschwindigkeit aufgetragen. 

Bei diesem Beispiel erfolgt der Lufteintritt in das System senkrecht von oben. — 
Eine weitere Leistungsverminderung fiir den Auftrieb erreicht man, wenn auch hier 
der Fahrtstaudruck ausgeniitzt wird, indem die Lufteintrittséffnung zum Geblise 
nach vorn gerichtet wird. Die Umlenkung um 90° nach unten hat aber einen zu- 
satzlichen Widerstand zur Folge, der durch vergrésserten Schub iiberwunden wer- 
den muss. 


Laminare Strémung in Rohren mit verschiedenen Querschnitten. 
Von V. STINGELIN, Ziirich’). 

Im Gegensatz zur turbulenten Rohrreibung, wo erfahrungsgemdss der Druck- 
abfallkoeffizient 4 fiir hydraulisch glatte Rohre nur von der Reynoldszahl abhangt, 
sofern diese auf den hydraulischen Durchmesser bezogen ist, ergibt die strenge 
Theorie fiir die laminare Strémung eine zusdtzliche Abhangigkeit von der Geo- 
metrie der Querschnittsform. 

In der Literatur wird diesem Gegensatz wenig Beachtung geschenkt. Entgegen 
den Erwartungen fallen die Messpunkte der Experimente von NIKURADSE fiir ver- 
schiedene Profile auch im laminaren Gebiet auf die theoretische Kurve fiir kreis- 
runde Rohre. SCHILLER hat schon auf diese Tatsache hingewiesen; doch sind die 
experimentellen Nachweise wenig iiberzeugend. 


Asympt.f Rechteck 


Asympt. t EIipse 


Ses i 


100 


90 


(oe Ce ie Geeta 9 S10 ait. 2 
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Figur 1 
Laminares Rohrreibungsgesetz A(Re) fiir verschiedene Querschnitte. 


%) Institut fiir Aerodynamik, ETH. 
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Das Ziel der nachfolgenden Arbeit war deshalb die experimentelle Bestatigung 
der schon bekannten exakten Lésungen der laminaren Stromung durch verschiedene 


Querschnittsformen. 
Das laminare Druckabfallsgesetz lautet: 
adele 
a Ren 


wobei f eine Konstante ist, die von der Geometrie des Querschnittes abhangt. 

In Figur 1 sind die theoretischen Werte von f in Abhangigkeit der Halbachsen- 
verhiltnisse der Ellipsen, der Seitenverhaltnisse der Rechtecke, und vergleichsweise 
die Werte des Kreisrohres und des gleichseitigen Dreieckes aufgetragen. Selbst bei 
kleinen Anderungen der Achsen- bzw. Seitenverhaltnisse werden die Abweichungen 
vom Kreisrohr schon bemerkenswert, daher war ein deutlicher experimenteller 
Nachweis zu erwarten. : 

Die Ubereinstimmung der Messungen mit der Theorie Figur 2 war dann auch 
gut, besonders beim kreisrunden und dem speziell konstruierten 1/10 Rechteckrohr 
aus Messing. Kleinere systematische Abweichungen der mittleren drei Rohre waren 
auf Ungenauigkeiten in der Geometrie zuriickzufiihren und zeigten die vermutete 
Empfindlichkeit der Zahl 6 von der Kontur des Querschnittes. 


ey ia SCO a, (GIT) (CME) 
11g 2° 0 04 01257 
s oM/; 0414 0159 
0 OM/2 069% 0489 
‘1 9521 0,439 
0366 0,4033 


8 
6 
fs f =Querschnittsriache 
U=benefzter Umfang 
9.1972 | | | | 
200 400 600 
U-d, 
Re= Y A ————S=_ 
Figur 2 


Laminarer Rohrreibungskoeffizient A fiir Rohre mit verschiedenen Querschnitten. 


Auswertung von Warmeleitfahigkeitsexperimenten mit Hilfe eines 
elektronischen Rechenautomaten. Von P. Wyner, Zurich), 

Bei Warmeleitfahigkeitsexperimenten im Gebiete tiefer Temperaturen haben 
Kohlewiderstandsthermometer gegeniiber den wesentlich unhandlicheren Gas- 
thermometern den grossen Nachteil, bei jedem Experiment mit Hilfe der Helium- 
dampfdruckskala neu geeicht werden zu mussen; sie sind nur innerhalb des Tief- 


10) Institut fiir kalorische Apparate und Kaltetechnik, ETH. 
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temperaturgebietes geniigend genau reproduzierbar. Warmt man sie bis zur 
Zimmertemperatur auf und kiihlt sie wieder ab bis zur Temperatur des fliissigen 
Heliums, so andert sich ihre Charakteristik um einige Prozente. 

> In dem uns interessierenden Temperaturgebiet von 1,5°K bis 4°K haben sich 
kaufliche Radiowiderstande von !/,, Watt Leistung der Allen-Bradley Company [1]}") 
als besonders geeignet erwiesen. Ebenfalls sehr giinstig sind nach DE VROOMEN 
hergestellte Widerstande aus einer Kohlepulver-Aquadag-Emulsion, die auf einen 


isolierten Kupferdraht aufgetragen wird [2]. Charakteristische Widerstandswerte 
sind in Tabelle 1 zusammengestellt. 


Tabelle 1 
ALLEN-BRADLEY DE VROOMEN 
Zimmertemperatur 100 Q 1000 Q 
4,2°K w= 2000 Q fw 2600 Q 
2g OS = 10000 Q ew 3600 Q 


Leider ist der Leitfahigkeitsmechanismus in diesen Kohleschichten noch nicht 
geniigend geklart. So gibt es nur eine zweikonstantige semiempirische Beziehung [3], 
die allerdings das Verhalten sowohl der Allen-Bradley- wie auch der de Vroomen- 
Widerstande sehr genau beschreibt: 


In R 


TR Nate Hie 


Dabei bedeutet T die Temperatur und R den Widerstand; a und 6 sind an- 
zupassende Konstanten. Wie man aus dieser komplizierten Formel sieht, sind 
genaue Auswertungen von Experimenten recht mtihsam und zeitraubend. DE VRoo- 
MEN [4] hat kiirzlich eine Methode angegeben, die wesentliche Zeitersparnisse mit 
sich bringt; trotzdem ist der Zeitaufwand fiir die Auswertuug immer noch sehr gross. 

Macht man verhdltnismassig viele Warmeleitfahigkeitsmessungen, so miissen 
fiir die Auswertung stets die gleichen Rechnungen wiederholt werden. Dies ent- 
spricht genau der einem elektronischen Rechenautomaten angebrachten Problem- 
stellung. Wir haben daher die Auswertung unserer Experimente fiir die ERMETH 
(elektronische Rechenmaschine der ETH) programmiert. 

Dem Automaten werden die Eich- und Messwerte, die Probenabmessungen und 
die zugefiihrte Heizleistung eingegeben, wahrend er direkt den Wert fiir die Warme- 
leitfahigkeit herausdruckt. Die Bestimmung der Temperaturdifferenz ist dabei der 
zentrale Teil des Problems. Zur Lésung kann man entweder den Eichpunkten 
eine Funktion einpassen, zum Beispiel mit Hilfe von Tschebyscheff-Polynomen, 
oder man beniitzt die angegebene Gleichung und bestimmt mit Hilfe einer Aus- 
gleichsrechnung die beiden Konstanten a und b. Wir wahlten den zweiten Weg, da 
er gewisse Extrapolationsméglichkeiten offen lasst. 

Beniitzt man die in der Ausgleichsrechnung iibliche Symbolik 


n 
(x, ¥] =x; Vi» 
ia 


11) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 534. 
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so erhalt man mit den Gleichungen 


In R,\1/2 
a+oblnR, ( T 4) =7; 
und 

In RF, \1/2 
A, = 10 Rye 2b = (=) 
fiir die beiden Konstanten die Werte 


[1, B] [4, A] — [A, B] (1, A] (4, Bit, =, By 
fief, (ea et ee ~ Tt, 1) (4, 4) = ia 


Zu Kontrollzwecken wird die Summe [y, 7] berechnet, othe bei einer guten Ausgleichs- 
rechnung moglichst klein sein soll. 

Die eigentliche Rechnung geht von der Voraussetzung aus, dass die Mess- und 
Eichwerte mittels Lochkarten gemass Tabelle 2 bereits in der Maschine gespeichert 
sind. Der Probendurchmesser d und der Sondenabstand Ax wird wahrend der 
Rechnung von Hand eingegeben. 


Figur 1 
Flussdiagramm fiir Ausgleichsrechnung. 
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Tabelle 2 
Bereits gespeichert: 
Physikalische Bedeutung Symbol Blockadresse 
Be ee ee a 
Eichwert Widerstand1... . RE . 9,0200 
Eichwert Widerstand2.... RS 4 9,0300 
Eichwert Temperatur ... . me 9,0400 
Messwert Widerstand1. . . . rit 9,0500 
Messwert Widerstand2... . ee 9,0600 
Messwert Heizung Strom . . . Je 9,0700 
Messwert Heizung Spannung . Vi 9,0800 


Berechnung des Warmewiderstandes 
1, Ax, Drucken 


k= 


ny) 


n($)? 
Ax 


JR6=>500 


in Re es In? 


MT 
6 


_ 


Drucken 1; 


Tits 72 
Is = 23, Drucken 


Al, =Tg"- Tf", Drucken 


W,=k a Drucken 


(a'+b' nrg ® = (924 b? ink) 
Las 


+200 YM5+300 


nein 
Tg’ = Q-Zanl 2 son Tgss0o= O-Zahl? > Nowz0d~ ZANT Wi Tal? 
ja ja ja Ja 
Teveo~O-Zanl? [={fin7) [find 
Ja 
Js-an0~4-Zahl? pea LE 
Jd 
=Q-Zahl2 P—Fing 
Vas sq O-Zanl? | 
Ja 
Figur 2 


Flussdiagramm fiir Warmewiderstandsbestimmung. 


qn 
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Die letzte Zahl jedes Blockes ist mit einem Q-Zeichen versehen, Das gibt eine 
gewisse Kontrolle dafiir, dass sich zueinander gehérende Widerstands- and Sones 
ratur-, oder Widerstands- und Leistungswerte auch entsprechen. In den Figuren 
und 2 sind die Flussdiagramme des Programms dargestellt. Der Mangel an Zyklen 
wird dadurch aufgehoben, dass durch die vielfache Verwendung das Programm 
«aussen herum zyklisch» wird. 


Twatr | scheinbar 


005, 
°K 
004 | 


Awe Fichpunkte 
Ac 0 Vergleichspunkte 


003 | 
002 | 


01 | 


000 
- 001 | 
-g02 | 
-903 | 
04. 


-005 | 


-006 * 
Figur 3 


Abweichungen in den nach der Ausgleichsrechnung bestimmten Temperaturen. 


Durch die Ausgleichsrechnung werden gewisse Abweichungen von den wahren 
Temperaturen induziert. Eicht man zwischen 2° und 4°K ungefahr jedes 4/,, Grad 
und verwendet alle Eichpunkte fiir die Ausgleichsrechnung, so ist diese induzierte 
Abweichung im Maximum 0,03°. Dies ist fiir unsere Experimente zu gross. Teilt 
man das ganze Temperaturintervall in zwei Gruppen von 4° bis 3° und von 3° bis 
2° ein und wahlt pro Gruppe je vier Eichpunkte, so kann man mit dieser Gruppen- 
einteilung erreichen, dass der induzierte Fehler im ganzen Bereich nie grésser als 
0,002° wird. Diese Abweichungen sind in Figur 3 dargestellt. 


Diese Arbeit wurde durch einen Arbeitsbeschaffungskredit des Bundes finanziell 
unterstitzt. 
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Vergleich hydraulischer | digitaler Schaltelemente mit elektronischen 
und elektromechanischen «Aquivalenten». Von H. GLArrrr und H. R. Mit- 
LER, Ziirich!2), ae 

Hydraulische logische Elemente unterscheiden sich von elektronischen und 
elektromechanischen Komponenten in mehrfacher Hinsicht: 

a) durch die mittels eines Elementes darstellbaren Funktionen, 

b) durch die bei Benutzung der verschiedenen Eingiange entstehenden Impedanz- 
verhdaltnisse, und 

c) durch die im Zusammenhang mit Verzweigungen entstehende Zunahme der 

Ansprechzeit. 

Figur 1 zeigt einen NPN-Transistor in Kollektorbasis-Schaltung und ein ein- 
faches Ventil sowie die dazugehérenden Schaltfunktionen. 


G 
m 
B 
agi 
‘ i 
B 
X=A+BC X=AB+BC 
Figur 1 


Vergleich Transistor - Ventil. 


In bezug auf das logische Verhalten sind beide Elemente als Bausteine fiir belie- 
bige Matrizen verwendbar, da in beiden Fallen ohne Beizug weiterer Komponenten 
Matrizenpunkte auf eine Ausgangsleitung zusammengefasst werden kénnen. 

Das Bild dndert sich aber, sobald man die Impedanzverhaltnisse in Betracht 
zieht. Das obenerwahnte Zusammenfassen besteht in einer Substitution 


A RG AX ome Ny By Cem ee Kg 3 Ba Go 


usw., wobei der damit verbundene Anstieg der Ausgangsimpedanz in hydraulischen 
Schaltungen gering bleibt, bei Verwendung von Transistoren aber sofort prohibitiv 
wirkt. 

Etwas ungiinstiger fiir das Ventil fallt der Vergleich mit einem Relais aus. Die 
Tatsache, dass, im Gegensatz zum Relais, ein einziges Ventil die «exclusive-or»- 
Funktion nicht darstellen kann, ist weniger gravierend als die im Zusammenhang 
mit Verzweigungen auftretende Verlangerung der Ansprechzeit bei hydraulischen 
Schaltmitteln. 

Besonders augenfallig wird dieser Effekt am Beispiel des Paralleladdierwerks, 
wo das Vorhandensein éiner mehrfachen Verzweigung ohnehin nicht auf den ersten 
Blick evident ist. 

Im Prinzip lasst sich eine hydraulische Schaltung analog zu bestimmten Relais- 
Schaltungen aufbauen, bei welchen das Ubertragssignal die Lage der Kontakte 

‘nicht dndert, sondern sich iiber vorbereitete Wege ausbreitet. Es geht nun aber 


12) [BM-Forschungslaboratorium. 
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nicht an, die Additionszeit aus der zum Stellen der Ventile notigen Zeit sowie aus 
Geometrie und Schallgeschwindigkeit zu bestimmen. Der Grund liegt darin, dass 


zum Zwecke der Weiterverarbeitung oder Anzeige der Information in den verschie- | 


denen Stufen endliche Volumina verdrangt werden miussen. 


Figur 2 gibt die Zeit zur Ausbreitung des Ubertragssignals in Funktion der — 


Anzahl Stufen fiir zwei verschiedene Viskositaten wieder : 


ay 
(msec) 8 7 


2+ 
iyee 
+++ ——— 
i ey Ga ts te) SW ail Se 
Figur 2 


Ausbreitungszeit des Ubertragungssignals in Funktion der Anzahl Stufen eines hydraulischen ~ 


Addierwerkes. 


Die beiden hier erwahnten Eigenschaften hydraulischer Komponenten und 
Schaltungen unterstreichen die Notwendigkeit, spezifische Schaltungen, speziell 
unter Verwendung der Matrizentechnik, zu entwickeln. 


Méthode de mesure de la stabilité 4 court terme de la fréquence 


d’oscillateurs 4 quartz. Par P. Karrascuorr!’), J. DE Prins!) et J. Bona- 
NOMI}4), 


1. Introduction 


Le but des mesures présentées dans cette communication est de déterminer la 
précision avec laquelle la fréquence des oscillateurs 4 quartz peut étre mesurée pour 
une durée de mesure donnée. Nous nous sommes surtout intéressés a la stabilité 
a court terme d’oscillateurs 4 quartz servant a la comparaison des horloges a quartz 
aux €talons de fréquence atomiques. En effet, seules les variations lentes pouvant 
étre corrigées a l’aide de servomécanismes, les oscillateurs employés doivent étre 
particuliérement stables. D’autre part, dans un résonateur a césium, les fluctuations 
de fréquence du signal excitant peuvent engendrer un élargissement de la raie de 
resonance observée. [1]1*). De méme, la précision de la comparaison de fréquence 
entre un oscillateur 4 quartz et un maser & NH, est limitée par les fluctuations 
rapides de la fréquence de l’oscillateur 4 quartz [2, 3, 4]. 

Nous avons mesuré ces fluctuations en faisant des séries de mesures de fréquence 
de courte durée. Ensuite, nous avons déterminé l’écart type o, de la fréquence 
mesurée pour des durées de mesure 7 allant de 10 A 100 ms. Comme les fréquences 


13) Laboratoire suisse de recherche horlogéres, Neuchatel. 
14) Observatoire de Neuchatel. 


1%) Les chiffres entre crochets renvoient A la Bibliographie, page 539. 
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sont comparées dans le domaine des microondes, la sensibilité du dispositif de 
mesure est élevée. 


2. Méthode de mesure 


La figure 1 montre le schéma de principe de l’appareillage employé. Les deux 
oscillateurs comparés sont identiques et oscillent a des fréquences de 8,5 MHz. Ces 
fréquences, multipliées par 2700 A l’aide de deux chaines de multiplication de 
fréquence, produisent deux signaux de 22 900 MHz dans le domaine microonde. Un 
récepteur superhétérodyne employant un klystron comme oscillateur local sert A la 
détection du battement de ces deux signaux. 

Le battement, dont la fréquence peut étre variée entre 0 et 50 Hz par un réglage 
approprié des oscillateurs 4 quartz est obtenu aprés démodulation du signal de 
moyenne fréquence. 

Les résultats des mesures sont obtenus de la maniére suivante: Les durées 1, des 
alternances positives du battement sont représentées par la longueur des impulsions 
sur l’écran d’un oscillographe cathodique et enregistrées par photographie. La 
hauteur de chaque impulsion est proportionelle 4 la valeur t; correspondante. La 
formation du train d’impulsions est illustrée par la figure 1. 


Addition 


Géneérareur 
Gein 
dent de scie 


vers Oscillographe 


22900MHZ 
MF 30MHZ 


Schmitt 
Trigger 


Figure 1 


Le nombre » de mesure par série est d’environ 150. Remarquons ici que les 
fluctuations t; sont petites par rapport a la valeur moyenne 


a ee) Se eo 
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ou v= 22900 MHz et Ay/r est l’écart relatif de la fréquence par rapport a sa 
valeur moyenne, pour un 7; donné. Ainsi pour tT = 10 ms, un écart de t, de 1% 
donne 4,v/v) = 2: 10-11. La hauteur des impulsions enregistrées a été mesurée a 


0,5% pres. 
3. Résultats 


Le dépouillement des enregistrements photographiques donne pour chaque 
série de résultats un écart-type des Tt; 


7 
| 
: 
: 
| 


avec 


Ensuite, on obtient pour l’ensemble des fluctuations relatives de fréquence des 
deux oscillateurs un écart type 
Vp il 


Oo oe Oh i == Se 
y Vo 14 b 2T 


L’écart type o, pour un seul oscillateur est donné par 


ceci en admettant que les fluctuations respectives sont décrites par deux distri- 
butions normales identiques. 

La figure 2 montre les résultats obtenus pour des valeurs de 7 allant de 10 a 
100 ms. 


10 x10" 
{ ait 
Bt ia is a 
| | i 
6 ~ | = 
5 — | 
RS 
4 aac r + 
re i 
3 { a le 1 = 
PS 
i const-T 2 
2 L | Ede | 
+ SQ 
a 
> eee 
- 20 30 40 50 700 200 msec 


Figure 2 
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4. Conclusion 


Dans le Soman’ Sheets, les valeurs de o, peuvent étre considérées comme 
proportionnelles a7r*,cela en tenant compte de l’imprécision de o, qui est estimé 
a 2 +1071, 

Il est, de plus, 4 remarquer que les fluctuations pouvant étre dues au bruit de 
fond dans les chaines de multiplication de fréquence sont négligeables. 
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Ein Ultrahochvakuum-System mit einer ausheizbaren H¢ -Diffusions- 
pumpe zur Erzeugung von Vakua im Gebiet von 10-!! Torr. Von E. B. Bas, 
Ziirich!*), 

In den UHV-Systemen mit Hg-Diffusionspumpen kann ein Endvakuum von 
etwa 2 X 107! Torr. erreicht werden. Dieses Endvakuum ist scheinbar in der 
Hauptsache durch die Gasabgabe der Diffusionspumpe selbst bedingt. Um noch 
Jkleinere Drucke erreichen zu kénnen, haben wir ein UHV-System entwickelt, bei 
dem eine einstufige Hg-Diffusionspumpe zusammen mit dem Rezipienten und 
den Kiuhlfallen ausgeheizt werden kann. Figur 1 zeigt dieses System schematisch 
dargestellt. Die im gestrichelten Rechteck 1 eingeschlossenen Elemente des Systems 
werden in einem Ofen bei 400-450°C ausgeheizt. Die einstufige Hg-Diffusions- 
pumpe 2 hat als Treibdiise eine Laval-Diise, und die Quecksilberverdampfung im 
Boiler 3 erfolgt in einer Lichtbogen-Entladung [1]1!"). Das Quecksilber wird der 
Pumpe iiber zwei barometrische Sadulen 4 und 5 zugefiihrt. Mit Hilfe der atmo- 
spharenseitig angebrachten Schwimmer 6 und 7 kann der Hg-Stand in diesen 
Saulen gehoben und gesenkt werden, womit wahrend des Ausheizprozesses das 
Quecksilber aus der Pumpe entfernt werden kann. Auch die Ziindung des Licht- 
bogens erfolgt mit Hilfe dieser Hebevorrichtung. Damit wahrend des Betriebes die 
atmospharischen Druckschwankungen das Brennen des Lichtbogens nicht beein- 
flussen, sind die Schwimmer in den hermetisch abschliessbaren Gefadssen unter- 
gebracht. Unterhalb des Bodens des Ausheizofens werden die Hg-Zufihrungslei- 
tungen mit Wasser gekiihlt. Bei allen UHV-Systemen mit Hg-Pumpen ist die Aus- 
_ frierung des Hg-Dampfes von grésster Bedeutung. Vor allem ist es wichtig, dass in 

den Kiihlfallen der Stand der fliissigen Luft konstant gehalten wird, hierfir sind 
verschiedene Anordnungen mit automatischer Nachfiillung der fliissigen Luft be- 
kannt geworden [2, 3, 4]. In unserem UHV-System haben wir zwischen der aus- 
heizbaren Hg-Pumpe und dem Rezipienten zwei Kihlfallen 8 und 9 vorgesehen. 
Fiir die dem Rezipienten nachstgelegene Kiihlfalle 9 haben wir in Abwandlung 
einer friiheren Konstruktion von VENEMa [5] eine neue Kiihlfalle entwickelt, deren 
Wirksamkeit fiir langere Zeit vom Absinken des Spiegels der fliissigen Luft unab- 
hangig ist. Wie aus Figur 1 ersichtlich, besteht diese Kiihlfalle aus vier ineinander 
geschachtelten Kugeln aus Pyrex-Glas. Der Durchmesser der dussersten Kugel 
betragt 165 mm, und die kleinste Kugel weist einen Durchmesser von 65 mm auf. 


16) Abteilung fiir industrielle Forschung, ETH. ye: 
17) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 542. 
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Figur 1 


Schematische Darstellung des UHV-Systems mit ausheizbarer Hg-Diffusionspumpe: 1 Ausheiz- 
ofen, 2 ausheizbare Hg-Diffusionspumpe, 3 Boiler mit Lichtbogenheizung, 4 und 5 barometrische 
Saulen fir die Hg-Zufiihrung, 6 und 7 Schwimmer zum Heben und Senter des Quecksilbers 
8 Vorkiihlfalle, 9 spezielle Kugelkiihlfalle, 10 Rezipient, 11 Bayard-Alpert-lonisationsmianemer a 


12 Hg-Sperre, 13 zu den Hg-Diffusions-Vorpumpen. 
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Die beiden dusseren Kugeln sind in einem zylindrischen Hals vereinigt. Wenn die 
Kihlfalle mit fliissiger Luft vollgefiillt wird, so dauert es etwa 1 See bis der 
Spiegel der fliissigen Luft das innere_ Kugelpaar erreicht. Wahrend dieser Zeit 
bleibt die Wiederverdampfung des Quecksilbers von den erwahnten Partien der 
Kiihlfalle fiir den Rezipienten ohne Bedeutung. Auch die weitere Verfliichtigung 
der fliissigen Luft ist wenig wirksam, da sich in der mittleren Kugel infolge sehr 
guter Warmeisolation die fliissige Luft sehr lange halt. 


p 
; [10-" Torr 


CES UTA ee gee | if 6 ee Gh MD 
—~— t[{min} 
Figur 2 
Einfluss des Standes der fltissigen Luft in den Kihlfallen auf das Endvakuum. 


Das Vorvakuum fiir das beschriebene System wird von zwei in Serie geschal- . 
teten einstufigen Hg-Diffusionspumpen erzeugt. Diesen ist wiederum iiber eine mit 
Trockeneis gekiihlte Falle die Rotationspumpe vorgeschaltet. Das UHV-System 
kann mit Hilfe eines Hg-Ventiles von den Vorpumpen getrennt werden. Nach 
griindlichem Ausheizen des UHV-Systems (zum Beispiel eine Ausheizung tber 
Nacht von 14 Stunden bei etwa 430°C) konnten zurzeit wiederholt Endvakua von 
ca. 10-11 Torr. erreicht werden. Fiir die Vakuummessung wurde ein verbessertes 
_Bayard-Alpert-Ionisationsmanometer verwendet. Die Ro6ntgenstrom-Grenze des 
neuen Manometers betragt 1 x 107! Torr. bei einer Anodenspannung von 150 V 
und 2,5 x 10-42 Torr. bei einer Anodenspannung von 50 V. Diese Verbesserung 
wurde durch die Verdoppelung der Manometerempfindlichkeit (C = 26 Torres an: 
= C-7_p) und durch Verkleinerung des Kollektordurchmessers auf 25  erreicht. 
Durch eine neuartige Halterung des Kollektorfadens , in dem das freie Fadenende 
in einer Kapillare gefiihrt wird, konnte diese Herabsetzung des Kollektordurchmes- 
sers auf 25 w betriebssicher realisiert werden. Es scheint, dass auf diese Art auch 
eine weitere Verkleinerung des Kollektordurchmessers noch méglich ist. 

Figur 2 zeigt uns den Einfluss des Standes der fliissigen Luft in den Kiihlfallen 


auf das Endvakuum. 
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Die Zustellung der Manuskripte an die Redaktion ZAMP von Seiten der 
Schweizerischen Physikalischen Gesellschaft erfolgte erst am 15. Juli 1960. Es war 
daher nicht méglich, die Berichte in einem der vorgesehenen Hefte (ZAMP 4 oder 
5/1960) zu ver6ffentlichen. 


International Federation of Information Processing Societies (IFIPS) 
Internationale Féderation fiir Datenverarbeitung 


Der Vorstand der Internationalen Féderation fiir Datenverarbeitung hat am 
16.-17. Juni 1960 seine erste Sitzung im Provisorischen Internationalen Rechen- 
zentrum in Rom abgehalten. 

Diese Féderation wurde anlasslich der durch die UNESCO abgehaltenen Ersten © 
Internationalen Konferenz fiir Informationsverarbeitung (ICIP), die im Juni 
1959 in Paris stattfand, gegriindet. Mehr als 1800 Delegierte aus allen Teilen der 
Welt haben dieser Konferenz beigewohnt, an welcher erfolgreiche Beziehungen ~ 
zwischen den Wissenschaftern auf diesem wichtigen Gebiet angebahnt wurden. — 

Der Vorstand setzt sich wie folgt zusammen: Isaac L. AUERBACH, Prasident, 
Auerbach Electronics Corp. National Joint Computer Committee, Philadelphia 3 — 
Pa. (USA), Prof. Dr. A. WALTHER, Vize-Prasident, Deutsche Arbeitsgemeinschaft 
fiir Rechenanlagen (DARA), TH Darmstadt, und Dr. A.P.SpEIsER, Sekretar- — 
Kassier, IBM Forschungslaboratorium Adliswil-Ziirich, als Vertreter der Schwei- — 
zerischen Gesellschaft fiir Automatik. 

Die Zweite Internationale Konferenz fiir Datenverarbeitung soll im September 
1962 in Deutschland stattfinden. Prof. Dr. A. WALTHER wurde zum Vorsitzenden — 
der Konferenz ernannt und Dr. Nrers Ivar Becu (Danemark) ist Vorsitzender 
des Programm-Komitees. 

Ausser der Organisation internationaler Konferenzen besteht die Tatigkeit 
dieser Féderation in der Koordination von Bestrebungen zur Schaffung von ~ 
Normen auf dem Gebiet der Datenverarbeitung, wofiir eine Kommission ins Le- 
ben gerufen werden soll. Auf Anregung der nationalen technischen Gesellschaften — 
werden weitere Aktivitaten aufgenommen werden, insoweit sie der Weiterent- 
wicklung der Datenverarbeitung, eines der wichtigsten Gebiete der modernen — 
Technik, forderlich sind. 

Nationale Gesellschaften, die der Féderation angeschlossen sind: 
Deutschland Deutsche Arbeitsgemeinschaft fiir Rechenanlagen (DARA) 
Belgien Association Belge pour l’ Application des Méthodes 
Scientifiques de Gestion (inklusive Association Belge 
de Calcul Numérique) 

Kanada The Computing and Data Processing Society of Canada 
Danemark Akademiet for de Tekniske Videnskaber 


ee 
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Spanien Consejo Superior de Investigaciones Cientificas 

Finnland The Finnish National Committee for Information Processing 
Frankreich Association Frangaise de Calcul (AFCAL) 

Japan Information Processing Society of Japan 

Holland Nederlands Rekenmachine Genootschap 

England The British Computer Society 

Schweden Svenska Samfundet fér Informationsbehandling 

Schweiz Schweizerische Gesellschaft fiir Automatik 


Tschechoslovakei Commission for Technical Cybernetics of the 
Czechoslovak Academy of Sciences 

Russland Computing Centre of the USSR Academy of Sciences 

S.A. National Joint Computer Committee 


A. P. SPEISER 


Seventh National Symposium on Reliability and Quality Control 

in Philadelphia, Pa. USA. 

In Philadelphia, Pa. (USA) findet vom 9. bis 11. Januar 1961 das 7th National 
Symposium on Reliability and Quality Control statt. Betriebssicherheit und Quali- 
tatsiiberwachung sind Gesichtspunkte, die in heutigen elektronischen Gerdten, 
welche tausende und oft zehntausende von Schaltelementen enthalten, grésste 
Bedeutung erlangt haben. Der Zweck dieser Tagung, welche vom Institute of 
- Radio Engineers (IRE) und anderen wissenschaftlichen Vereinigungen durchge- 
fiihrt wird, ist, Spezialisten, welche auf dem Gebiet der Betriebssicherheit und 
Qualitatsiiberwachung iiber Erfahrung verfiigen, zusammenzufiihren. Programme 
und Anmeldungsformulare sind zu beziehen bei J. H. Goopman, Publicity Chair- 


man, Burroughs Corporation, Research Center, Great Valley Laboratory, Paolt, 
Penna., USA. A. P. SPEISER 


Errata 


On the Scattering of Waves by a Disk. By Avsert E. Heins and RICHARD 
C. MacCamy. ZAMP 77, 249-264 (1960). 


P,. 249 

| Equation (1.2) should read 
Uty,. ¢;.0) = 0. fone < 1. 
The definition of R? in the line below Equation (1.3) should read 
R2= 72+ 9?— 2recos(@— 4) + 2. 

mP. 251 
Line 2 from the bottom: ‘(1.4)’ should read ‘(1.3)’. 
mE. 252 
The operator N,,(v, z, g) is identically zero and should be deleted. 


In Equation (I) replace (—)” by (—)"/2. 
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PeZ53 | 
Line 8 from the bottom: Replace (—)"/4 by (—)”/2. i 


PP, 254 
Equation (3.1): Replace u(r, z) by a w(r, 2). 


IPA255 
Delete (—)” in Equation (4.1). 
Delete (—)” in Equation (4. 3). 


Replace U by u in Equation (4.1). 


ees oe ae eS 


IPF 256; 
Delete (—)" in the definition of H (1 — y). 
The upper limit of the second integral in the definition of H (1 — y) should be y. — 


For ‘H (1 — y)’ read ‘H” (1 —y)’. 

P. 257 
Line 14 from the top should read + 

1 é be 

ieee (RK, on eee F 

Ey m\() =e ae i VB ae : 

iB Ate : 
Line 11 from the bottom: Replace (5.4) by (5.5). 

f 


Pe Z62 


Equation (6.2): In the second integral ‘cosy’ should read ‘cosh’. 


P2603 
Equation (6.4) should read 


t oe 1 
RV — y2 i 20 72 
ela aa c*(t) dt +i / om Bye c2(t) dt 


Ve — y2 Vy? — pf 
aa ‘ y? coshk, y 3 y sinhk 3 coshk 
= 48 sinty | 2 zo 7 eV 4 ha U4 AE+ Aby?. 


Line 7 from the bottom: Replace ‘ (6.2) and (6.3)’ by ‘(6.3) and (6. 4)’. 


Line 5 from the bottom: The sentence which begins ‘These are evaluated ..’ 
should read ‘These are evaluated by the requirement that c(t) and c?(¢) are finite 
at t= 0 and obey the regularity conditions which we described earlier. 


ee 
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Addendum 


Ae ri SSDS fondamentale d’une membrane vibrante: évaluations 
par défaut et principe de maximum. Par JosEpH Hrerscu, ZAM 
387-413 (1960): : ae Tee 


Ad § 4.4 et § 5: Ce résultat central se trouve essentiellement déja dans 
E. Picarp: Traité d’ Analyse, t. II (1te éd.: 1893), p. 25-26, et dans T. Boacero: 
Sull equazione del moto vibratorio delle membrane elastiche, Atti Accad. Lincei, ser 5 
vol. 16 (2° sem.), p. 386-393 (1907), notamment p. 390. — Leurs hypothéses de 
continuité sont les mémes que celles de M. H. PRotrer, que j’ai critiquées au 


§ 5. — Contrairement & Protrer, Picarp et Boaceto semblent avoir sous-estimé 
importance de cette inégalité : elle apparait dans la démonstration de simples 
propriétés de monotonie. J. HeRscH 
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Grundziige der Elektroakustik. Von F. A. FiscHEer (Fachverlag Schiele & 
Schén GmbH, Berlin 1959). 210 S., 141 Abb., 7 Tabellen; DM 24.-. 

Nach genau 10 Jahren kam diese 2. Auflage heraus. Sie wurde in einigen 
Punkten wesentlich gedndert. Wertvoll ist, dass auch auf der mechanischen Seite 
die mechanische Innenimpedanz dem Wandler als Vierpol zugeordnet bleibt. Dieser 
ist dann definiert als ein Vierpol mit elektrischen Eingangsklemmen und einer 
Ebene, in welcher der Zusammenhang zwischen Kraft und Schnelle aus der gege- 
benen Belastung hervorgeht. 

Das leichtverstandlich geschriebene Buch gliedert sich in die Kapitel elektrische 
und mechanische Schwingungsgebilde, in welchem die Analogien behandelt werden, 
dann die Kraftwirkungen elektrischer und mechanischer Felder auf Materie, auf- 
geteilt in elektrische und magnetische Felder. Es folgen in Kapitel 3 die elektro- 
- mechanischen Umwandlungsprinzipien und ihre Grundgesetze, ebenfalls wieder 
aufgeteilt nach Wandlern mit elektrischen Feldern und solchen mit magnetischen 
Feldern. Darin diirfte sich ein Irrtum eingeschlichen haben, weil die normalen 
elektromagnetischen Telephone nicht zu den behandelten elektromagnetischen 
Wandlern geh6éren, sondern zum elektrodynamischen Wandler. In einer Tabelle des 
Kapitels 4 ist die Systematik der elektroakustischen Wandler zusammengestellt. 
Es folgt in Kapitel 5 die Grundgleichung der Ersatzschaltbilder der vier Wandler- 
klassen. Das Kapitel 7 ist der Abstrahlung des Schalles gewidmet, wobei die ge- 
richtete Strahlung eine wesentliche Rolle spielt. Wichtig sind die Fragen des 
Wirkungsgrades eines Wandlers als Schallsender (Kapitel 8), die Empfindlichkeit 
eines Wandlers als Schallempfanger (Kapitel 11) und die elektrischen und mecha- 
nischen Eingangsimpedanzen (Kapitel 9 und 10). Sehr kurz ist der thermische 
- Schallwandler und der Relais-Empfanger gestreift. Den Schluss bildet ein Kapitel 
iiber die elektroakustische Messtechnik. 

Auf beschranktem Raum wird das Wesentliche gesagt und griindlich behandelt. 
Der Elektroakustiker wird deshalb gerne nach diesem Buch greifen, das ihm hilft, 
seine Probleme richtig anzupacken. Sollte er keinen geniigenden Aufschluss im 
- Buche finden, so weist ihn ein Schrifttumsverzeichnis auf eine ausfiihrlichere Be- 


~ handlung hin. H. WEBER 


- ZAMP X1/35 
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IS-Tafel fiir Luft und Verbrennungsgase. Von O. Lutz und F. WoLrF. 
2. Auflage, neu bearbeitet und erweitert von Prof. Dr.-Ing. O. Lutz und Dipl.-Ing. ; 
W. Louse (Springer-Verlag, Berlin 1959). 18 S., 4 Abb., 5 Kurventafeln; DM 7.50. 

Die IS-Tafel fiir Verbrennungsgase von Lutz und Wo r, die bereits jahrzehnte- 
lang vielen Ingenieuren wertvolle Dienste geleistet hat, liegt nun in 2. Auflage vor. 
Wie der Verfasser im Vorwort sagt, hat er sich zundchst die Frage gestellt, ob die 
Herausgabe dieser Tafel heute tiberhaupt noch vertretbar sei, oder ob man nicht 
vielmehr in Zukunft auf das graphische Arbeiten anhand solcher Tafeln vdllig ver- 
zichten werde. Man muss ihm aber sicher zustimmen, wenn er zum Schluss gelangt, 
dass auch weiterhin die IS-Tafel ein unentbehrliches Hilfsmittel des Ingenieurs 
bleiben wird. Der Ingenieur sollte sich nicht zu sehr daran gewohnen, sogleich zu 
komplizierten automatischen Rechengeraten seine Zuflucht zu nehmen, sondern 


er sollte durchaus versuchen, mit einfachen anschaulichen Verfahren rasch zum ~ 


Ziel zu gelangen, dort wo der Charakter des Problems dies zulasst. : 
Die neue Tafel ist gegeniiber der alten ins Gebiet wesentlich héherer Tempera- 
turen erweitert, was einen zusdtzlichen Kunstgriff zur Beriicksichtigung des dann 


wesentlichen Einflusses der Dissoziation notwendig gemacht hat. Sonst ist die — 
Tafel in ihrem Aufbau gleich geblieben wie die alte. Im beigefiigten Text sind die — 


Grundlagen klar und knapp zusammengefasst, und es ist auch angegeben, wie im © 


Falle abnormaler Kraftstoffe (die der Verfasser «exotische» Kraftstoffe nennt) zu — 


verfahren ist. 


Die neue IS-Tafel wird sicher wie die alte eine weite Verbreitung finden. An ~ 
sich méchten wir bedauern, dass man bei dieser Uberarbeitung noch beim alten — 


technischen Ma8system geblieben ist und nicht den Ubergang zum MKSA-System — 


vollzogen hat. Man versteht aber leicht, dass dies im heutigen Zeitpunkt in bezug — 


auf die Verbreitung der Tafel vielleicht doch noch etwas gewagt gewesen ware. 


W. TRAUPEL ~ 


Méthodes Numériques; Interpolation, Dérivées. Von J. KUNTZMANN © 


(Dunod Editeur, Paris 1959), 272 S., 60 Fig.; fFr. 3600.-. 


Dieses Buch behandelt in ausfiihrlicher und leichtverstandlicher Darstellung — 
die verschiedenen Methoden der Interpolation und der numerischen Berechnung © 
von Ableitungen. Die stoffliche Beschrankung gestattet es dem Verfasser, die Teil- — 


gebiete einschliesslich deren Querverbindungen sehr griindlich zu behandeln. Sorg- 
faltige Fehlerabschatzungen und viele numerische Beispiele zu allen behandelten 
Methoden machen das Buch zu einer wertvollen Hilfe fiir praktische Anwendungen. 


M. ENGELI 


_Advanced Mechanics of Fluids. Herausgegeben von HuNTER RousE (John 
Wiley & Sons, Inc., New York 1959). 144 S., 142 Fig.; $9.75. 


EE 


Diese Fortsetzung der «Elementary Mechanics of Fluid» des als Hydrauliker — 


bekannten Autors beschrankt sich auf die Behandlung der inkompressiblen Fliissig- 
keiten. Sie ist vor allem als Hilfe fiir Vorlesungen iiber héhere Strémungslehre 
gedacht, und dieser Zweck ist zweifellos gut erfiillt. Nach einer interessanten Dis- 


kussion iiber die Wechselwirkung zwischen Theorie und Experiment wird die Me- 


thode der Dimensionsanalyse behandelt, und dann werden die klassischen Theorien 
der Potentialstrémung und ihre mathematischen 
verfahren und konforme Abbildung, 
iiber Navier-Stokessche Gleichungen, 


eine kurze Diskussion des Turbulenz i i i 

) problems. Die Theorie der laminaren und 
a Grenzschicht und der Mischvorgange bei freier Turbulenz wird aus- 
uhrlich behandelt. Dem Zweck des Buches entsprechend wird das Hauptgewicht 


ausfiihrlich dargestellt. Es folgt ein Kapitel 


Hilfsmittel, wie Singularitaten- — 


laminare Strémung und ihre Stabilitat und 
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auf die mathematische Behandlung der Probleme gelegt; es sind nur wenige, zur 
Hauptsache altbekannte experimentelle Resultate zur Illustration der Theorie 


wiedergegeben. Viele gut ausgelesene Probleme und numerische Beispiele erganzen 
das sehr sch6n gedruckte Buch. P. DE HALLER 


Glossary of Terms Used in Vacuum Technology. Herausgegeben vom 
Committee on Standards der American Vacuum Society (Pergamon Press, London 
1958). 63 S.; $2.00. 

Durch die rasche Entwicklung der Vakuumtechnik entstand in vielen Landern 
das Bediirfnis, eine auf diesem Fachgebiet allgemeingiiltige Nomenklatur einzu- 
fiihren. Das vorliegende Heftchen versucht, diesem Bediirfnis im angelsachsischen 
Sprachgebiet entgegenzukommen. Insgesamt werden 428 vakuumtechnologische 
Begriffe definiert und naher umschrieben. Man findet unter anderem den begriissens- 
werten Vorschlag, auch in Amerika die Druckeinheit 1 mm Hg-Saule mit «Torr» zu 
bezeichnen und damit fiir eine internationale Verwendung dieses Begriffes den Weg 
zu ebnen. E. Bas 


Les problémes aux limites de la physique mathématique (Introduction 
& leur étude générale). Von H. G. Garnir. Lehrbiicher und Monographien aus dem 
Gebiete der exakten Wissenschaften, Mathematische Reihe (Birkhauser Verlag, 
Basel/Stuttgart 1958). 234 S., 5 Abb.; sFr. 25.-. 

Die vorliegende Monographie behandelt nicht etwa, wie es der Titel angibt, 
simtliche Randwertprobleme der mathematischen Physik, sondern nur das «Dirich- 
let-Neumannsche Problem» fiir die Wellen- und Warmeleitungsgleichung. Der Ver- 
fasser hat sich zum Ziel gesetzt, den Leser in die moderne funktionalanalytische Be- 


_handlungsweise von Randwertproblemen einzufiihren. Es handelt sich hier also um 
ein Werk, das sich an diejenigen Mathematiker wendet, welche die klassische 


Theorie dieser Randwertprobleme schon kennen. Dementsprechend fehlen Moti- 
vierungen und Beispiele vollstandig. Der Einsatz der Funktionalanalysis ergibt als 
Hauptresultat die Méglichkeit, die Randwertprobleme ohne irgendeine Annahme 
iiber den Rand der Gebiete, in denen das Problem gelést werden soll, zu formulieren 
und zu «lésen». Insbesondere ist es méglich, die Randbedingungen, welche die 
Normalableitung in den Randpunkten betreffen, auch in diesem ganz allgemeinen 
Fall zu formulieren. Die «Lésungen» sind dann allerdings auch keine Funktionen, . 
sondern Distributionen. Diesen Distributionen kann nicht ohne weiteres ein phy- 
sikalischer Sinn zugesprochen werden. Der Autor unterscheidet deshalb zwischen 
«mathematischen» und «physikalischen» Lésungen des Problems. «Physikalisch» 
heisst dabei eine Lésung, wenn die entsprechende Distribution Integralform hat. 
Das Buch setzt eine sehr griindliche Kenntnis der Lebesgueschen Integrations- 


_ theorie voraus, hingegen wird die Theorie der Hilbertraume sowie der Distributionen 


soweit wie nétig entwickelt. Etwas seltsam mutet die ausfiihrliche Behandlung 
elementarer Begriffe der linearen Algebra (Skalarprodukte usw.) an; wer die 
Lebesguesche Theorie kennt, wird auch lineare Algebra kennen. Wenn man sich 
aber schon entschliesst, derart elementare Dinge ausfiihrlich zu behandeln, so 


-sollten auch an anderen Orten Motivierungen gegeben werden, so fiir den Begriff 


_ondes et de la diffusion. 


der schwachen Lésung usw. Die Ausstattung des Buches entspricht derjenigen der 

iibrigen Publikationen dieses Verlags. 
Inhaltsiibersicht: Livre 1: Théorie des espaces fonctionnels hilbertiens; Livre ie 

Probléme de DirICHLET-NEUMANN pour l’opérateur méta-harmonique; Livre III: 


L’opération J; Livre IV: Probléme de Dir1IcHLET-NEUMANN pour les opérateurs des 
E. STAMM 


; i 
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Handbuch der Physik - Encyclopedia of Physics. Herausgegeben von — 
S. FLricce und E. Creutz, Band 41/1: Kernreaktionen 11: Theorie (Springer- § 
Verlag, Berlin 1959). 571 S., 61 Fig.; DM 145. i 

tiber den Wert bzw. Unwert von Enzyklopddien ist an anderer Stelle schon so ~ 
viel geschrieben worden, dass es nicht verlohnt, diese Frage erneut aufzugreifen. & 
Dagegen drangt sich oft der Vergleich zwischen Handbiichern und riickblickenden ~ 
Zeitschriften, wie etwa der «Review of Modern Physics», auf. Solche Zeitschriften — 
haben gegeniiber Handbiichern zweifellos Vorteile: Sie verzichten darauf, Artikel — 
in einem gewissen Zusammenhang zu ordnen, das heisst, die jeweiligen Aufsatze — 
erscheinen zu einem Zeitpunkt, der ihnen durch den momentanen Entwicklungs- 
stand der Physik zugewiesen ist. Damit bilden sie im besten Sinne des Wortes ein — 
Kontinuum, denn es ist ohne weiteres méglich, ein bereits behandeltes Gebiet zu ~ 
einer spateren Zeit wieder aufzugreifen, um es den neusten Erkenntnissen gemass 
zu erganzen und zu berichtigen. Auf der andern Seite gelingt es in guten Hand- ~ 
biichern, dem Leser das momentane Wissen um ein bestimmtes Sujet in tibersicht- — 
licher und konzentrierter Form vorzulegen. 

In dem vorliegenden Band wird das eben geschilderte Ziel besonders durch das ~ 
erste Kapitel tiber die Theorie von Resonanzreaktionen (G. Breit, 407 S.) in héchst ~ 
vollkommener Weise erreicht. Allgemeine Uberlegungen iiber Streuungen an Zen- 
tralfeldern fiihren den Leser in den Wignerschen -Matrix-Formalismus ein, dem 
ein besonders breiter Raum gewidmet wird. Das Kapitel klingt aus mit Hinweisen ~ 
auf die Anwendung der Theorie auf spezielle Probleme, wie Reaktionen nahe der — 
Schwelle, Riesenresonanzen bei Kernphotoreaktionen, durch schwere Ionen aus- — 
geloste Prozesse usw. Das Kapitel von H. Hutz und G. Breir (57 S.) tiber Coulomb- © 
Wellenfunktionen verzichtet auf Hinweise iiber Anwendbarkeit dieser Funktionen, — 
sondern bietet dem Leser ausschliesslich eine recht konzentrierte und tibersichtliche ~ 
Darstellung der mathematischen Handhabung dieser wichtigen Funktionen (Ap- ; 
proximationen). In einem weiteren Kapitel diskutieren G. Breit und I. S. McIn-— 
TOSH (29 S.) die Polarisation von Nukleonen durch Streuung an Atomkernen, ~ 
wahrend sich das letzte Kapitel (G. BReir und R. L. GLucksTERN, 64 S.) mit der — 
Coulomb-Anregung beschaftigt. Dieser Artikel diirfte als Erganzung des bekannten 3 
Aufsatzes tiber Coulomb-Anregung von ALDER ét. al. [Rev. Mod. Phys. 28, 432 © 
(1956)] zu verstehen sein. 

Samtliche Artikel erscheinen in englischer Sprache; das Inhaltsverzeichnis wird ; 


ee tee 
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deutsch und englisch gefiihrt. Das Buch wendet sich in erster Linie an den theo- — 


retisch interessierten Experimentalphysiker, dem es zweifellos von erheblichem — 
Nutzen sein wird. F. HEINRICH 


E 
i 

Plasma Physics - Accelerators - Thermonuclear Research (Journal of | 
Nuclear Physics: Part C.) Herausgegeben von J. V. DuNworTH (Pergamon Press, 
London 1959). Preis: $30.— jahrlich. 

Zur kurzen Charakterisierung der Aufgaben, die sich diese neue Zeitschrift stellt, 
diirfte die Angabe einiger Beispiele aus dem Inhalt des ersten Heftes am geeignet- 
sten sein: Das Plasmabetatron, Gleichgewichtskonfigurationen toroidaler Plasmen, 
Zirkularbeschleuniger mit azimutaler Magnetfeldvariation, ‘Colliding beam’-Sy- 
steme mit hohen Protonintensitaten, Ionenquellen fiir Tandem-Bandgeneratoren 
allgemeine Charakteristika von selbsterhaltenden Fusionssystemen und vieles mehr. 
Besonders zu begriissen ist es, dass dem auf diesen Gebieten interessierten Leser 
auch die Ubersetzung russischer Artikel geboten wird. Ferner enthalt jede Nummer 
ein recht volistandiges Verzeichnis kiirzlicher Publikationen aus etwa 50 bis 100 


anderen Zeitschriften iiber Themen der Plasma- und Fusionsphysik und iiber den 


Acceleratorbau. F. HEINRICH 


Lehrbiicher und Monographien aus dem 
Gebiete der exakten Wissenschaften 


Mathematische Reihe 


Band 1: Versicherungsmathematik. Von E. Zwinggi, 2. Auflage (1958). Fr./DM 36.-. 
Band 2: Analytische Geometrie der Ebene und des Raumes. Von R. Fueter (1945). Fr./DM 24.-. 


Band 3: Statistische Methoden fiir Naturwissenschafter, Mediziner und Ingenieure. Von A. Linder, 3 
umgearbeitete und stark erweiterte Auflage (1961). Fr./DM 54.-. 


Band Vorlesungen tiber Differential- und Integralrechnung. Von A. Ostrowski. Band I. 2., neubearbei- 

4/5/7: tete Auflage (1960): Funktionen einer Variablen. Fr./DM 35.—. Band II: Differentialrechnung 
auf dem Gebiete mehrerer Variablen. 2., neubearbeitete Auflage (1961). Ca. Fr./DM 40.-. 
Band III (1954): Integralrechnung auf dem Gebiete mehrerer Variablen. Fr./DM 75.—. 


Band 6: Lehrbuch der darstellenden Geometrie. Von E. Stiefel (1947). Fr./DM 28.50. 


Band8,9: Funktionentheorie. In zwei Banden. Von C. Carathéodory. Band I (1950): Fr./DM 39.50, 
Band II (1950): Fr./DM 28.-. 


Band 10: Ma8 und Integral und ihre Algebraisierung. Von C. Carathéodory (1956). Fr./DM 39.-. 
’ Band 11: Uber Kurven und Flichen in allgemeinen Raumen. Von P. Finsler (1951). Fr./DM 16.-. 


Band 12: ag een Laplace-Transformation. Von D. Voelker und G. Doetsch (1950). 
Yr. 48.-. 


Band 13: Theorie der geometrischen Konstruktionen. Von L. Bieberbach (1952). Fr./DM 19.—. 


Band Handbuch der Laplace-Transformation. Von G. Doetsch. Band I (1950): Fr./DM 84.~. Band II 
14/15/19: (1955): Fr./DM 57.-. Band III (1956): Fr./DM 42.-. 


Band 16: Analytische Geometrie. Von W. Blaschke (1954). Fr./DM 19.-. 
Band 17: Projektive Geometrie. Von W, Blaschke (1954). Fr./DM 19.-. 
Band 18: Ubungen zur projektiven Geometrie. Von H. Herrmann (1952). Fr./DM 17.-. 
_ Band 20: Elementare Dijferentialgeometrie. Yon W. Haack (1955). Fr./DM 22.-. P 
Band 22: Die Theorie der Gruppen von endlicher Ordnung. Von A. Speiser (1956). 4. Auflage. Fr./DM 26.—. 
Band 23: Les problémes aux limites de la physique mathématique. Par H. G. Garnir (1958). Fr./DM 29.-. 


Band 24: Einfiihrung in Theorie und Anwendung der Laplace-Transformation. Von G. Doetsch (1958). 
Fr./DM 39.40. 


Band 25: Vorlesungen iiber Funktionalgleichungen und ihre Anwendungen. Von J. Aczél. Erscheint 1961. 
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